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ABSTRACT 


This  report  presents  an  analysis  of  the  dynamics  of  supercritical  shafts  on 
many  flexible  supports  in  terms  of  a  so-called  "transmission  line 
analogy".  The  primary  objective  of  the  study  is  to  develop  a  direct 
analytical  approach  for  optimizing  the  support  conditions,  in  terms  of 
minimum  flexural  vibration  behavior,  for  supercritical  shafts  flexibly 
supported  on  mass- spring-damper  units  at  different  locations  along  the 
shaft.  The  approach  is  based  on  the  traveling  wave  concept  as  used  in 
electrical  transmission  line  theory.  The  governing  differential  equation 
used  in  this  analogy  includes  terms  which  account  for  the  effects  of 
rotating  inertia,  gyroscopic  motion,  and  shear  deformations. 

If  the  solution  of  the  governing  differential  equation  is  manipulated  by 
means  of  rather  complicated  matrix  algebra,  the  dynamic  response  of 
the  rotating  shaft  can  be  expressed  in  traveling  wave  form,  which  leads 
to  the  criterion  for  impedance  matching  and  optimized  support  conditions. 
The  impedance  at  each  supporting  location  corresponding  to  minimum 
vibration  response  must  equal  the  characteristic  impedance  of  the  shaft. 
This  condition  is  termed  a  "matched"  condition. 

A  weaker  than  optimum  form  of  impedance  matching  ir  the  "quasi- 
matched"  condition,  in  wh’ch  only  the  predominant  term  of  the  reflection 
matrix  for  a  support  is  made  to  vanish. 

A  rotating  shaft  with  three  supports  is  vised  to  illustrate  the  matched 
impedance  concept  for  determining  optimum  support  conditions  for 
multisupported  hypercritical  shafts. 

The  study  has  led  to  the  following  specific  conclusions: 

1.  The  transmission  line  analogy  can  be  extended  to  shafts  having 
any  number  of  interior  supports.  However,  the  solution  is 
considerably  more  complicated  for  shafts  having  more  than 
one  interior  support. 

2.  For  the  shaft  with  both  end  impedances  matched  with  the 
characteristic  impedance  of  the  shaft,  no  intermediate 
support  is  needed  to  assist  in  the  minimization  of  vibration 
response. 

3.  If  it  is  physically  impossible  or  impractical  to  terminate  a 
shaft  in  its  characteristic  impedance,  quasi- matched  end 
impedances,  or  quasi- matched  interior  supports  when  end 
conditions  are  not  available  for  optimization,  should  provide 
good  performance. 

4.  When  the  shaft  and  end  support  impedances  are  not  matched 
and  only  one  interior  support  is  used,  two  approaches  may 

be  employed  to  assist  in  the  minimization  of  vibration  response: 
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a.  If  both  ends  have  the  same  configuration,  the  matched 
interior  support  may  be  placed  at  the  mid- span. 

b.  If  one  end  support  is  different  from  the  other,  the  closer 
the  matched  intermediate  support  is  placed  to  one  of  the 
ends,  the  more  effectively  it  will  minimize  the  vibration 
response  of  the  shaft. 

5.  When  the  shaft  and  end  support  impedances  are  not  matched, 
the  use  of  two  matched  intermediate  supports  placed  closely 
to  the  ends  of  the  shaft  is  recommended. 
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CHAPTER  1 


OBJECTIVES 


INTRODUCTION 


For  any  rotating  shaft,  there  exists  a  series  of  discrete  speeds  at  which 
centrifugal  forces  resulting  from  mass  unbalances  cause  progressively 
greater  shaft  deflections.  The  elastic  restoring  forces  developed  as  the 
shaft  deflects  are  overcome  by  these  ever-increasing  centrifugal  forces. 
Extremely  large  deflections  and  even  destruction  of  the  shaft  and  its 
bearings  can  result  from  operation  at  these  speeds,  called  critical  speeds. 
For  this  reason,  designers  of  power-transmission  equipment  normally 
avoid  the  problem  by  operating  shafts  below  their  first  critical  speed. 

There  are,  of  course,  disadvantages  to  restricting  operation  to  below  the 
first  critical  speed.  For  transmitting  a  given  horsepower,  torque  and, 
consequently,  shaft  size  must  be  increased  as  operating  speed  is  reduced. 
In  the  case  of  long  shafts,  the  shaft  size  must  be  increased  above  the  size 
required  to  transmit  the  torque  simply  to  raise  the  first  critical  speed 
above  the  operating  speed  range;  alternatively,  the  shaft  size  may  be 
determined  by  the  torque  loading,  but  additional  bearings  must  be  installed 
to  support  the  shaft  and  thereby  to  raise  its  first  critical  speed.  The 
major  disadvantage  of  these  conventional  practices,  especially  as  applied 
to  aircraft,  is  the  weight  penalty. 

It  has  been  shown  that  shafting  can  be  operated  consistently  far  above  its 
first  critical  speed,  with  consequent  savings  in  shaft  and  support  weight 
(references  5,  11,  and  21).  In  short,  supercritical- speed  shafting,  with 
its  associated  advantage  in  weight,  is  a  very  practical  and  feasible  means 
of  transmitting  power.  With  today1  s  high-speed  power  sources,  it  is 
especially  attractive,  sin^e  considerable  weight  could  be  pared  from 
engines  and  bearings  by  transmitting  power  at  the  same  speed  as  it  is 
produced. 

Transmitting  large  horsepowers  with  small- diameter  shafts  presents 
the  problem  of  controlling  shaft  vibration  at  the  critical  speeds.  Success¬ 
ful  operation  of  flexible  shafts  is  usually  achieved  by  balancing  to  reduce 
dynamic  forces  and  also  by  introducing  support  conditions  (spring, 
damping,  mass)  which  tend  to  minimize  runout  amplitudes  and/  or  bearing 
loads  at  the  important  critical  frequencies  (references  11  and  21).  Both 
of  these  techniques  should  be  employed  simultaneously  to  bring  about 
smooth  shaft  operation  through  the  critical  speeds.  Although  the 
provision  of  appropriate  spring,  damping,  and  mass  coefficients  (imped¬ 
ance)  at  the  supports,  alone,  may  permit  the  rotor  to  negotiate  the 
critical  speeds  in  an  acceptable  manner,  the  ease  with  which  this  may 
be  accomplished  will  be  greater  for  the  better  balanced  (less  crooked) 
shaft.  Balancing  of  a  supercritical  shaft  may  be  achieved  by  either  the 
proper  attachment  of  counterweights  or  the  placement  of  greater  restric¬ 
tions  on  the  fabrication  tolerances  of  the  shaft. 
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From  a  theoretical  poi  it  of  view,  different  mathematical  formulations 
(reference  8)  have  to  be  developed  to  implement  the  two  techniques.  The 
analysis  of  a  supercritical  shaft  supported  at  various  points  along  the 
length  of  the  shaft  by  flexible  damping  bearings  is  represented  with  good 
accuracy  by  the  steady-state  solution  of  the  equation  of  motion  for  the 
beam  vibrating  in  two  mutually  perpendicular  planes.  This  solution  has 
been  used  to  simulate  the  performance  of  the  rotating  shaft  under  actual 
running  conditions  in  which  the  shaft  defects  (unbalance,  and  initial 
crookedness)  provide  distributed  forcing  functions.  The  approach  provides 
a  means  to  study  the  effects  of  imperfections  and  balancing  on  shaft 
performance. 

On  the  other  hand,  the  solution  can  be  interpreted  also  in  terms  of  "a  trans¬ 
mission  line  analog".  This  approach  is  based  on  the  recognition  of  the 
existence  of  an  analog  between  the  amplitude  response  of  the  vibrating 
beam  and  the  voltage  amplitudes  in  electrical  transmission  lines.  The 
shaft  runout  is  treated  as  a  series  of  deflection  waves  (voltage  waves) 
traveling  along  the  shaft  (transmission  line).  These  waves  are  in  part 
absorbed  and  in  part  reflected  at  the  supports  (loads).  In  other  words, 
the  dynamic  responses  are  expressed  in  traveling  wave  form  along  the 
shaft  in  a  manner  analogous  to  the  treatment  of  electrical  response  waves 
in  transmission  line  theory  (references  3,  13,  and  16).  The  transmission 
line  analogy  solution  is  particularly  useful  for  the  direct  establishment  of 
the  support  conditions  needed  for  optimum  shaft  operation  through  the 
critical  speeds. 

The  purpose  of  this  report  is  to  present  the  development  of  a  general 
analysis  of  the  multisupported,  supercritical  shaft  in  terms  of  the  trans¬ 
mission  line  analogy  and  to  indicate  the  usefulness  of  these  tools.  The 
conditions  of  optimized  supports  corresponding  to  minimized  dynamic 
responses  will  be  established. 

SCOPE 


The  basic  mathematical  concepts  upon  which  the  transmission  line  analogy 
proposed  ir.  this  report  is  based  were  originally  developed  by  Nelson  in 
terms  of  a  shaft  on  end  supports  only  (reference  14).  Liu,  Friedericy, 
and  Eppink  extended  this  work  to  apply  to  shafts  having  one  additional 
intermediate  support  (reference  8).  The  studies  discussed  in  this  report 
are  concerned  with  an  extension  of  Nelson' s  electrical  transmission  line 
analogy  for  supercritical  shafts,  to  include  the  effects  of  any  number  of 
intermediate  supports.  The  derivation  of  the  equations  of  motion  of 
supercritical  shafts  with  respect  to  fixed  reference  axes  incorporates  the 
effects  of  rotational  inertia,  gyroscopics,  and  shear  deformations. 

Independent  of  the  University  of  Virginia  work,  Voorhees  and  coworkers 
(reference  5)  formulated  their  version  of  a  transmission  line  analogy.  In 
this  version  the  fourth-order  differential  equation  of  motion  for  beams  has 
been  reduced  to  an  approximating  second-order  equation  which  is 
completely  analogous  to  the  second-order  equation  which  governs  electrical 
network  behavior.  This  reduction  in  order  requires  that  a  one-to-one 
relationship  shall  exist  between  moments  and  deflections  of  the  shaft.  Such 
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a  relationship  can  only  be  brought  about  by  the  introduction  of  compromise 
boundary  conditions  at  the  supports.  However,  the  approach  has  the  ad¬ 
vantage  that  all  the  terminology  and  computational  aids  developed  for  the 
electrical  transmission  line  problem  can  be  utilized  directly  in  the  design 
of  supports  for  supercritical  shafts.  The  support  optimization  formulas 
in  the  Design  Manual  for  Supercritical- Speed  Power-Transmission  Shafts, 
prepared  by  the  Battelle  Memorial  Institute  for  the  U.  S.  Army  Transpor¬ 
tation  Research  Command  and  the  U.  S.  Air  Force  Research  and 
Technology  Division  (reference  2)  are  based  on  this  direct  approach,  and 
they  are  extremely  convenient  for  designing  supercritical,  power  trans¬ 
mission  shafting  because  of  their  simplicity  and  straightforwardness. 

In  the  transmission  analogy  of  this  report,  the  fourth-ordev  differential 
equation  which  governs  supercritical  shaft  behavior  is  solved  in  an  exact 
manner  and  the  various  component  terms  to  the  solution  are  worked  into 
standing  wave  forms  which  are  analogous  to  voltage  wave  forms.  The 
optimization  of  support  conditions  is  then  performed  in  the  same  manner 
as  loads  are  optimized  in  the  electrical  transmission  line  problem. 
Usually,  the  supports  are  formed  with  rotational  and  translational  mass- 
spring- damper  units  and  these  are  interpreted  as  impedances  to  the 
standing  waves. 

The  development  of  simple  formulas  for  the  determination  of  optimum 
support  values  has  been  attempted;  however,  due  to  the  complexity  of 
the  solution  of  the  multisupported  shaft,  a  shaft  with  only  one  interior 
support  has  been  used  as  a  specific  illustration.  The  various  formulations 
developed  in  the  illustration  should  be  useful  in  the  design  of  optimum 
support  conditions  of  supercritical  shafts  on  three  flexible  supports  and 
should  complement  the  equations  and  results  in  the  Design  Manual  men¬ 
tioned  earlier.  The  three  support  formulations  of  this  report  allow  for  the 
determination  of  optimum  interior  support  parameters  in  terms  of  end 
supports  which  are  not  just  fixed  and  simply  supported,  but  may  consist 
of  rotational  and/  or  translational  mass- spring- damper  units.  Specific 
results  can  be  found  in  Chapter  4  of  this  report. 

ASSUMPTIONS  AND  NOTATIONAL  CONVENTIONS 


It  is  shown  in  Figure  1  that  the  mathematical  model  of  a  pris¬ 
matic  shaft  is  embedded  in  a  right-hand  orthogonal  normalized 
space  coordinate  system  with  fixed  axes  Xj,  X^,  X^.  Let  Sq  represent 

the  shaft  with  elastic  curve  S.  Let  A^  and  Bq,  where  the  end  supports 

are  attached,  represent  the  end  bodies  with  elastic  curves  A  and  B, 
respectively.  A,  S,  and  B  coincide,  at  rest,  with  the  intervals  [a,  0], 

(0,  f  ),  and  [  J?n,  b],  respectively,  on  the  X^-axis,  where  a.<0<i^<b. 

A  closed  interval  and  an  open  interval  with  end  points  Ct ,  P  on  the  X^-axis 

are  represented  by  [ O' ,  P]  and  (ot  ,  P),  respectively.  Except  for  their 
common  connection  (the  shaft,  Sq),  Aq  and  Bq  are  dynamically  independent. 
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Figure  I.  Pictorial  Representation  of  A„,  S0,  B0  Embedded  in 
*1*  X2,  X,  With  a  Driving  Force  P*(aii.). 


The  intermediate  supports  are  attached  at  x  =  ,  j  =  1,  2,  .  .  .  ,  n-1, 

where  0  ^  .  Force,  torque,  and  motion  along  the 

Xj-axis  are  assumed  to  be  zero.  Constant  angular  velocity,  small 

transverse  motion,  axial  symmetry,  and  linearity  are  assumed  through¬ 
out. 


Where  C  denotes  complex  numbers,  let  denote  the  n- dimensional 

vector  space  over  a  field  of  complex  numbers  with  the  elements,  or 
vectors,  of  C  being  thought  of  as  "column  vectors".  In  the  following, 
the  appearance  of  a  bar  (possibly  together  with  other  symbols)  over  a 
quantity  indicates  that  it  is  an  element  of  C  for  some  n>l  (usually  n  =  2). 

Let  ASB  denote  collectively  the  bodies  A^,  Sq,  Bq.  Let  d(x)  be  that 

portion  of  ASB  which,  at  rest,  has  Xj  -  coordinate  x  on  [a,  b]  . 

A  "mass  element  on  [a,  b]  "  is  that  portion  of  ASB  which  lies  between 
d(x  -  dx /Z)  and  d(x  +  dx/  2)  where  dx  i6  small.  The  normal  to  a 
differential  mass  element  is  the  vector  normal  to  d(x)  at  its  intersection 
with  the  elastic  curve,  with  positive  orientation  being  the  same  ar*  that 
of  the  X^  -axis.  (See  Figures  2  and  3. ) 

For  pc  ints  on  elastic  curves  A,  S,  B,  "position"  as  the  vector  is  defined 
as  follows: 


Y*  = 

i 

*-< 
•-* 
_ i 

— 

'  y*  +  i  Y* 

1  11  12 

r - 

i _ 

Y*  +  i  Y* 

1 21  1  1 22 

This  is  a  continuous  function  of  a 

and  t 

to  C i.  e.  ,  2- 

vector  space  over  a  field  of  complex  numbers,  where  x  is  on  [a,  b]  . 

"*  "  denotes  the  unnormalized  variables.  YTl  andY12  are  the  pro¬ 
jections  of  the  deflections  at  a  general  point  of  the  elastic  curves  A,  S,  B, 
on  the  XjX^— ,  X^X^-planes,  respectively.  The  normal  to  a  differential 

mass  element  has  projections  on  the  planes  XjX^,  XjX^,  which  are 
inclined  at  acute  angles  with  respect  to  the  Xj-axis.  These  angles  of 
inclination  are  Y^j  ,  Y*^.  respectively,  where  the  positive  senses  of 
rotation  are  from  X^,  X^,  respectively,  to  Xj.  (See  Figures  2  and  3.) 
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Similarly,  "force"  as  a  vector  is  defined  as  follows: 


’QT  ‘ 

Q11 

+  1  Q?2  ' 

0*  = 

= 

-Q*_ 

i 

D 

tNJ* 

*  iQh. 

This  is  also  a  function  of  x  and  t  to  C^,  and  x  is  on  [a,  b]  .  Q*j  and 
Q*2  are  the  components  of  force  along  the  X^- ,  X^-axes,  respectively; 
and  Q*j  and  are  components  of  torque  along  the  -X^-,  X^-axes, 

respectively. 

External  forces  on  (0,  H  )  will  be  restricted  to  a  series  of  concentrated 

'  n 

driving  forces  which  are  caused  by  the  mass  eccentricities.  The  force 
vector 


P* 


p* 

ii 

P11  +  1  P12 

p*_ 

_P!l  +  iP22_ 

denotes  the  "driving  fore.",  a  function  of  x  and  t  to  where  x  is 


on  (0,  Q-n),  which  vanishes  at  a  finite  set  of  points.  These  points  are 


located  at  x  =  a^;  j  =  1,  2, 


;  n  and  k  =  1,  2, 


;  k(j),  where 


k(j)  means  that  k  may  be  a  different  integer  for  different  spans.  (See 
Figure  1 .  ) 


Similarly,  "internal  force"  as  a  force  vector  i~  defined  as  follows: 


F* 


'F1  ' 

F* 

2  _ 

F*  +  i  f* 
f 11  1  r  12 


F21  +  1  F22 


This  force  vector  is  equal  to  the  force  applied  to  the  differential  mass 
element  at  x  -  dx/  2.  F  is  a  piecewise  continuous  function  of  x  and 
t  to  C^,  where  x  is  on  [a,  b]  ,  and  the  magnitude  of  each  upward  jump 

at  a  point  of  discontinuity  equals  the  external  concentrated  force  at  that 
point.  (See  Figures  2  and  3.  ) 


Throughout  the  following  equations,  a  quantity  with  "#"  above  it  is,  unless 
otherwise  specified,  a  2  X  2  matrix  with  complex  entries,  and  the  inverse 
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of  such  a  matrix  will  be  denoted  by  attaching  the  superscript  "-1".  A 
tilde  (perhaps  together  with  other  symbols)  above  a  quantity  denotes  the 
Laplace  transform  of  that  quantity;  e.  g.  , 


s)  = 


°o 

J  g(x,  t)e 
0 


-  st 


dt 


where  g  is  a  function  of  x  and  t  to  C ^  where  x  is  on  [a,  b]  . 

Assume  throughout  that  ¥  (x,  t)  =  0  and  dY  (x,  t)/  dt  =  0  at  t  =  0 
(ft  is  the  null  vector  in  C Physically,  this  means  that  the  system 


is  initially  at  rest  or  that  the  initial  transverse  deflection,  angular 
deformation,  transverse  velocity  and  angular  rotation  are  zero;  the 
initial  position  of  the  shaft  is  coincident  with  the  Xj  -  coordinate  of  the 


fixed  reference  frame.  Thus,  the  solution  will  be  limited  to  steady- 
state  conditions  by  applying  the  Laplace  transform  technique  to  solve 
the  governing  differential  equation,  and  the  transient  state  is  assumed 
to  vanish  automatically. 

Let  Q  be  an  external  force  vector  applied  to  d(x)  as  defined  before; 

let  Y  be  the  position  vector  of  d(x).  Because  of  axial  symmetry,  the 

relation  between  the  force  vector,  Q  ,  and  the  position  vector,  Y  ,  may 
be  written  in  either  of  the  two  following  forms: 

Q  =  s  Z  Y 


or 


#  ~ 
M  Q 


where  "  ~  "  denotes  the  Laplace  transformation  variable. 


# 

Z 


# 


and  M 


are 


called,  respectively,  impedance  and  mobility  in  2  X  2  matrices  with  com¬ 
plex  entries.  If  ^  and  &  are  non- singular,  $  M  =  Li -A,  where  f 
is  the  2X2  identity  matrix.  Impedance  and  mobility  are  functions  of  s, 
and  their  evaluation  depends  solely  on  the  parameters  of  the  dynamic 
system  itself.  Impedance  is  usually  referred  to  as  a  transfer  function 
in  mathematics.  It  is  assumed  that  the  impedance  or  mobility  at  every 
support  location  is  a  known  quantity  or  can  be  calculated. 
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CHAPTER  2 


TRANSMISSION  LINE  ANALOGY  SOLUTION  OF  THE 


SU PERCRITICAL  SHAFT 

WITH  ANY  NUMBER  OF  INTERMEDIATE  SUPPORTS 


CONFIGURATIONS 


A  physical  model  of  the  shaft,  which  corresponds  to  the  mathematical 
model  in  Figure  1,  is  presented  in  Figure  4.  All  the  parameters  used 
in  the  following  theoretical  analysis  will  be  assumed  to  be  positive  in 
the  sense  indicated  in  the  figure.  The  shaft  end  conditions  and  the 
conditions  of  all  intermediate  supports  may  be  described  as  impedances 
or  mobilities;  however,  for  convenience,  only  impedances  will  be  used 

# 

throughout.  The  impedances  of  all  supports  are  denoted  as  Z  (0) , 

#  #  U 

Z  .(£.),  .  .  .  ,  Z  ((  )  where  0,  i  f.  .  .  .,  C  indicate  the  locations 
11  n  n  12  n 


of  supports  measured  from  the  left  end  support.  All  the  support 
impedances  are  assumed  to  be  given. 


P  <ajk) 


. ,  ,  th 

is  the  k 


driving  force  in  the 


.th 

J 


span  where  j  =  1 ,  2,  .  . 


n 


and  k  -  1,  2 . k.  The  value  of  k  may  be  different  in  different  spans, 

depending  on  each  span  length  and  on  how  many  differential  mass  elements 
are  assumed  to  exist  in  each  sp  n.  Hence,  k  will  be  written  as  k(j), 
which  means  that  the  value  of  k  depends  on  the  span  number. 


EQUATIONS  OF  MOTION  FOR  THE  SHAFT 


The  equations  which  govern  the  behavior  of  rotating  shafts  include  the 
effects  of  gyroscopic  motion,  rotational  inertia,  and  shear  deformations. 
The  projections  of  an  infinitesimal  element  of  the  shaft,  of  length  dx 
and  with  all  the  forces  acting  on  it,  are  shown  in  Figures  2  and  3.  The 
symbols  used  in  these  two  figures  are  listed  as  follows: 

*  =  the  unnormalized  variable 

co*  =  angular  velocity  of  the  rotating  shaft 

I  =  moment  of  inertia  of  the  shaft  cross-sectional 

area  about  a  diameter 

p  =  mass  density  of  the  shaft  material 

A  =  area  of  the  shaft  cross  section. 

Under  certain  conditions,  such  as  at  very  high  rotating  speeds,  not  only 
the  centrifugal  forces  of  the  infinitesimal  rotating  mass  but  also  the 


9 


X 


10 


movements  of  the  axes  of  the  rotating  mass  are  of  importance  and  should 
be  taken  into  account.  The  mathematical  derivations  of  the  expressions 
for  the  moments  due  to  angular  movements  (rotatory  ar  gyroscopic 
effects)  and  the  expressions  of  inertia  forces  are  included  in  Appendix  A. 

If  the  projection  on  the  X^X^-plane,  which  is  treated  as  a  real  plane,  is 


considered,  the  following  two  equations  can  be  obtained  by  applying 
equilibrium  conditions  and  by  neglecting  higher  order  terms  (see  Figure 
2): 


By  summing  moments, 


I  Y* 

1  1  2  ltt 


-  2 


P*“o  Y22t 


=  0 


(1) 


By  summing  forces, 


Y* 

1  ltt 


=  0 


(2) 


By  applying  shear  relationship, 

F*,  =K.AEs(-Y»ix-  Y*i>  ; 


where 

K'  =  numerical  factor  depending  on  the  shape  of  cross  section 


average  shearing  stress  on  cross  section 
shearing  stress  at  neutral  axis 


E 

s 


-YlVY21 


q  ,  where  Q  is  the  moment  of  cro6s-sectional  area 

of  shaft  above  the  neutral  axis  with  respect  to  this  axis, 
and  d  is  the  thickness  of  section  at  the  neutral  axis,  in 
this  case  the  diameter  of  the  circular  section 

shear  modulus  of  shaft  material 

slope  of  deflection  curve 

slope  of  deflection  curve  when  shearing  force  is  neglected 

change  of  slope  deflection  curve  due  to  shearing  force 
only. 


The  preceding  equation  may  be  rearranged  as  follows: 


By  applying  the  bending  relationship, 


E  1 

y 


moment  on  the  differential  mass  element 

change  of  slope  of  deflection  curve 

Young' s  modulus  of  shaft  material. 

The  above  equation  may  be  rearranged  as 

F*  +  E  IY+  =0.  (4) 

l.  1  y  2 1  x  '  ' 

Similarly,  if  the  projection  on  the  XjX^-plane,  which  is  treated  as  an 

imaginary  plane,  is  considered,  and  if  the  preceding  arguments  are  used, 
the  following  four  corresponding  equations  are  obtained  (see  Figure  3): 

F12  "  F22x  '  pIY22tt  +  2pI“0Y21t  =  ° 


where 


F* 

r  21 

Y* 

21x 


FTzx  +  A  YTa,t  =  0 

f*2  +  k- aes(y*2x  +  Y*2)-o 


F* 

*  22 


E  I  Y*  =  0. 
y  2Zx 


If  each  pair  of  equations  is  considered,  in  complex  variable  form  the 
corresponding  four  equations  are  obtained: 


(1*1  +  1 1  *2)  -  (f  *lxT  iF*2x)  -  pI(Y*Ut+  lY*2tt)T  2lPI“*  (Y*lt  +lY*2t)  ° 


(FllxtiF12x)  +  PA(YT.UtiY?«,>  =  ° 

<FTl  +  iF?2>  +  K'  AES  I  <YTlx  +  ‘Y?2x>  +  <Y*1  1  iY?2> 

^F*l  +  iF22^  +  Ey  '^Y21x^'Y22x^  • 


=  0 


If  the  notations  for  position  vector  and  force  vector  as  defined  on  page  8 
in  Chapter  1  are  applied,  the  equations  may  be  written  in  more  compact 
forms  in  referring  to  the  framework  of  a  2- dimensional  vector  space 
over  a  field  of  complex  numbers  as  opposed  to  the  more  general  space 
of  four  dimensions  over  a  field  of  real  numbers. 
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F*  -  F* 


2x  P,YStt  +  2iPI“0Y*t 


F?X  +  p  AVTtt  - 0 

F*  +  K'  AE  (Y*  +  Y*)  =  0 


F*  +  E  IY*  =  0  . 

Z  y  2x 

For  normalizing  these  equations  of  motion,  R^.  E^A,  an<3  R^/  Cs  are  u8ec*' 

they  correspond  to  unit  length,  unit  force,  and  unit  time,  respectively. 

R^>  cg  are,  respectively,  radius  of  gyration  in  bending  and  the  velocity  of 

sound  for  the  shaft  material.  Note  that  cg  =  "\JE^/p.  By  applying  the 

standard  procedures  for  normalization,  the  differential  equations  of 
motion  may  be  written  as  follows: 

F.'F2x-  Y2tt+2i“oY2t  =  ° 


F,  +  Y„,  =  0 
lx  ltt 


F2  +  Y2x  *  0 


e1  Fj  +  Yjx  +  Y2  =  ° 


(5) 


where  e'  =  ( l/  K'  )(E Eg),  which  is  a  numerical  coefficient  depending  on 
the  shape  of  cross  section  and  the  shaft  material,  and  where  is  the 

angular  velocity  about  Xj-axis,  as  before  Yj,  are,  respectively,  the 

transverse  position  and  the  inclination  of  normal  of  the  elastic  curve  of 
shaft.  Since  all  the  initial  values  of  the  parameters  appearing  in  Eqs.  (5), 
i.  e.  ,  the  right-hand  limits  at  t  =  0 ,  are  null  vectors,  Eqs.  (5)  in  the 
Laplace  transform  are  as  follows: 

/v  rw 

F1  •  F2x  +  s<2i“0  -  s>  Y2  =  0 


F,  +  8  Y ,  =  0 

lx  1 


F2  +  Y2x  =  ° 


e'  F,  +  Y,  +  Y,  =  0 
1  lx  Z 


(6) 
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In  Eqs.  (5),  if  Fj,  F^,  and  are  eliminated,  the  governing 

differential  equation  in  terms  of  Yj  (transverse  deflections)  of  the  dynamic 
system  may  be  obtained  in  the  following  form: 


Ylxx*x  +  YU.-<1+e 


>Ylxxtt+e' 


ltttt 


+  2io> 


,<Y1 


xxt 


-  e 


Ym.>  * 


(7) 


If,  as  has  been  noted,  all  the  initial  values  for  Yj,  Y^,  Yj  ,  and  Y 

are  zero,  and  if  the  Laplace  transform  is  used  again,  Equation  (7)  may 
be  transformed  as  follows: 


Y,  +  s 
lxxxx 


2ia>  -  s(l  +  e'  ) 


Y.  +s 
lxx 


1  +  se'  (s  -  2ico  ) 


Yl=0 


(8) 


The  solution  of  Eq.  (8)  corresponds  to  the  steady-state  solution  of  Eq.  (7), 
in  which  a  solution  is  assumed  in  the  form  of 


!  =  y(x)< 


iwt 


where  w  is  the  frequency  of  the  driving  force  and  y(x)  is  of  the  form 

y(x)  =  AemlX  +  Bem2X  +  Cem3X  +  Dem4X 

where  m^,  m^,  m^,  and  are  roots  of  the  characteristic  equation  of 

the  governing  differential  Eq.  (7).  A,  B,  C,  and  D  are  arbitrary 
integration  constants  which  may  be  determined  by  applying  boundary 
conditions. 


SOLUTIONS  FOR  THE  EQUATIONS  OF  MOTION 


With  reference  to  Eq.  (8),  let 


Yj  =  Aemix  +  Bem2X  +  Cem3x  +  Dem4X 


Then,  if  the  above  expression  is  substituted  into  Eq.  (8),  the  following 
characteristic  equation  is  obtained: 


m^  +  s 


2iu>  -  s(l  +  e'  ) 


2  , 

m  +  s 


1  +  se'  (s  -  2i«n) 


=  0 


which  yields 


m  j ,  m 


2 


±  ie  j  V^" 


where 


a)  =  -is  or  s  =  iw  (w  is  the  frequency  response  which 
is  the  same  as  the  frequency  of  driving  force,  i.  e.  , 
the  angular  velocity  of  the  rotating  shaft  itself) 
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=  j  1/  +  e'  <0 


1  +  Qf^  +  a 


ot  =  wQ  +  ~  a)  (e'  -  1) 


m  m  =  ±  e2Y“ 


where 


/e^  -  e'  cj 


The  replacement  of  s  by  iu  corresponds  to  the  steady- state  solution. 

If  the  relationships  in  Eqs.  (6)  are  applied,  the  following  general  solutions 


rsj  rss 


for  Y^,  Fj,  are  obtained: 

Yj  =  Aemlx  +  Bem2x  +  Cem3x  +  Dem4x 


y2  =  . 


i_Aemlx+— i_  Bem2X--l,-Cem3xt  II  De”1^ 


ele3 


ele3 


Fj  =  <*>' 


F2  =  u 


3/2  ,-i_Aem‘x  +  — i— Bem2X+— Cem3x--i-  Dem4x 


—  Aem|x - —  Bem2x  +  e,Cem3x  +  e,Dem4x  . 

e^  33 


If  the  following  set  of  arbitrary  constants  (q^,  q^,  r^,  r^)  is  introduced: 


1  p+1 

— w  r  e , e,r , 
s  13  1 


„  i  p+ 1 

B  =  -V“  ele3ql 


C  =  -  —  u>  p+1e,r, 
s  2  2 


i  p+ 1 
1" U  e2^2 
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where  p  is  an  arbitrary  real  number  which,  once  chosen,  remains  the 
same,  and  where  q^,  q^,  r^,  r^  are  complex  numbers;  therefore,  the 


expressions  for  Yj,  Y^,  Fj,  can  be  rewritten  as  follows: 

\£>x  -e2'V« 


0  1  p 

Y  =  — w  h 
1  s 


-  ie  l^wx  -  e£vw  x 

e^uie  1  q j  -  le^w  e  q2 


iei^x  ep  'Jw  x 

ie  j e^e  1  r] '  ie2We  a  r^ 


rv  i 

Y,  =  -LUP 

2  s 


3/2  -iel^  x 


w  '  e  qj  -  le^  a 


3/2  -  ez’Ju  x 


U) 


3/2eie1^xri+ie^3/2ee2'^x^ 


F  =  « 


3/2  -ieiN/ux  ,,3/2  -ep'fu  x  , 
xe^w  '  e  1  q,  +  «  e  £  q  + 


3/2  ieM*>x  3/2  x 

le^w  e  1  r  j  -  w  '  e  4  r^ 


F-,  =  to 


1  'Vu 


-xei^u)  x 

elWe  qi-e2e3u>e 


L  *2  + 


ie  1^ 


n/u) 


_  j  1  vu  x  e?  Vu>  x 

elu  e  ri  "  e2e3U>  e  r2 


(9) 


By  defining 


_<*1  “ 

and 

r  = 

'  rf 

_  q2  . 

r2_ 

as  two  arbitrary  constant  vectors  in  a  2- dimensional  vector  space  over  a 
field  of  complex  numbers,  and  by  recalling  the  expressions  of  the  position 
and  the  force  vectors  as  defined  under  "Assumptions  and  Notational 
Conventions"  in  Chapter  1,  Eqs.  (9)  may  be  rewritten  in  a  much  more 
compact  form,  as  follows: 
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(10) 


at 

Y 

1 

6  y 

ft 

0 

i(x) 

# 

0 

q 

at 

# 

# 

f 

# 

it 

#  ,  t 

F 

0 

cf 

C21 

C22 

0 

R(-x) 

r 

tilt*  f 

where  C  ,  C^,  C j  ^ ,  Cj^,  ^21’  anc*  ^22  are  as  s^own  *n  Appendix  B,  and 
$  is  the  2X2  null  matrix. 


The  "propagation  matrix",  R(a)f  is  a  function  on  the  real  numbers 
defined  by 


R(a) 


cis 


a) 


where 


cis  0 


i9 

e 


exp  0  =  e 

#  ,  v  # 

R(0)  =  I 


c  os  0  +  i  sin  0 


0 

exp(-e2^>  a) 


Simply  by  expanding  the  matrix  form  of  Eq.  (10),  this  equation  may  be 
verified  as  being  identical  to  Eqs.  (9). 


BOUNDARY  CONDITIONS 


By  applying  the  boundary  conditions  at  x  =  0,  i  j,  J?^,  .  .  .  ,  and  a^ 
on  the  Xj-axis,  the  integration  constants  appearing  in  the  general  solution 
of  the  governing  differential  Eq.  (8)  can  be  determined.  (See  Eq.  10.  ) 

At  Leit  End  Support 

The  boundary  condition  at  x  =  0-0,  i.  e.  ,  the  left  end  support,  may  be 
obtained  by  considering  the  impedance  at  this  point.  (See  Figure  5. ) 

at 

Since  F  at  x  =  0-0  is  in  the  negative  direction,  according  to  the  sign  con¬ 
vention  defined  in  Chapter  1,  the  boundary  condition  of  the  left  end  is  ex¬ 
pressed  as: 

-  “(0)  =  s  Z0(0)  “(0)  (11) 

which  should  be  satisfied  at  x  =  0-0;  -0  denotes  the  left-hand  limit. 
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*S  °'°*:C 


I - 1 


F(0)  F(0) 


Figure  5.  Left  End  Support 


At  Right  End  Support 

The  boundary  condition  at  x  =  i^+0,  i.  e.  ,  the  right  end  support,  may  be 
obtained  in  a  similar  way,  as  follows: 

”<<„>  =  »  "<f  »>■  <12> 

at 

Note  that  F  is  in  the  positive  direction  at  x  =  f  n+0;  also,  +0  denotes  the 

right-hand  limit.  Eq.  (12)  is  the  condition  that  should  be  satisfied  at 
x  =  ^n+0.  (See  Figure  6. ) 


*  *  8  n-0  x  *  i„  ♦  0 


77777 


# 

Zn<«n> 


F<in>  mn> 


Figure  6.  Right  End  Support. 
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At  Intermediate  Supports 


The  boundary  conditions  at  x  =  H.,  j  =  1,  2,  .  .  .  ,  n-1,  may  be  obtained 

J 

by  considering  the  continuity  of  the  deflection  curve  and  the  force  equilib¬ 
rium  at  every  intermediate  support.  (See  Figure  7.) 


x  =ii  -0  *  =£j  +  0 

I  A  I  >  I 


# 

zi'« 

J» 

CZZZ1 


3 


F(fi  J  -o)  Fit  j  +0) 


Figure  7.  Intermediate  Supports. 


If  the  impedance  is  isolated  at  x  =  C. .  and  if  the  net  force  applied  at  x  =  i. 
r-  -i  J  J 


is  [r(^-o)-  f(*.+o) 


iind  if  the  deflections  at  x  =  £.-0  and  at  x  =  t.  +0 

:  j 


are  the  same,  the  boundary  conditions  at  x  =  i.  are  as  follows: 

3 

“(£.-0)  =  Y(£.+0)  =  Y(£.) 


F(£.-0)  -  F(£.  +0)  = 


(13) 


j  =  1,  2,  .  .  .  ,  n-  1 

where  +0  and  -0  denote  right-hand  and  left-hand  limit,  respectively. 

At  Driving  Forces 

The  boundary  conditions  at  the  locations  where  driving  forces  are  being 
applied  may  be  obtained  in  a  similar  manner.  If  it  is  assumed  that  a 
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single  concentrated  driving  force,  P  (a.,  )  is  applied  at  x  =  a.,  (i.  e.  ,  the 

jk  jk 

th  th 

k  driving  force  in  the  j  span),  the  corresponding  boundary  conditions 
at  x  =  a^  are  as  follows  (see  Figure  8): 


Y  V°>  =  Y  (V0)  =  Y  V 


(14) 


F  (ajk-0)  +  P  (ajk)  -  F  (ajk+0)  =  0 


REFLECTION  MATRICES  AT  SUPPORTS  IN  TERMS  OF  IMPEDANCES 


If  all  the  boundary  conditions  except  those  at  x  =  a.,  are  applied,  as 

jk 

mentioned  in  the  preceding  section,  some  useful  functions  can  be  found 
during  the  process  of  evaluating  boundary  conditions.  (See  the  complete 
mathematical  derivations  in  Appendix  C.  ) 


At  x  =  0  on  the  Xj-axis,  the  reilection  matrix  (looking  to  the  left  at 


x  =  0)  is  defined  in  terms  of  end  impedance  by 


K 


,(»)  = 


« 


-  CZ2+“0(G>C1Z 


-1 


A 


c2rV0)cn 


(15) 


where 


A  # 

Z  =  E  Z  E, 


1  0 

0  - 1 


# 

z 


# 

F  Z  F, 


0 


\Ju> 


At  x  =  i.  j  ,  H  ^ , 


the  left  are  defined  as  follows: 


^n-1’  re^ec^^on  matrices  looking  to 
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Left  End  Support 


P(ai|()  (assumed  in  a 
i  positive  direction) 


^l£  [’^22+*U^1^12]  [^2l“  *11^1^11 

^21^2^  ['C22+  Z2f^2^iz]  [C2l“  ^24^2^11 

i  %  r  i  a  ,*  j  i'1n  a  j  i  U16) 

rG-l)^j-l)=L’£22+  *(j- lK^j-  1)C12J  [C2l'  Z(j-0<^j-l)CllJ 

•  • 

•  • 

^(n-l)£^n-l^=  "C22+  ^(n- 1)«  (fn- 1  *C  1 2  C2l"  ^(n-  l)£^n-  1^C1 1 


where  subscripts  1,  2,  .  .  .,  j-1,  .  .  .  ,  n- 1  indicate  the  locations  of 
supports;  subscript  I  indicates  that  the  reflection  matrices  are  being 
looked  at  to  the  left  of  x  =  t  j,  j,  ...»  I  |  , 

®l|(*l).  ^2£^2^’  •  •  •  ’  *  ■ 

^(n  l)£^n  1^’  rePre8ent»  respectively,  the  total  impedance  looking  to  the 

left  at  x  =  £j,  . £j,  ^2'  •  •  ••  £-_j.  •  •  »  1’  which  are 

defined  as  follows:  , 


* 
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z]£(*j)  "  2  l(^  1 )  + 

z  ZI^t)  ~  z2^ 2)  +  z  l£^2^ 

?  (17) 

Z(j-i)je(£j-i)  =  *j-i^j-^  + 

Z(n-l)C^n-l^  =  Zn-l^n-P  +  Z(n-2)£^n-P 


where  the  terms  ^  Z  1^2^ . Z  (j-2)f^j-  1^ . Z(n-2)i^fn-l^ 

will  be  explained  in  the  next  section. 

At  x  =  fn  on  the  Xj-axis,  the  reflection  matrix  (looking  to  the  right) 
is  defined  in  terms  of  the  end  impedance  by 

#  r#  #  #  1  - 1  r #  #  # 

r  (*  )  =  -  c77+  z  (i  )c,,  c71-z  (£  )c,,  . 

n'  n '  (_  22  n'  n'  12J  |_  21  n'  n'  1 1 

At  x  =  (■  , ,  f  the  reflection  matrices 

n- 1  n- c.  j  i 

looking  to  the  right  are  defined  as  follows: 
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# 

r,  (f  ,)= 

(n-  1  )r'  n- 17 


"  #  #  # 
C22+  Z(n-  l)r^n-l^C12 


"#  #  # 
C21 "  Z(n-l)r^nVCll 


r(n-2)r^n-2^= 


"  #  #  f 

_C22+  Z(n-2)r'^n-2^C12 


"#  #  # 
C2 1  ”  Z(n-2)r^n-2^Cll 


# 

"  #  #  #  " 

-  1 

"#  §  #  ‘ 

r.  (£.)  = 

jr'  j 

_C22+  Zjr^j^C12 

C,.-  z.  ( £.)C ,  . 

L  21  Jr  J  nJ 

>(19) 


# 

#  #  » 

-  1 

#  # 

►1 

•-1 

It 

'C22+  zlr^fl>C12_ 

C2 1 "  Zlr^l^Cll_ 

where  subscript  r  indicates  that  the  reflection  matrix  is  being  looked  at 
to  the  right  of  x  =  j,  ...,  ly  ... ,  i  ^ . 

Z(n- l)r^n- 1^  z(n-2)r^n- 2^’  •  zjr^j^  •  •  •*  zir^i)  represent, 


Jr 

respectively,  the  total  impedance  looking  to  the  right  of  x  =  £ 


n-  r 


f  £ 

n-2’  •  •  •’  V 

.  .  ..  *j. 

which  are  defined  as  follows 

i,  , V  (£  .)  = 

(n- 1 )r'  n-  1 7 

&  ,(£  ,) 
n- 1 '  n-  1 7 

+ 

i  (£  ,) 

nr'  n- 1 

i,  (1  ,)  = 

(n -2)r'  n-2' 

f  ,<e  ,) 

n-2'  n- 27 

+ 

i,  ..  (e  ,) 

(n-l)r'  n-27 

>(20) 
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1  (£.)  =  t.(£.)  +  (£.) 

jr  J  j'  y  G  +  1)r  J 


=  ^  i  i )  +  ^2r^V 

where  the  terms  z  (£  , ),  I ,  , »  (£ 

nr'  n-1'  (n-l)r'  n-2" 

^2r 1 )  be  explained  in  the  next  section. 

0.  « 


^(j  +  O^V . 


For  each  closed  interval 


n- 1  n 


t  £ 
r  2 


£  £ 
*• :  l  » 


#  •  •  •  i 


j-Jl  jj 

on  the  Xj-axis,  the  corresponding  reflection  matrices  looking 
to  the  left  are,  respectively,  defined  as  follows: 


#  #  #  # 
r0(x)  =  R(x)r0(o)R(x) 

#  #  #  # 
rlJf(x)  =  R(x-£j)r 


#  #  #  # 
r(j-l)£(x)  =  R(x"*j- l^rG-l)£  (*j-l)R(x’*j_i) 


>  (20 


r(n-l)£^X^  '  R^X-ln-l^r(n-l)£^n-PR^x"fn- 1^ 


For  each  closed  interval 


o,  £ 


- 

“ 

r  >  i 

£  .,  £ 
n- 1  n 

f 

£  ,,  £  , 
n-2  n- 1 

• 

•  •  l 

i 

i 

* 

on  the  Xj-axis,  the  corresponding  reflection  matrices 
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looking  to  the  right  are,  respectively,  defined  as  follows: 


rn(x)  =  R(«n-x)rn(tn)R«n-x) 

r(n-l)r^X^  =  R^n-rX^r(n-l)r^n-l^R^n-rX^ 


#  #  #  # 
r.  (*)  =  r(*.  )r.  (£.)r(I.-x) 

jr'  '  '  j-x'  jr'  y  '  j  ' 


P  (22) 


#  #  #  # 

rlrM  =  R(|rx)rlr(<1)R(<1-x)  . 

IMPEDANCES  IN  TERMS  OF  REFLECTION  MATRICES  AND  SHAFT 
CHARACTERISTIC  IMPEDANCE 


It  is  shown  in  Appendix  D  that  £nf(x),  the  total  impedance  looking  to  the 

0,1,1  ,  is  related  to  ^Tq(x) 


Ol 

left  of  a  generic  point  x  where  x  is  on 
by  the  equation 


™  — 

A 

r#  #  #  -| 

# 

#  i  1 

On  0,  fj 

•  *0«w  • 

C21  +  C22*VX^ 

_  - 

cn 

+  c,2r0W 

Similarly, 


<-# 


#  # 


r# 


#  #  1-1 


r  .  n 

A 

r#  #  #  -j 

r#  #  #  -| 

On 

1 

1 

_ J 

’  z(j-l)^x)  = 

C21  +C22r(j-l)^X^ 

C11  +C12r(j-1)£^ 

•  * 

"#  #  # 

'#  #  # 

-  1 

On 

n- 1  n 

’  Z(n-  l)£^X^ 

C21  +C22T 

C11  +C12r(n-1)£M 

• 

/ 

It  is  also  shown  in  Appendix  D  that  the  other  similar  set  of  equations 
which  represent  the  total  impedance  at  x  looking  to  the  right  is: 


r  i 

u 

"#  #  # 

"#  #  #  1 

On 

i 

c 

1 

C 

_ i 

,2  (x)  = 

nr' 

C21  +C22Fn^x^ 

C11  +C12rn^ 

"#  #  # 

"  #  #  # 

-1 

On 

£  £ 

L  n-2*  n-lj 

•  !(„-«)rW  = 

_C2 1  +C22r(n-l)r^X^_ 

.Cll+C12r(n-l)rW 

# 

■#  #  # 

"#  #  # 

On 

1.  ,,  i. 

L  j-1  jJ 

,  z.  (x)  = 
jr' 

C21  +C22rjr(x) 

c. ,  +  c. 7r.  (x) 

|_  11  12  jr  J 

u 

'#  #  # 

'#  #  # 

On 

0, 

•  *i,W- 

C21  +C22rir^X^ 

C11  +C12rir^ 

The  expressions  for  ^(fj),  £j£(*2).  •••  >  •  •  •• 

^(n  2)1  ^n  1^  *n  0^)»  anc*  the  expressions  for  ^nr(^n  j)> 

^(n-l)r<<n-2) . ^(j+l)r(£j) . in  E19'  (20)  can  be 


obtained  from  Eqs.  (23)  and  Eqs.  (24),  respectively.  Eqs.  (23)  and 
(24)  will  also  enable  us  to  evaluate  the  impedance  looking  both  to  the 
left  and  to  the  right  at  any  location  along  the  shaft. 
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Let 


#  i  #  -l 
Z8-  U21U11 


which  is  expanded  in  complete  matrix  form  in  Appendix  B,  and 


i 


# 

F  Z  F, 
s 


F 


(25) 


where  Z  is  called  the  characteristic  impedance  of  the  rotating  shaft.  Note 

#  8 

that  Z  is  a  function  of  frequency.  The  significance  of  the  characteristic 
impedance  will  be  discussed  in  Chapter  3. 

COMPLETE  SOLUTION  IN  MATRIX  FORM 


If  all  the  boundary  conditions  beginning  on  page  17  of  this  chapter  are 
used,  all  the  integration  constants  appearing  in  the  general  solution  can  be 
determined.  (See  the  complete  mathematical  derivations  in  Appendix  E.  ) 
The  complete  solution,  from  which  one  may  calculate  the  dynamic  response, 
i.  e.  ,  deflections  and  internal  forces  caused  by  the  single  driving  force 


P(a.,  ),  at  any  location  along  the  shaft,  may  be  written  in  the  following  matrix 
JK 


form: 


On 


£  £ 
1’  2 


Y 

. 

1  # 

—  C 
s  y 

# 

0 

#cll 

*12 

#  # 

R(x)r  (0)  r2 

F 

i 

A, 

*22. 

1 

id*: 

i 

X 

rv 

i _ 

(26) 
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On 


£  £ 

j-2’  j-1 


On 


a.,, 

J 


L 

J 

~ 

J_ 

# 

s 

c 

0 

- 

s 

y 

# 

# 

- 

F 

0 

cf 

1 

*. 

J 

-  1  ’ 

"jk 

t 

r* 

- 

_ 

r%j 

_1_ 

# 

# 

Y 

-c 

0 

~ 

8 

y 

# 

# 

F 

0 

cfJ 

T 

- 1 

o*, 

J 

i 

iH 

1 

1 

o 

«»< 

rvT 

i _ 

#  # 

C21  C22 

R(-x)r.  j 

# 

# 

#  ,#  , , . 

C11 

C12 

R(x)r(j-i)f(°)  rj< 

# 

C21 

# 

C22  ^ 

i 

ja=* 

i 

►<i 

w. 

l 

Y 

1  # 

—  C 
s  y 

o5*? 
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rj-2~R^j-2^  C1  l^(j  -  3)£^j-  2^+Cl  2  Clir(j-2)je^j-2^+Cl 2  R^“^j-2^rj-  1 

rj-l=^(^j_l)  ^ll^(j-2)£^j-  1^+C12  Cll^(j-l)f^j-  P+C12  R(‘£j-  l^rjjf 


rjr^(ajk) 


-1 


r#  .  J  # 

T.  a.JC.+C 

jr'  jk'  + 


#  _  1  - 
cf  PV 


qjr=R<-ajk> 


?-f'(j-»((aj^Jr(aJk)]‘?r(J-0«(V*-+*C+ 


#  -  1-  ,  . 
Cf  P<ajk> 


^ll+C12^(j  +  l)r*V 


-  1 


Cll+^12f'jr^j^ 


#  i. 


R(JP.  q. 

J  jr 


V2=R(-ij  +  l) 


^ll+^12^(j  +  2)r^j  +  ^  Cll+^12^(j  +  l)r^j  +  l^ 


R(^j  +  l)qj  +  l 


#  #  £ 


q  ,  -R(-J0  J  C.^+C.X  n  (£  ?) 

nn- 1  '  n-2'  11  12  (n-l)r'  n-2' 


-  1 


Cll+C12^(n-2)r^n-2^ 


R^n-2)V2 


#  #  # 


-1 


'#  #  # 


q  =R(-4  .)  C-.+C.^r  (*  ,)  C11+C1,r/  n  (4  ,)  R(*  ,)q  .  . 

nn  '  n- 1  l_  11  12  n'  n-1'  11  12  (n-l)r'  n-1'  '  n-l'm-1 

(26a) 

If  a  superposition  technique  is  used,  the  total  dynamic  responses  caused 
by  all  the  driving  forces  distributed  along  the  shaft  can  be  obtained. 


SOLUTION  IN  WAVE  FORM 


Proceeding  with  a  lengthy  matrix  algebraic  manipulation  of  the  resulting 
Eqs.  (26),  as  shown  in  the  preceding  section,  and  applying  the  super¬ 
position  technique,  one  may  express  the  total  dynamic  response  caused 

by  all  driving  forces  P(a.,  ),  j  =  1,  2,  .  .  .  ,  n  and  k  =  1,  2 . k(j) 

jk 

along  the  shaft  in  the  following  traveling  wave  forms.  (See  the  complete 
mathematical  derivations  in  Appendix  F.  )  Explanations  of  the  traveling 
wave  form  solutions  will  be  given  in  detail  in  Chapter  3. 
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In  the  span,  or  on 


Y 

F 


±1 
s  y 

0 


£  £ 
i-l’  i 


# 

o 

# 

c, 


# 

c 

# 

c 


,  i  =  1,  2,  .  .  .  ,  n 


# 

C. 


11  12 

ft 

21  C22 


>» 

U;(X) 

»> 

where 

»- 


U  (x)  =  total  wave  traveling  to  the  right  in  the  i^  span 

"  >» 

•I  Vx) 

j=l 


«< 


th 


U  .(x)  =  total  wave  traveling  to  the  left  in  the  i  span 


JL  «< 

*Z  Du« 


j=l 


If  j^i,  j  =  1,  2,  .  .  i-l;  then 


» 


U.  (x)  =  U..(x)  +  U..(x) 
JJ  JJ  JJ 


«< 


<s 


<r 


U..(x)  =  U . . (x)  +  U..(x) 

JJ  JJ  JJ 

->  k(j )  ^  ^  ^  ~ 

°j|W  =ZR(x-ajk)Cfc/lp(ajk) 

k=l 


H(X-ajk)-H(X-fj) 


k(j) 


#  a 


°jj(x)=Z  R(ajk-x)c.Cf'1P(ajk)  H(a.k-x)-H({.  rx) 


k=l 


(27) 
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»  #  r#  #  #  ]-i# 

uj.(x)=R(x-£j.1)[i-r(j.1),(C..1)r.r(£j.1)j  r0.1)4(*  ) 


X 


r#  §  >  < 

Rfvv.vvwvv^ii^-v^-^-v 


<<  # 

U.  ,(x)=R(f  .-x) 
JJ  J 


'#  #  # 

'-Wo-DftVJ 


-1# 

r.  (£.) 

Jr  j 


X 


[4(v+«(vfj-i)4-o«(Vi4(v.) 


H(x-f.  ,)-H(x-C.) 


»> 


nO+l)j(x)=R(x-£j) 


■##_,## 

I+Cl  lC12r(j  +  l)r^V 


-II 


»>  #  .  #  «< 

Jj'  J  11  12  Jj'  J7 


«< 


UG  +  l)j^'R(Vl‘x)rG  +  l)r^j  +  l,R(Vl"V 


#  #,#  # 

I+CHC12r(j  +  l)r^V. 


-1 


>» 


#_1#  «< 


[Wcmci*W. 


>>>  # 

U..(x)=R(x-^.  .) 
ij  J- 1 


##,## 

i+c:}c,,r.  (£.  .) 

11  12  ir'  l- 1 


-1 


»> 


#  _  j  #  «< 


[Vi)j(£i-i)+cncuV*)A-i)J 


«< 


#  # 

U..(x)=R(f.-x)r.  (*.-*.  ,) 
ij'  ’  '  i  '  ir'  x  l-l 


'#  #  #  # 

i+c:  c,,r.  (i.  ,) 

11  12  ir'  i-l' 


-1 


U*  <<K 


^(i-l)j^i.l)+CHC12^(i-l)j^i-l^ 
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If  j  =  i,  then 


>£> 


«< 


» 


TJ..(x)  =  U..(x)  +  U..(x) 
11'  11'  '  n'  ' 


« 


U..(x)  =  U..(x)  +  U..(x) 
11'  '  11'  '  11'  ' 


>  k(i)#  #  #  ^ 

t5i.(x)  =[R(wik|c(c;1P(.iki 
k=l 


H(x-a.k)-H(x-*.) 


<  #  #  =*  r 

Uii(x)  =  L  RKk-x)C-Cf  P<aik^  H(aik-x).H(f..1-x) 
k=l  L 


»  # 


U..(x)=R(x-l_1) 


#  # 


1-1# 


_I_r(i-  l^ir^i-l^J  r(i-  l)£^i-  1* 


r#  #  >  < 

r(£.-£.  ,)r.  (£.)u..(*.)+u..(£.  ,) 

'  i  i-l'  ir'  v  u'  v  ii'  i-r 


Hfx-^jJ-Hfx- 


«  # 

U..(x)=R(f .  ) 

ii'  '  '  i-x' 


’#  #  # 

/■rir^d-D^V 


-i# 


r.  (£.) 

ir' 


H(x-f.  ,) 


If  j  >  i,  j  =  n,  n-1,  .  .  .  ,  i+1;  then 

»>  »  > 

U..(x)  =  tJ..(x)  +  U..(x) 

JJ  JJ  JJ 


«< 


« 


u..(x)  =  tr..(x)  +  u..(x) 

jj  jj  jj 
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V0  -  L  R(-»Jk)c+Cf'  P(ajk) 

k=l 


H(x-ajk)-H(x-f.) 


6JjW  =  I  R(ajk-x)C.Cf-,p(a.k)['H(ajk.x).H(f..,-x) 
k=l 


»  # 


U.  .(x)=R(x-  ,) 

JJ  J"1 


#  #  # 


-1# 


R(*.-£.  ,)r.  (4.)U..(*.)+U..(4.  ,)  H(x-£.  ,)-H(x-f.) 

'  J  J-1  Jr  J  jj  J  jj  J- 1  J L  '  j-i7  '  j7J 


<<  # 

U.  .(x)=R(f.-x)| 
JJ  '  J  1 


i-r.  (f.)r,.  )£(f.) 

jr'  j7  (J-  1)KV  j7 


-1# 


r.  (*.) 

j  r'  j 


>  #  #  <- 

U..(f.)+R(*.-«.  jr,.  ..-(I.  ,)u..(£.  ,) 
jj'  J  '  J  J-1'  U*1)*  J-1  jj'  j-1 


H(x-£.  p-HCx-f.) 


»> 


l)j  =R^X‘£j-2)r(j-Z)£^j-2)R^j- l^j-2^ 


I+Cll  C12F(j-2)£^j- 


'#_!#  »> 


<« 


11  12  Jj'  j-17  jj'  j  -  1 7 


<<<  #  #  #  #  # 

U(j  -  1  )j =R^ j - 1  " . 1+C 1 1 C 1 2r(j  -  2)1 j  - 1  \ 


-1 


#  #  >>>  «< 

c:!c1?u..(£.  ,)+u..(i.  . 

11  12  jj'  j-17  jj'  j-1 
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where 

# 

0  when  p  <  0 

# 

I  when  p  >  0 

>  < 

U,  U  can  be  thought  of  as  "incident  waves"  traveling  to  the  right  and  left, 
respectively.  They  are  independent  of  shaft  support  conditions. 

»  « 

U,  U  can  be  thought  of  as  "reflected  waves"  traveling  to  the  right  and  left, 
respectively.  They  are  dependent  on  shaft  support  conditions. 

»>  «< 

U  and  U  are  the  resultant  of  the  incident  and  reflected  waves  traveling 
to  the  right  and  left,  respectively.  They  may  be  thought  of  as  "total 
waves". 
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CHAPTER  3 


MANIPULATION  OF  THE  SOLUTIONS 
TRAVELING  WAVE  CONCEPT 


The  steady-  state  wave  interpretation  of  the  solution  of  Eqs.  (27)  provides 
a  means  toward  visualizing  the  effects  of  support  parameters  (mass, 
spring,  damping)  and  shaft  characteristics  on  hypercritical  shaft  behavior. 

For  the  purpose  of  illustration,  only  one  driving  force,  P(a.,  ),  will  be 

jk 


considered  in  the  following  discussions. 

By  considering  the  jth  span  first,  i.  e.  ,  on 


£  £ 

j- 1*  j 


,  incident  waves 


traveling  along  the  shaft  are  considered  as  being  initiated  by  the  action 

~  >  < 
of  the  single  driving  force,  P(a^),  located  at  x  -  a.^.  U..^(x)  and  U..t(x) 


jjk' 


j  jk' 


are  defined,  as  before,  respectively,  as  incident  waves  traveling  to  the 
right  and  left  from  x  =  a  The  first  subscript  incidates  the  span  in  which 


the  wave  is  traveling;  the  second  subscript,  the  span  from  which  the  wave 
originates  or  the  span  in  which  the  driving  force  is  located;  and  the  third 


subscript,  the  location  of  the  driving  force.  For  example,  U  (x)  means 

that  this  wave  is  traveling  to  the  right  on  the  span;  also,  that  the 
original  wave  originated  in  the  jth  span,  and  that  the  driving  force  is 
located  at  x  =  a^.  It  can  be  seen  from  the  following  expressions  that 

the  incident  waves  are  independent  of  support  conditions.  (See  Eqs.  (108) 
and  (109)  in  Appendix  F.  ) 


# 


-  i~ 


U.jk(x)-R(x-a.k)C+C;  P(a.k) 


jk' 


jk' 


» 


#  ft 


-i: 


Ujjk(x)=R(ajk-x)C.c;  P(ajk) 


H(x-a.k)-H(x-C,) 


H(a.k-x)-H(«.  rx) 


(28) 


exists  only  on 


a.,,  £. 

Jk  J 


exists  only  on 


£  .  , ,  a., 
J-1  Jk 
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With  reference  to  Figure  9  and  Eqs.  (28)  and  (29),  the  propagation  of  the 
wave  from  the  point  of  application,  x  =  a.,  ,  of  the  driving  force 

-  Jk  «  # 

P(a^)  is  described  by  the  "propagation  matrix",  R(x-a^)  or  R(a^-x),  as 

the  wave  is  traveling  to  the  right  or  left,  respectively.  For  instance, 

>  > 

at  x  =  £.,  TJ. ..  (x)  =  U...  (£.)  has  traveled  a  distance  (i.-a.  )  to  the  i**1 
J  Jjk'  jjk'  j'  '  j  jk'  J 

support  of  the  system.  The  quantities  H(x- a.,  )- H(x- £ .)  and 

-i  J*  J 

H(a.,  -x)-H(£ .  i-x)  simply  indicate  the  regions  in  which  the  waves 

J  J  #  #  #  # 

exist.  The  quantities  C+C^.  and  C  C ^  perform  elementary  transforma¬ 


tions  or  operations  on  the  applied  driving  force  and  are  associated  with 
frequency  and  shaft  characteristics  as  defined  in  Appendix  B. 

» 

A  similar  interpretation  may  be  given  to  the  quantities  U.  ,  (x)  and 
«  JJk 

x),  which  are  regarded  as  reflected  waves  traveling  to  the  right 


and  left,  respectively.  They  may  be  expressed  as  follows  (see  Eqs.  (114) 
and  (115)  in  Appendix  F): 
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*V^>VVVV+Vfi-i>J 


H(x-< .  j)-H(x-C  ) 

J  J  _ 


(30) 


« 


Ujjk(x)  =  R(^-x) 


#  #  #  1-1# 
i-r.  (t)r,.  ...it j) 1  r.  (f.) 

jr  J  G'1)*  J  J  jr'  y 


X 


> 


%k<  v+«<  v  v .  4  - .  )c(£i- .  >^(v  i  >)  iH(x- f  j- .  h(x- v. 

(31) 


Both  exist  only  on 


£  a 
j-1'  j 


The  second  term  in  Eq.  (30)  or  Eq.  (31)  may  be  related  to  an  infinite 
series: 


•#  #  # 


1 


-11  ■ 


111' 


^j-i)«(£j-i4(Vi,4-.)«(£3-i)i(Vi4-i)‘(£j-i)i(£j-i)  + 


IV 


4o«(£j-i4(Vi4-iK(£j-i)i(£j-i4-iu(£j-i)i(Vi)  +  -  •  • 


# »  # 


-1 


rW - 11 


111 


#  #  #  #  #  #  # 

=i+r.  (*.)iy  ,,/)(£. )+r.  (£.)r,.  ,^(je.)r.  (*.)r,.  .*,(£.)  + 
jr  J  G"1)*  J  Jr  J  G"1)*  J  jr'  r  G-1)*'  r 


(32) 


■IV 


#  #  #  #  #  # 

r.  (*.)r,.  n-(#.)r.  (*.)ry  n,(*.)r.  (*.)r,.  ..,(*.)+ 
jr  J  G-U£  J  jr  J  G'1)*'  J  jr'  j'  G'1)*'  J 


(33) 


If  the  first  term  is  called  i,  the  second  ii,  etc.  ,  as  noted  above,  the 

» 


expression  for  U^.^(x) 


can  be  obtained.  Again, 


restricts  the  range  of  the  waves  to  |  1 

L  J"  J  J 

Consider  the  first  term  i  of  the  expansion  for 


j-1-  •  y 


» 


j  #  #  #  § 

U..jx)  .=r(x-5.  ,)r,.  Q. .  ,)r.  («.)u. ..(*.) 

jjk'  '  !  1  '  J-1'  (j-l)t'  j-17  '  j  j-r  jr'  y  jjk'  j 


#  #  xif  , 

i-r,  ,  r  (*.  , 

(j-  1)£V  J-  1  jr  J-  1 

> 


-  1 


+ 


R(x"Vi)ro-iK(Vi)VVi)  ■ 

< 

The  last  quantity,  U..^(£.  ^),  represents  an  incident  wave  traveling  to  the 

left  at  the  point  x  =  C.  .,  i.  e.  ,  the  left  end  of  the  span.  If  the 
#  J 

quantity  l)f^j  1^  ('^^ich  is  a  function  of  the  configuration  parameters 

of  all  the  supports  located  to  the  left  of  x  -  j  +0  and  also  of  the 

associated  shaft  characteristics)  signifies  that  this  wave  has  been 
reflected  at  x  =  i.  j,  then  the  reflected  wave  now  traveling  to  the  right 

^  if 

has  propagated  a  distance  (x  -  .  .)  as  indicated  by  R(x-£.  .).  Simila  * 

J  “  *  J  ■  * 

the  last  quantity  of  the  first  term,  represents  an  incident 

which  is  traveling  to  the  right  and  is  located  at  x  =  i.  e.  ,  the  right  >  a 
th  $ 

support  in  the  j  span.  The  quantity  IT\ r(^j)»  a  function  of  the  configura¬ 
tion  parameters  of  all  those  supports  located  to  the  right  of  x  =  jf.  -  0 
and  the  associated  shaft  characteristics,  signifies  reflection  of  this  wave 
to  the  left  at  x  =  The  .ext  term,  R(f  -f  j),  signifies  propagation 

of  the  reflected  wave  to  the  left  through  a  distance  ( S..-H .  j),  i.  e.  ,  length 
t  h 

of  the  j  span,  until  the  wave  reaches  x  =  f.  ..  The  wave  is  now 
reflected  again  to  the  right  as  indicated  by  Tq  l)jp^j  j)’  anc*  R(X_J^-  j) 
signifies  propagation  through  a  distance  (x-H.  .)  to  the  point  under 

»  J "  1 

investigation.  Thus,  the  term  U^^(x)  .  represents  the  contributions  of 

two  reflected  waves  to  the  right  immediately  following  the  initiation  of  the 
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I 


> 


incident  waves,  U_^(x)  and  U_^(x).  This  is  illustrated  in  Figure  10. 
Consider  the  second  term,  ii,  of  the  expansion  for 


a  # 


# 


i-r,.  I)c(t  ,)r.  (t  .) 

L  U-  1)*'  J-  1  jr'  j-  r 


» 


tJ...  (x) 
Jjk' 


ii 


,)rjr(V  i4-  ,)(<V  .>*<V V>> 


x  VV  VV 4 

R(x‘ 1  j  - 1  )r(j  - 1  )f  ( 1  j  -  1  )rj  r({  j  -  1  )r(j  -  1  )t  (£j  -  I  )Cjjk(  V  1 ) 


which  may  also  be  written  as  follows  (see  Eqs.  (22)): 


» 


#  #  #  #  # 

..-R(x-f.  ,)r,.  iU(jC.  , )r(jC ,)r.  (t)R(je.-*.  ,) 

u  j-F  (j-l)i'  J-F  J  J-F  jr*j  j  j-F 


#  #  #  > 

x T(.  nf(C.  ,)R(i.-e.  ,)r.  (i.)u...(«.)  + 

(j-l)f'  J-F  J  j-F  jr'  y  Jjk'  y 


< 


R(x-  cj.  i)r(j.  i)t j)R(f  j-*j-  l)rjr(VR(V  V  1  )r(j-  l)^£j-  l)t7jjk(V 


» 


If  the  same  reasoning  is  used  for  the  analysis  of 


the  added 


terms  correspond  to  a  wave  propagating  and  reflecting  an  additional  two 
more  times,  with  the  final  reflected  wave  again  traveling  to  the  right. 
This  is  shown  in  Figure  1 1.  Similar  reasoning  may  be  applied  to  the 

» 

remaining  terms.  These  remaining  terms  in  the  series,  i.  e.  ,  U..,(x) 

»  I  Vi  . 


» 


V 


etc.,  together  with  tL^fx) 


» 

.  and  tJ. ..  (x) 

i  jjk' 


in 


ii 


account  for  all 


of  the  reflected  waves  traveling  to  the  right.  The  summation  of  all  of  the 
waves  traveling  to  the  right  is  complete  if  the  incident  wave  is  added  to  the 
above;  i.  e.  , 


»> 


» 


» 


V(x)'njik(x)1Vx)ii+ujji<(x) 


» 


» 


=ujjkW+ujjk(x) 


..+0...  (x) 
u  jjk' 
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.  .  .  + 

111 


(34) 


Figure  10.  First-term  Propagation  of  Reflected  Waves 
in  the  j"/ j - j •  j€ j]  Portion  of  the  Shaft. 
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P(ojk) 


Figure  II.  First*  and  Second* term  Propagation  of  Waves 
in  the  [J'j-T  l  j]  Portion  of  the  Shaft. 
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Similar  results  may  be  obtained  for  ^^(x).  That  is, 


« 


« 


tijjk(x)-Ujjk(x) 


.V 


« 


..+c...  (x) 
11  jjk' 


. .  .  . 

111 


(35) 


« 

The  summation  of  tJ...  (x),  the  reflected  wave  traveling  to  the  left  in  the 
<  JJk 

j**1  span,  with  Uj^(x),  the  incident  wave  traveling  to  the  left,  accounts 

for  all  the  waves  traveling  to  the  left  as  they  are  observed  passing  point 
x  in  the  jt^1  span  of  the  shaft. 


The  complete  response  on 


£  £ 

j-1’  j 


at  point  x  can  be  determined  by 


properly  combining  all  of  the  waves  traveling  past  point  x  in  both  the 
left  and  right  directions  (see  Eq.  (127)  in  Appendix  F).  That  is, 


On 


£.  1. 

j- 1  j 


Y 

F 


1  * 

—  c 

s  y 

# 

0 


# 

0 

* 

c, 


# 

c 

# 

c 


1 1 
21 


# 

c 

if 

c 


12 

22 


1 

r 

>» 

V* 

l  v*>J 

or 


afa  »>  a  «< 


# 

F(x)=C. 


»> 


«< 


C21tJjjk^x*+C22tIjjk^x* 


(36) 


Eq.  (36)  indicates  a  relatively  direct  manner  in  which  the  waves  are 
combined  to  obtain  the  total  response  at  point  x. 


By  examining  the  portion  of  the  shaft  on 

_  »>  «< 


£.,  £.+1 

LJ  j+l. 


in  terms  of  a  single 
driving  force  Pfa^),  U^  +  1^k(x)  and  U^.  +1)jk(x)  may  be  thought  of  as 


"modified  total  waves"  which  are  initiated  by  the  driving  force  at  x  =  a.,  . 

>»  «<  pi  Jk 

I*.  £. 

j-1  j 

x  =  £.,  the  waves  are  modified  or  refracted  by  the  j  n  support  and  pass 
J 


As  soon  as  the  total  waves,  U„k(x)  and  tJ_k(x),  on 


travel  to 


through  as  indicated  by  ^(j  +  i)jk(x)  anc*  ^(j  +  l)jk^  traveling  in  the 


(j  +  1)^  span  (see  Eqs.  (138)  and  (139)  in  Appendix  F).  Thus, 
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>:> 


#  #_j  #  if 


u(j«jjkW=R(-<J,rcnci2r0ii)r(‘j). 


- 1  »>  #  .  #  «< 
u...  (*.)+c:  (1.) 

L  JJk'  J  11  12  JJk  J 


(37) 


«< 


Vu,kW-“<Vi-x,^«)r(1j+i,R(VrV 


#  #  ,#  # 

.1+cnci2r(,«)r<V. 


-  1 


r»> 


#  .  a  <« 

V‘j)+Cii  c«wv 


(38) 


where,  as  before,  the  first  subscript,  j  +  1,  indicates  that  the  wave  is 

t  h 

traveling  in  the  (j  +  l)  span;  the  second  subscript,  j,  indicates  that 

the  original  wave  originates  in  the  span  where  the  driving  force  is 
located;  the  third  subscript,  k,  indicates  the  exact  location  of  the 
driving  force. 

An  explanation  of  traveling  wave  behavior  similar  to  the  previous 

»>  <« 

discussion  may  be  applied  to  ^(j  +  l)jk^X^  and  +  in  the  (j  +  l) 

span.  Note  that  the  second  term  of  ^(j  +  l)jk^  can  be  restated  as  an 
infinite  series: 


th 


-  1 


i  r _ n  . 

V  #  ,#  a 


in 


#  ,#  # 


=  H-c-1Jcu)r(.M)r(t.H-c;;c12,r(.+I)r(CK-c;Jcu)r(.+I)r(C)  + 


IV 


#_1#  # 


#  j  # 


*  !#  # 


^~C1  lC12^r(j  +  l)r^j^_  C1 1C 1  +  C1 1C1  2^r(j  +  l)r^j^  + 


(39) 


■A  ga  th  irder  of  ^e  terms  is  indicated  by  i,  ii,  iii,  etc.  Considering 
,  i,  -nd  taking,  for  the  time  being,  only  the  first  term  i  of  the 


•  b ^  e  ” t.u  db  Sion , 


»> 


n0<i)jkM 


r»> 


#  _ j  a  «< 


.*R<*-<AW+cnc«V<J, 
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I 


th 

where  the  last  term  signifies  that  the  total  wave  in  the  j  span  passes 

through  the  j**1  support  and  enters  into  the  (j  +  l)**1  span  as  a  modified 

wave.  R(x-£.)  signifies  that  the  modified  wave  has  propagated  to  the 
J 

right  a  distance  (x-£.).  This  is  illustrated  in  Figure  12. 

J 

Consider  now  the  second  term,  ii,  ol  the  expression  for 

r»  #  ,#  #  1-1 

I+cncurG+i)r<V 


»> 


0(j  +  l)jk(x) 


#  #  .#  # 
ii=R(x-fj)(-c11c12)r(j+1)r(i) 


>» 


#  ,  #  «< 
u...  (£.)+c;  c1?u..,(£.) 

jjk'  y  11  12  jjk  j'. 


which  may  be  rewritten  as  follows  (see  Eqs.  (22)): 

»>  #  #  ,#  #  #  # 

U(j+i)jk^x^  ii=R(x-fj)(-cnci2)R^j+rVr(j+i)r(Vi)R(Vi'V 


»> 


#_j#  «< 


LW+ciicuV<J> 


If  the  same  reasoning  as  before  is  followed, the  added  terms  correspond 
to  the  wave  propagating  to  the  right  after  being  reflected  twice,  first  at 

the  (j  +  l)^  support  and  then  at  the  j*1*1  support.  Appendix  G  shows  that  the 

#_1# 

term  ("C]iC]2^  corresponds  to  a  reflection  matrix  of  a  fixed  end  support. 

This  means  that  the  j**1  support  acts  as  a  one-way  fixed  support  which 
completely  reflects  modified  waves  coming  from  the  right  after  being 

reflected  back  from  the  (j  +  l)**1  support,  but  which  permits  waves 
traveling  from  the  left  to  pass  through.  This  is  illustrated  in  Figure  13. 
The  same  analysis  may  be  applied  to  the  remaining  terms  of  the  series. 

«< 

Also,  similar  reasoning  may  be  applied  to  tj(j  +  i )jk(x)*  the  modified  wave 
traveling  to  the  left  as  indicated  in  Eq.  (38). 


The  complete  response  at  point  x  on 


£ 


j*  V 


due  to  P(ajk)  can  be 


determined  by  properly  combining  all  of  the  waves  traveling  past  point  x 
in  both  the  left  and  right  directions  (see  Eq.  (140)  in  Appendix  F).  That 
is, 


46 


P(Ojk) 


Total  wav*  traveling  to 
the  right  In  the  j til 
span 


(x) 


I 


Figure  12.  First"  term  Propagation  of  Modified  Total  Wave 
Traveling  to  the  Right  in  the  Tj  + 

Portion  of  the  Shaft.  L  J 
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Figure  13.  Second-term  Propagation  of  Modified  Total  Wave 
Traveling  to  the  Right  in  the 
Portion  of  the  Shaft.  L  J 


on 


'■‘V 

Ic' 

s  y 

-1 

# 

1 

# 

»> 

Y 

0 

C11 

C12 

D0+i)jk<x) 

# 

# 

t 

# 

«< 

F 

0 

crj 

C21 

C22 

J 

%Djk<x) 

or 


Y(x)= 


"#  »>  #  «< 

C 1 1  U(j  + 1  )jk(x)+C  1 2U(j  + 1  )jk(x) 


F(x)= 


#  r#  »>  #  «< 

Cf  C21U(j  +  l)jk(x)+C22U(j  +  l)jkW 


>  (40) 


The  same  manipulation  of  traveling  wave  concept  used  for  the 
explanation  of  the  dynamic  response  in  the  (j  +  l)™  span  may  be  extended  to 
explain  the  dynamic  response  in  the  (j  +2)t^1  span  as  well;  it  also  can  be 

th 

extended  to  explain  the  response  in  the  (j+3)  span,  and  so  forth,  until  the 

n**1  span  is  reached.  Similar  argument  also  can  be  used  for  the  dynamic 

response  in  the  (j- I)***  span,  in  which  the  modified  waves  (see  Eqs.  (135) 
and  (136)  in  Appendix  F)  are  as  follows: 


>>>  #  #  # 

u0-i)jk(x)=R(x_V2)r(j-2)je(V2)R(£j-i 


V2> 


'#  # 
I+C 


-1 

11 


#  # 

C12I(j-2)Je^j-l). 


-1 


#  #  >>>  <<< 

Lciici2ujjk(Vi)+V(lj-i). 


«< 


u0.l)jk(x)=R(C.i-x) 


'#  #  .#  # 

I+cilcl2r0-2)£({M)J 


-1 


X 


#  #  >>>  <<< 

Lc'uCizVWm1 


(41) 


(42) 


The  term 


J+Cn*i24Wi) 


infinite  series  as 


1 


can  be  restated  in  terms  of  an 
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■##,## 

1+cnci2ro-2)je(£j-i). 


-i 


it  if  §  if  if  if  if  if  if  if 

■H-c-i;c12)ijj.2){(tj.1H-ci!ci2)r0.2,t(«j.1)(-c-1Jc12)r{..2)4(Jj.1)  + 


*  i.#  • 


#  1#  if 


#  i#  # 


^CllC12^I(j-2)^£j-l^’CnCl2^I(j-2)4^j-l)(“CnC12)I(j-2)l^j-l^  + 

(43) 

Again,  the  order  of  the  terms  is  indicated  by  i,  ii,  iii,  etc.  Considering 

»> 

U^.  +  j)j^(x)  and  taking,  for  the  time  being,  only  the  first  term  i  of  the 
above  expression. 


»> 


VujkM 


#  ' 


i‘R(x' -  2)rG -  2)J?( V  2)R ( V  J  '  V  2) 


<*< 


r#  .#  >» 

cnc12V(Vi)+V(Vi) 


where  the  last  term  signifies  that  the  total  wave  in  the  j**1  span  passes 

th  til 

through  the  (j-l)  support  and  enters  into  the  (j-2)  span  with  modified 

if 

wave.  R(f .  j-f.  ^ )  signifies  that  the  modified  wave  has  propagated  a 

J  J 

distance  (H.  ,)  and  then  is  reflected  to  the  right  by  the  (j-2)**1  support, 

J  3  if  if 

as  indicated  by  the  reflection  matrix  2)f(£j  2)’  R(x_J^j  2)  a8a*n 

signifies  that  the  reflected  modified  wave  has  propagated  a  distance 

(x-4.  P),  and  then  the  wave  is  at  the  point  under  investigation.  This  is 
}~c 

illustrated  in  Figure  14. 


Consider  now  the  second  term,  ii,  of  Eq.  (41),  corresponding  to  the  ii 
term  of  the  expansion  of 


50 


Figure  14.  First-term  Propagation  of  Modified  Total  Wove 
Traveling  to  the  Right  in  the  j-2'  ^  j-|] 
Portion  of  the  Shaft.  *■  ^ 
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#  #  J#  # 


I  +  CllC12r(j-2)£^j-l^- 


-  1 


>» 


U(j-l)jkW 


11 


-r(x-?j  2)r^j  2)2(i  2)r(^j  j-fj  2)(-cHci2)r(i 


#_j#  # 


I  ”  “  T  T 


«< 


[cnc.2Vfj-.)+Wi>J 


which  may  be  rewritten  (see  Eq.  (22))  as  follows: 


»> 


0(i-i)jkw 


ii=R(x-^_2)I(j -2)je(fj - 2>R(£j-  r£j-2)(‘CllC12)R(V1"^j-2) 


Xr(j-2)^(ij-2)R(£j-r  fj-z> 


#-l#  »> 


«< 


cnci2V(lj-*)+V*J-i) 


If  the  same  reasoning  as  before  is  used,  the  added  terms  correspond  to 
the  wave  propagating  to  the  right  after  being  reflected  two  more  times. 
Again,  the  (j-l)**1  support  acts  as  a  one-way  fixed  support.  This  is 
illustrated  in  Figure  15.  The  same  analysis  may  be  applied  to  the 
remaining  terms  of  the  series.  Similar  reasoning  may  be  applied  also 
«< 

to  Uq  jjj^(x),  the  total  wave  traveling  to  the  right. 

The  complete  response  at  point  x  on 


B  B 

J-2’  j-1 


due  to  P(a^)  can 


be  determined  by  properly  combining  all  of  the  waves  traveling  past 
point  x  in  both  the  left  and  right  directions  (see  Eq.  (137)  in  Appendix  F). 
That  is, 
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B  B 

j-2’  j-1 


“  “1 

[  1  # 
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‘  # 

#  1 

>» 
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C  1  2 

°0-HjkW 
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«< 
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C2I 

C„ 

22  J 

Y(x)  = 
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s  y 

#  r# 


»> 


#  <« 


LC1  ig(j-  l)jk(x)+C12U(j-  i)jk(x)J 

F(x)=c^C2itJ(j_  1)jk(x)+C22U(j .  1)jR(x) 

52 


(44) 


Figure  15.  Second-term  Propagation  of  Modified  Total 
Wave  Traveling  to  the  Right  in  the  .. 2  ^i-|l 
Portion  of  the  Shaft.  ’  ** 
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The  same  manipulation  of  traveling  wave  concept  used  for  the  explanation 

of  the  dynamic  response  in  the  (j-l)**1  span  may  be  extended  to  explain  the 

dynamic  response  in  the  (j-2)^^1  span  as  well;  it  also  can  be  extended  to 

explain  the  response  in  the  (j-3)^  span,  and  so  forth,  until  the  1  8*  span  is 
reached. 

The  above  analysis  is  based  on  the  assumption  that  a  single  driving  force, 

P(allr).  acts  on  the  shaft.  For  driving  forces  distributed  along  the  shaft, 

the  traveling  wave  solution  for  each  individual  driving  force  is  combined 
to  give  the  total  response  at  any  point  x  by  direct  superposition.  Com¬ 
plete  solutions  are  given  in  Chapter  2  under  "Solution  in  Wave  Form". 


IMPEDANCE  MATCHING 


The  dynamic  response  of  the  shaft  system  shown  in  Figure  4  in  terms  of 


distributed  driving  forces,  P(a.,  ),  j=  1,  2 . n,  k  =  1,  2,  .  . 

is  expressed  by  Eq.  (27)  in  Chapter  2.  In  the  ith  span,  or  on  [f . 
i  =  1,  2,  .  .  .  ,  n  1_ 
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■# 
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^(x) 

.  k(j) 
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By  proper  control  of  the  conditions  at  both  ends  and  intermediate  supports, 
the  reflected  waves  and  modified  waves  may  be  eliminated,  or  at  least 
minimized,  such  that  the  dynamic  response  consists  solely,  or 
principally,  of  incident  waves  which  are  independent  of  support  config¬ 
urations.  The  mathematical  representation  of  the  above  statement  is  as 
follows: 

>>>  >  «<  < 

tJ.(x)  and  U.(x)=  tJ.(x) 


or 

»  « 

=  fr(x)=  0  (45) 


where  i  =  1 ,  2,  .  .  .  ,  n. 

It  can  be  seen  from  examination  of  the  reflected  and  modified  wave  terms 
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»  « 

in  Eq.  (27)  that  the  condition  U.(x)  =  U.(x)=  0,  i  =  1 ,  2,  .  .  .  ,  n,  is 
satisfied  by  setting 


#  #  # 

r0(°)=rit(fi>=r2£(C2); 


#  # 

r(n-  1  )f  ^n-  ° 


and 

#  # 

r  (£  )=r,  n  (* 

nx  n7  (n-  1  )r' 


n- 1 


# 

)=r, 


(n-  2)r^n-  2 


)  = 


# 

=r 


(46) 


Physically,  this  is  equivalent  to  setting  the  impedance  at  each  support 
looking  both  to  the  left  and  to  the  right  cf  the  support  point  equal  to  the 
shaft1  s  characteristic  impedance.  This  manipulation  is  analogous  to 
load  matching  in  electric  transmission  line  theory,  in  which  a  line  is 
terminated  by  a  load  equal  to  its  characteristic  impedance;  i.e.,  there 
is  no  reflected  voltage  wave.  The  proper  ratio  of  input  to  output  voltage 
in  the  incident  wave  is  satisfied  at  the  boundary,  such  that  there  is  no 
need  for  the  presence  of  a  reflected  wave  to  satisfy  the  boundary 
conditions.  The  existence  of  reflection  waves  stems  from  a  need  for 
satisfying  boundary  conditions.  A  mathematical  verification  of  the  above 
statement  may  be  obtained  from  Eqg.  (23)  and  (24),  which  are  expressions 
of  impedance  in  terms  of  reflection  matrices  and  shaft  characteristic 
impedance  as  shown  in  Chapter  2  on  pages  26  and  27  .  These  relation¬ 
ships  show  that  by  letting  all  reflection  matrices  vanish,  i.  e.  ,  conditions 
of  Eqs.  (46),  the  following  set  of  equations  holds: 


4(n-l)£*C 


#  #  . 
„.i>=c2ic'n 


and 


I  (*  )4,  n  (£  ,)=!,  (*  ?)= 

n'  n7  (n-l)r'  n-17  (n-2)r'  n-27 


#  #  #  i 

4lr<fl>=C21Cn 


(47) 


M  H  1  # 

where  zg  =  ^^l^ll’  cluant^Y  zs  called  the  characteristic 

impedance  of  the  rotating  shaft.  In  this  case,  all  impedances  appearing 
in  Eq.  (47)  are  termed  "matching"  or  matched  impedances. 


If  Eqs.  (46)  or  Eqs.  (47)  are  satisfied  by  proper  control  of  all  support 
conditions,  the  solutions  of  the  dynamic  response  of  hypercritical  shafts 
in  wave  form  due  to  distributed  driving  forces  as  expressed  in  Eq.  (27) 
may  be  simplified  as  follows: 


In  the  i^ 


span, 


or  on 


=  1,  2, 


n 
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(48) 
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where 


>  r  > 

M=L  Vx) 

j=i 


n 


ui(x)=I,  Vx)  • 


Note  that  the  expression  for  U^j(x)  and  U..(x)  are  given  in  Eq.  (27), 


< 


> 


except  that  each  wave  function  U^.+  })j(x) . U.^(x),  and  U^.  j^.(x), 

DISCUSSIONS  ON  IMPEDANCE  MATCHING 


Examination  of  Eq.  (48)  in  the  preceding  section  permits  the  following 
observations  to  be  made.  If  the  impedances  looking  both  to  the  left  and 
to  the  right  at  eacn  support  are  matched  to  the  shaft' s  characteristic 
impedance  or,  in  other  words,  if  all  the  reflection  matrices  looking 
both  to  the  left  and  to  the  right  at  any  support  are  equal  to  null  matrices, 
the  dynamic  response  will  contain  only  the  incident  waves  due  to  driving 
forces.  This  means  that  incident  waves  simply  pass  through  all  supports 
without  any  modification.  This  is  consistent  with  the  mathematical 
explanation,  since  'he  incident  waves  are  independent,  of  support  conditions. 
Hence,  Eqs.  (46)  or  Eqs.  (47)  may  serve  as  a  criterion  for  the 
determination  of  the  support  conditions  for  shafts  which  are  subjected  to 
incident  waves  only.  Supports  which  satisfy  Eqs.  (46)  or  Eqs.  (47)  are 
termed  optimized,  since  the  dynamic  response  is  minimized;  no  amplitude 
buildup  due  to  reflected  waves  can  take  place. 


It  can  be  concluded, then,  from  the  above  observations  that  if  both  end 
impedances  are  m  itched  with  the  shaft  characteristic  impedance,  i.  e.  , 

r0(°)=rn(Cn)=  0 ,  or  z0(0)=in(£n)=C21cj  J  =  zs  ,  the  effects  of  intermediate 

supports  on  the  dynamic  response  of  the  rotating  shaft  are  redundant. 
However,  the  presence  of  intermediate  supports  does  raise  the  frequencies 
at  which  the  critical  speeds  occur,  including  ihe  fundamental  frequency. 
Actually,  one  could  place  a  predetermined  number  of  intermediate  supports 
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along  the  shaft  sufficient  to  raise  the  fundamental  frequency  above  the 
operating  frequency  of  the  shaft.  Obviously,  this  would  involve  the 
acceptance  of  a  p  malty  in  added  support  weight. 

Further  inspection  of  Eq.  (48)  shows  that  even  if  the  intermediate 
support  conditions  are  not  optimized,  their  effect  is  redundant.  Thus,  if 
end  impedances  can  be  matched  to  the  characteristic  impedance  of  the 
shaft,  nothing  is  gained  from  the  optimization  of  the  intermediate  supports 
in  terms  of  minimum  vibration  control.  Thus,  one  or  two  matched  end 
supports  should  more  than  adequately  suppress  any  excessive  vibration 
response  in  hypercritical  shafts  of  any  length.  The  impedance  values  for 
the  supports  of  the  shaft,  whether  supported  at  one  end  or  at  both  ends, 
are  the  same;  the  difference  lies  in  the  deflection  amplitudes,  where 
those  of  the  single- ended  shaft  are  twice  those  of  the  double-ended  shaft. 

The  matching  of  end  impedance  with  the  shaft  characteristic  impedance  is 
a  formidable  task,  not  only  because  impedances  are  a  function  of  frequency 
but  primarily  because  of  the  uniqueness  of  all  the  conditions  required  to 
control  these  impedances.  Moreover,  in  practical  applications,  the  choice 
in  end  support  configurations  is  limited,  since  they  usually  are  governed 
by  such  factors  as  transmission  gears,  couplings,  and  unwieldy  mountings. 
Thus,  the  intermediate  supports  have  to  be  employed  for  optimization 
purposes,  since  the  matched  end  impedances  are  not  available. 

If  all  supports  are  optimized  except  those  at  both  ends,  i.  e.  ,  if  all 

#  #  # 

reflection  matrices  are  equal  to  null  matrices  except  ^(0).  Fn(^n)^  0, 

the  dynamic  responses  due  to  distributed  driving  forces,  Eq.  (Z 7 ) ,  may  be 
expressed  in  the  same  form  as  Eq.  (48)  except  for  those  pertaining  to  the 


1st  and  n^  spans.  Those  responses  may  be  listed  as  follows: 


«<  "  «< 

Uj(x)=  2_,  ^ijW 

j=i 
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»> 


«< 


The  expressions  for  Ujj(x)  and  Uj^(x)  can  be  obtained  from  Eq.  (27)  by 
letting  i  be  1  for  j=i,  and  then  letting  j  =  n,  n-  1 ,  .  .  .  ,  2  for  j>  1.  It 


« 


should  be  noted  that  the  expression  for  the  reflected  wave  function  U,  ,  (x) 

#  >  > 

=0,  and  that  the  terms  Uq  •  •  •*  ^^(x)  vanish  also. 

On  j,  or  the  n**1  span, 


F 

where 

»> 


— c 

s  y 

# 

0 


# 

0 

# 

c. 


c, 


# 

c 


11 


21 


# 

c 

# 

c 


12 


22 


“1 

»> 

U  (x) 
n'  ' 

«< 

L  J 

(50) 


n 


>» 


u  (x) = y  u  .(x) 

n'  '  Z—  nj'  ' 

j=l 


«< 


n 


u  (x)  =  Y  ^  (x)  • 

n'  Z_,  nj' 

j  =  l 


»> 


«< 


In  this  span,  the  expression  for  U  .(x)  and  FT  .(x)  can  be  obtained  from 

Eq.  (27)  by  letting  i  be  n  for  j  =i,  and  then  letting  j  be  j  =1,  2,  .  .  .  , 

»  #  < 

n- 1  for  j  <n;  also,  the  term  TT  =0,  and  each  term  U/.,,%.(x) . 

<  ^  nn  (j  +  ljj 

Vi)j(x,=0- 

It  may  be  observed  from  the  above  expressions  that  only  incident  waves 
can  exist  in  every  interior  span  except  for  the  two  end  spans,  in  which 
incident  waves  originated  by  driving  forces  are  partially  reflected  at  the 

»>  «<  »>  «< 

end  supports  as  indicated  by  the  U^(x),  ffjj(x),  FJ^fx),  and  U^(x)  terms. 

The  existence  of  partially  reflected  waves  in  both  end  spans  is  an 
unavoidable  situation,  since  both  left  and  right  end  supports  are  not 
optimized,  i.  e.  ,  are  not  matched  with  the  shaft  characteristic  impedance. 
Thu6,  incident  waves  traveling  from  the  interior  spans  to  x=£j  and 

x=^n  j  simply  pass  onto  £(),  £  j  j  and  j,  £nJ  ,  respectively,  never  to 

return  to  £  £j,  $n  jJ  again.  These  waves  will  then  propagate  back  and 

forth  in  the  ^0,  and  j,  intervals,  to  be  damped  out  eventually 
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by  the  unmatched  impedance  at  x  =  0-0  or  at  x  =  £^+0,  respectively.  It 

should  be  remembered  from  the  first  section  of  this  chapter  "Traveling 
Wave  Concept"  that  for  the  reflected  waves  in  the 


and 


intervals,  the  supports  at  x  =  and  x  =  H  j 


act  as  fixed  supports  which  will  reflect  completely  all  incoming  reflected 
waves,  respectively,  from  the  left  and  right.  Thus,  since  reflected 
waves  and  modified  waves  vanish  along  the  shaft,  except  for  both  end 


spans,  the  shaft  system  is  optimized  on 


and 


£ 


n-r 


but  not  on 


The  minimization  of  vibration  response  in  a  shaft  having  end  supports 
that  do  not  lend  themselves  readily  to  optimization  may  be  accomplished 

by  letting  the  1st  and  the  (n-1)^^1  intermediate  supports  approach, 
respectively,  the  left  and  right  ends  of  the  shaft  as  closely  as  possible. 
In  this  case,  the  reflected  waves  will  be  restricted  to  very  short  end 
spans  and,  if  all  intermediate  supports  are  matched  with  the  shaft' s 


impedance,  only  incident  waves  will  exist  in  the  £j,  jj  portion  of 

the  shaft.  Mathematically,  the  above  statement  is  indicated  in  Eqs.  (49) 

»>  «< 

and  (50),  where  as  £.  and  (£  -  £  .)  approach  zero,  U.  .(x),  tJ,  .(x), 

»>  «<  1  n  n- 1  ij  ij 

U  .(x),  and  tJ  .(x)  also  approach  zero.  It  can  be  said,  then,  that  the 
nj'  nj'  '  rr 

closer  the  matched  l6t  and  (n-l)^  supports  are  placed  to  the  left  and 
right  end  supports,  respectively,  the  more  effective  the  amplitude 

suppression  will  be  in  the  end  portions,  j^O,  £jJ  and  j*  ^nJ  •  of  the 

shaft. 


It  has  been  observed  earlier  that  matched  or  unmatched  interior  supports 
would  be  redundant  if  the  end  support  impedances  are  matched  to  the 
characteristic  impedance  of  the  shaft.  By  the  same  reasoning,  inspection 
of  Eq.  (48)  should  show  also  that  for  the  case  in  which  the  two  outermost 
interior  supports  are  matched,  all  supports  interior  to  these  two  supports 
are  redundant.  From  all  this  discussion,  the  conclusion  may  be  made 
that  in  the  case  for  which  the  end  supports  are  not  available  for  optimization, 
one  or  two  matched  intermediate  supports  placed  very  close  to  the  end 
supports  should  more  than  adequately  suppress  excessive  vibration  response 
for  shafts  of  any  length.  The  impedances  of  the  interior  supports  should 
be  the  same  for  the  shaft  with  one  interior  support  placed  close  to  one 
end  and  the  shaft  with  two  interior  supports  placed  close  to  the  ends.  The 
only  difference  in  behavior  between  these  two  shaft  systems  exists  in  the 
amplitude  response.  It  is  larger  in  the  case  of  the  shaft  with  one  interior 
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support;  how  much  larger  this  value  of  the  amplitude  response  will  be 
depends  on  the  amount  of  wave  reflection  that  will  take  place  at  the  far 
end  support.  In  the  case  of  the  shaft  with  two  interior  supports  placed 
close  to  the  end  supports,  no  such  reflection  of  incident  waves  will  take 
place  at  either  end,  and  the  resulting  deflecting  amplitude  response 
should  be  less  than  that  for  the  single  interior  support  case. 


DETERMINATION  OF  THE  MATCHED  IMPEDANCES  OF  INTERIOR 
SUPPORTS 

# 

The  optimum  impedances,  Zj(iL),  j  =  1,  2,  .  .  .  ,  n- 1 ,  of  interior 
supports  are  determined  from  the  following  conditions; 


#  # 

rjC(^)=0,  j  =  1,  2,  .  .  .  ,  n-J 


and 


(51) 


#  # 

r.(U=0,  j  =  n-  1,  n-2 . 1  . 

J  J 

It  may  be  shown  as  before  (see  Eqs.  (47))  that  the  above  expressions  are 
equivalent  to  setting 


Vi)=c2icn= 


#  #  ,  # 
z 


and 


#  #_j  # 


z.  (d  .)=C.. .  C.  ,  =  z 
jr'  y  21  11 


j=l,  2,  .  .  .  ,  n- 1 


j  ("i"*li  n_  2,  .  .  .  ,  1 


(52) 


A  fr 

The  details  of  z^  (V  and  z^(^)  are  summarized  in  Eiqs.  (17)  through 
(24). 


From  these  equations,  a  set  of  supports  can  be  uniquely  determined  such 
that  the  shaft  will  have  minimum  dynamic  response  characteristics. 
However,  it  should  be  noted  that  the  shaft1  s  characteristic  impedance  is 
a  function  of  frequency,  w  ;  hence,  all  matched  support  impedances  must 
necessarily  be  functions  of  frequency.  To  provide  these  supports,  in 
reality,  with  matched  impedances  is  a  formidable  task,  not  only  because 
they  are  a  function  of  frequency,  but  primarily,  as  has  been  stated  before 
in  the  preceding  section,  because  of  the  uniqueness  of  the  requirements  for 
these  impedances.  Each  support  must  be  a  translational  and  rotational 
mass- spring-damper  unit  with  a  built-in  frequency  dependency  to  meet 
these  requirements.  This  is  a  physical  situation  that  is  not  easily  attained. 
Hence,  from  a  practical  point  of  view,  the  possibility  of  imposing  weaker 
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conditions  on  support  requirements  for  minimum  vibration  response  of 
hypercritical  shafts  should  be  investigated. 

QUASI- MATCHING 

A  compromise  approach,  "quasi- matching",  to  minimizing  critical  speed 
vibration  may  be  realized  in  the  form  of  weaker  conditions  on  impedances 
at  support  locations.  By  examining  the  expressions  of  reflected  and 
modified  waves,  some  interesting  observations  may  be  made.  First,  by 

> 

writing  out  the  expression  for  the  incident  wave  U„(x)  in  Eq.  (27)  in 
complete  matrix  notation, 


k(i) 

iiM=  I 

k=l 


-iejVw  (x-aik) 


-  e 


,VT  (*-aik) 


#  #  .=* 
c+c;Ip(aik) 


The  term 

-i (x-a.,  )  _ 

e  1  =cos(eiVw  (x-a.k))-i  sin  (ejVw  (x-c.k)) 

is  a  complex  entry  which  affects  only  the  phase  relation  of  the  wave  as  the 

-e2V^  (x-aik) 

argument  (x-a^k)  increases.  The  other  term,  e 


is  a 


function  for  which  the  numerical  value  decays  exponentially  as  (x_ajk) 

-e2Y^(x-aik) 

increases.  Generally,  this  term. _ e  ,  is  very  small  when 

(x"aik) 

compared  to  the  first  term,  e  ,  which  does  not  decay.  Hence, 

it  may  be  neglected  with  very  little  effect  on  the  solution.  Second,  by 

examining  the  expression  of  reflected  wave  U..(x)  in  the  i  span  as  given 


11 


« 


by  Eq.  (27),  the  term,  which  corresponds  to  !J..(x) 


ir 


(see  Eqs.  (33) 


and  (35)),  can  also  be  written  out  in  complete  matrix  notation: 
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#  #  > 

'  i  7  ir'  i7  11'  i7 


k(i) 

■  E 

k=i 


-iejfiTfi.-a.  )  -ierfT(l-x) 

(e  1  1  lkr.  (1.)..)®  1  1 

'  ir' 

-ie.Vw^.-a  )  -e2V^(^-x) 

(e  1  1  lkr.r(£.)21)e  2 

-e7~\/Z  (l.-a...  )  -i e.VcE  (£.-x) 

(e  2  1  lkrir(«i>i2)e  1 


( 


-e^(Vaik) 


r.  (£.)„)e 

ir'  i7227 


-e2V^  (*fx) 


u  u  ~ 
c+cf  lp(aik> 


In  the  reflection  matrix  I\  (£.),  the  only  element  not  modified  by  an 


ir'  l 


exponential  decaying  function  is  F.  (£.)jj.  Thus,  for  cases  in  which 

these  elements  become  negligible  with  increasing  w  ,  only  the  element 
which  is  not  modified  by  an  exponential  decaying  function  need  be 
considered.  This  leads  to  the  simpler  expression: 


R(«i-x)rir(fii)0.i(£i) 


k(i) 

E 


-ie  V^C-a.J  -i«,VT (C-x) 

(e  Fir  i  1 1  e  0 


#  t  .=* 

C+Cf  P<aik> 


The  same  manipulation  may  be  applied  to  all  the  other  reflected  and 
modified  waves.  Hence,  as  a  logical  conclusion,  the  shaft  matching 
conditions  as  stated  in  Eqs.  (51)  can  be  replaced  by  weaker  ones  as 
follows: 

ri£(£l>ll=r2£<Vir*  •  '  =  r(n-l)£^n-l^ll  =  ° 
and 

r,  ,%  (l  i)n=r/  (*  J,  ,=  ...=  T,  (£.)=0 

(n-l)r'  n-l'll  (n-2)r'  n-2'11  Ir'  l7 


r  53) 


where  the  subscript  "11"  indicates  the  element  of  the  first  row  and  the  first 
column  of  each  reflection  matrix.  If  all  interior  supports  are  designed 
such  that  the  above  conditions  are  satisfied,  then  the  shaft  is  termed 
quasi- matched.  By  examining  Eqs.  (23)  and  (24)  and  their  related 
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expressions,  it  ia  evident  that  under  the  condition  of  quasi- matching,  each 


A 

term  z 


£ i ) »  •  •  ••  ^(n-l)r^n-l^’  ^(n-2)r^n-2^* 

^  #  W  ff  _  j 

.  .  .,  | ),  approaches  zg  =  C^jCj^Eqs.  (47)),  which  relates  to  shaft 


characteristic  impedance  (i.  e,  ,  matched  conditions)  as 

lik 


or  expf-e^V^” (4.-x)),  i  =  1,  2,  .  .  .  ,  n,  tends  to  zero,  where  a.,  or  x 


is  on 


[°-  ‘n]  • 


The  application  of  the  criteria  given  in  Eqs.  (53)  for  the  quasi- matching 
of  interior  support  impedances  to  the  characteristic  impedance  of  the 
shaft  to  the  design  of  several  support  configurations  will  be  illustrated 
in  Chapter  4.  For  convenience  sake,  a  shaft  with  only  one  interior  support 
and  two  end  supports  will  be  studied. 
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{ 

i 

i 

i 


i 


expressions,  it  is  evident  that  under  the  condition  of  quasi- matching,  each 

term  z2£^2^ . z^-lj^n-P’  Z(n- l)r^n- 1^'  ^(n- 2)r^n- 2^’ 

.  .  approaches  zg  =  C^iCjjfEqs.  (47)),  which  relates  to  shaft 

characteristic  impedance  (i.  e.  ,  matched  conditions)  as  exp(- 
or  expf-e^V^ (f.-x)),  i  =  1,  2,  .  .  .  ,  n,  tends  to  zero,  where  a^  or  x 
is  on  £o,  J 


is  on 


The  application  of  the  criteria  given  in  Eqs.  (53)  for  the  quasi-matching 
of  interior  support  impedances  to  the  characteristic  impedance  of  the 
shaft  to  the  design  of  several  support  configurations  will  be  illustrated 
in  Chapter  4.  For  convenience  sake,  a  shaft  with  only  one  interior  support 
and  two  end  supports  will  be  studied. 
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CHAPTER  IV 


TRANSMISSION  LINE  ANALOGY  SOLUTION  OF  THE 
SUPERCRITICAL  SHAFT  WITH  ONE  INTERMEDIATE 
SUPPORT- -USED  AS  AN  ILLUSTRATION  PROBLEM  FOR 

QUASI- MATCHING 
DESCRIPTION  OF  PROBLEM 


A.  uniform  rotating  shaft  with  one  interior  support  will  be  used  in  this 
chapter  for  the  development  of  equations  for  different  support  config¬ 
urations  in  terms  of  quasi- matched  conditions.  These  expressions 
should  be  useful  to  the  designer  in  that  impedance  values  termed  "best11 
for  the  supports  may  be  established  directly  without  an  excessive 
amount  of  computations. 


The  criterion  for  quasi- matching  of  the  interior  support  impedance  with 
the  shaft  impedance,  looking  to  the  left  and  right  of  x  =  f  j  is,  referring 


to  Eqs.  (53). 

ru^ih r° 

rir<Vir° 


#  # 


(54) 


In  general,  ^(C  j)L  ^  r(f  j  )j  j ,  or  1^(4  j^r^^j);  the  condition 

£  # 

exists  only  if  the  interior  support  is  placed  at  the  midpoint 
(£^  =  1/2  d  the  shaft  span,  and  both  end  supports  have  the  same 


configuration. 


A  physical  model  of  the  shaft  used  for  the  analysis  in  this  chapter  is 
shown  in  Figure  16.  This  corresponds  to  the  general  case  shown  in 
Figure  4  by  letting  n  =  2.  If  a  shaft  is  used  here  having  identical  end 
support  configurations, 


20(0).Z2(I2)  .  (55) 


By  direct  substitution  of  Eq.  (55)  into  the  expressions  of  f,»P  ,)  and 

# 

rir<V  and  assuming  that  ij=l/2  it  can  be  shown  that  (see  Appendix 


H) 
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Figure  16.  Uniform  Shaft  With  One  interior  Support  at  the 
Midpoint  of  Total  Span. 


(56) 


#  # 

Wl*  rir^l^ 

rU^lhl*rir^l^ll  ‘ 

The  conditions  of  Eqs.  (55)  and  (56)  will  be  used  in  the  development  of 
a  number  of  example  problems  of  this  chapter. 

REFLECTION  MATRICES  AND  SUPPORT  IMPEDANCES 


Reflection  Matrices  at  Supports  in  Terms  of  Impedances 

A  set  of  corresponding  equations  for  the  shaft  system  shown  in  Figure  16 
(V1/2 1  2 )  can  be  obtained  from  Chapter  2  simply  by  letting  n  =  2  in 

Eqs.  (15)-(22),  as  follows: 


# 

*  ,  * 

-1 

'#  A  * 

r0(°)= 

;cz2^0(°)c12j 

.c2rV°>cn 

# 

~  #  A  # 

-  1 

"#  # 

C 2 1  ’  Z 1  f  ^1^C11 

ri£<V= 

_C22+Zl£^l^C12 

# 

*  #  #  #  “ 

-  1 

'#  #  #  ■ 

rlr(ci>= 

'C22+zlr^l^C12 

C2l"zlr^l^Cll 

i 

"  #  #  #  " 

-1 

■#  «  i  • 

rz(i2)= 

"C22+Z2^2^C12 

c2,-z2(£2)cu 

where 


(57) 


#  #  #  f 

r0(<1)=R(£1)r0(o)R(^1) 

#  #  #  # 

r2(£i)=R(€2-Ci)rz(£2)R(£2-£i)  . 
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Impedances  in  Terms  of  Reflection  Matrices 


A  set  of  corresponding  equations  for  the  shaft  system  shown  in  Figure 
16  (J^l/2  £  can  be  obtained  from  Chapter  2  simply  by  letting  n  =  2 


in  Eqs.  (23)  and  (24),  as  follows: 

Onjo,  , 

A  r#  ##  ir#  ##  I-1 

z0£(x)xC2l  +  C22r0^x^  Cll  +  C12r0^X^ 

°n[^l,  «2]  , 

a  "#  #  ^  ip#  #  #  i-i 

zl^x)=  C21  +  C22ri<^x)  Cll  +  C12ri^X^ 


°n [ ^ 1 ’  ^2  ’ 


m  f#  #  #  ir  #  ##  i-i 

z2r^  =  C21  +  C22r2^x^  C1  l  +  C12r2^x^J 

On  [0,  £j  1  , 

*  r#  #  #  ir#  #  # 

zlr(x)=  C21  +  <^22rir^x^  Cll  +  C12rir^X^ 


SOLUTION  IN  WAVE  FORM 


From  Eq.  (27)  in  Chapter  2,  the  corresponding  solution  in  traveling 
wave  form  for  the  shaft  system  shown  in  Figure  16  (£  ^  l/2  £^)  can  be 


obtained  as  follows: 
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where 


»>  »>  >;> 

^1(x)=^11(x)+U12(x) 


«<  «<  «< 

u1(x)=un(x)+ui2(x) 


>»  »>  »> 

U2(x)=U2i(x)+U22(x) 


«<  «<  «< 

U2(x)=U21(x)+U22(x) 


»>  »  > 

U^x^U^xJ+U^x) 


> 


«<  «  < 

U1 1 (x)  =  u i  i(x)+un(x) 


>  kW#  #  #  * 

Uu(x)=  £  R(x-alk)C+Cf"1P(alk)[H(x-alk)-H(x-f1)] 


k=l 

>  #  #  « 

uu(x)=  2]  R(alk-x)C_C'1P(alk)H(alk-x) 

k=l 


<<  # 

Uu(x)=R(x) 


« 


"#  #  # 

-1# 

r 

i-r0(o)rlr(o) 

1  V°>1 

er: 

k) 

r#  #  # 

1 

-1 

-*) 

L^lr^lWJ 

#  # 

< 

1 

#  #  >_  <_ 

R(ll)rir(ll)Ull(ll)+Ull(°) 


lr 


Un(x)=R(lrx)[ 


69 


»>  »  > 

U22(x)=U22(x)+U22(x) 

<3C<  «  < 

U22^=U22^+U22^X) 

>  k(2)  0  |  |  ~ 

S22(x)=  £  R<x'a2k)C+Cf  lp(a2k)H(x'a2k) 
k=  1 


<  k(2)  §  §  §  ~ 

«22<x>=  I  R(a2k-x)C-C^P(a2k>[H(a2k'x)-H<trx>] 

k=l 


>>  ff  if  ff  if  -Iff 

u22(x)=R(x-i1)Li-rll(i1)r2(i1)J  rljg(f j) 

LR^2“^l^^l^U22^2^+U22^1^.  H^x“^ 


« 


ff  ff  if  ff 

U22^=R^2-x^_I-*2^2^l£  ^2^. 


-1# 

W 


»> 


>  #  #  < 
U22^2^+R^2"^l^l£^  1^U22^2^J  H(x_^ 

ff  # 


#  #  ,#  ff  -1  #_,#»>  «< 

U12(x)=R(x)ro(°)R(V|  I+Cl  1C12^0^1^  C11C12^22^1^+^22^1^ 


XH(fj-x) 


«< 


«< 


ff  ff  ff  ff  ff  - 1  ff  ff  »> 

U12(X)=R(£1'X) -I+Cl  1 C 1 2ro(^  1 )  -  -C11C12U22^1^+U22^1^ 


»>  ff 

^21  (X)=R(X“  ^  i ) 


■##.## 

I+CllC12r2^1^ 


-ir>»  ff  ,  ff  «< 

^n(^1)+c11c12^11(f1) 


H(£rx) 

H(x-fj) 


«<  ff  ff  ff 

tr21(x)=R(£2-x)r2(£2)R(£2-f1) 


#  #_,#  # 
I+CIlC12r2^1^ 


-  1 


The  preceding  Eqs.  (57),  (58),  and  (59)  give  the  pertinent  formulas 
required  for  the  optimization  of  support  conditions  of  the  two- span 
shaft  with  the  interior  support  at  position  H  j  and  two  general  support 

configurations  at  the  ends.  Also  from  Eqs.  (59),  the  dynamic  response, 
i.  e.  ,  the  induced  deflections  and  forces  due  to  distributed  driving  forces, 
of  the  shaft  under  investigation  can  be  calculated  directly. 

DETERMINATION  OF  THE  QUASI- MATCHED  IMPEDANCE  AT  THE 
INTERIOR  SUPPORT 

In  the  subsequent  set  of  example  problems  for  the  shaft  on  three  supports 
it  is  assumed,  for  convenience  sake,  that  both  end  supports  have  the 
same  configuration  so  that  Eq.  (55)  holds.  In  addition,  the  interior 
support  is  assumed  to  act  at  the  mid- span  (£.=  l/2  &?)  of  the  shaft  such 
that 

^lJ0^l)  =  ^lr^P’ 

or  in  terms  of  the  quasi-matching  condition  for  the  impedance  of  that 
support  to  the  characteristic  impedance  of  the  shaft, 

ru<li>ii"rir<*i>n=0  •  <60) 


t  »  » 

Since  ^(jP^s  Ij  (fj),  only  the  condition  rir(4j)=0  needs  to  be  considered 

in  the  following  analysis  to  determine  design  criteria  for  the  interior 
support  meeting  the  requirements  for  quasi-matching.  Thus  the  proper 

#  A 

values  will  be  obtained  for  both  Zjr(^j)  and  which  are  the 

impedance  quantities  the  interior  support  must  provide,  respectively, 
looking  to  the  right  and  left  of  the  support  for  a  minimum  vibration 
response  of  the  symmetrical  shaft. 

The  matrix  forms  for  the  impedance  at  each  support  are  as  follows 
(see  Appendix  H  for  explanations): 

#  Z°(U)  ’Z 

z0(°)  = 

-Z0(21)  Z 

#  [Zl(U) 

Z1<*1>= 

0  z 


At  x=0, 


At  x=ij, 
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ZE(11) 

Z2(12) 

s 

Atx=C2,  Z2(£2), 

_Z2(21) 

Z2(22) 

where  ZQ(0)s  I2(£2),  i.e.,  Z0(i j)(0)=Z2(.j)(£2),  i. 

j=l,  2. 

If  the  following  equation, 


t 

'  #  #  #  1 

-  1 

■»  *  #  ■ 

W= 

-C22+Zlr^l^C12 

C21-*lr<VCll 

and  its  related  expressions  (see  Eqs.  (57)  and  (58))  are  used,  then 

rir(*i)n  =  0  leads  to  the  following  expression  which  corresponds  to 

a  quasi- matched  interior  support  (see  the  mathematical  derivations 
in  Appendix  I): 

('  +  e3K^r(22)««1e^ir(Ii))+(e|e2  +  ie2)(zlr(12)zlr(2!)-Zir(u)Zir(22)) 


-  e3(e  2+ie  j  )(z ,  r(i  r(2 1))- (e  j +1^13)- 0 

where 


(62) 


lr<*l>= 


'lr(ll)  Z 1  r(12) 


I  Z 1 r(2  1)  Z 1 r (2  2 ) 


lr(lir2l(ll)+ 


(Hej) 


det 


’#  #  # 

Cll  +  C12r2^1^ 


„  -*  ) 
^(H)(^2)e  -1 
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Zlr(12)  =  2 1  r(  2  1) 


det 


#  #  " 

Cll+Cl2r2^l) 


X 


“(e2+iel^r2(  1  l)^2^e  +(e2~ie  1 ) 


ie  l  e2(  1^e3> 

‘1  r(22)  "  161(22)t  f?  ¥  F 


det  Cll+Cl2r2(fl> 


r2(l  l/f  2^e 


-2ie1N/co(jf  2  j ) 


+  1 


det 


’#  #  # 

Cll+C12r2^l) 


=  “e1 e3 


ie. 


-2ie1  Vu>(jP 2"f  j) 


- 2 i e  (  jP  ^  ~  j) 


r2(n)<f2>  +  1 


r2(U)<f2)  -  1 


r2(ll)<*2>  ‘ 


1 


det 


r  #  u  #i 

"C22  +  22^2^C12 


(1+e3)(z2(22)+ieie2z2(n))  + 


det 


r  #  M  #  i 

"C22  +  22^2^C12 


(e  1  e3+i«2)( z2(  i2)Z2(21)_Z;2(  1 1)Z2(22)^ 

e3(  eZ+ie  1 )(  ZZ(  1 2)  +Z2(2  l))~(el+le2e3) 


vie3+ie1e3z2(11)  +  z2(12))(e2e3-ie2z2^1)+ie3Z2(22))  + 


(e1+ie1e3z2(21)  +  z2(2Z))(l  "ie2z2(11)+ie3z2^12^) 


Tt  appears  from  the  preceding  discussion  and  equations  that  the  example 
problems  have  been  severely  limited  to  several  simple  cases  of  the 
symmetrically  supported  shaft.  That  this  is  not  so,  however,  will  become 
apparent,  if  the  following  principles  are  understood: 

Consider  the  case 


then, 
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Uv^i) : 

this  means  that  besides  establishing  fjr(j£j),  a  second  set  of  conditions  for 
r(£  i )  must  be  met  by  the  interior  support.  Equations  for  evaluating 

^^(f)  may  be  developed  in  accordance  with  the  steps  indicated  in  Appendix 

I.  To  control  the  shaft' s  minimum  vibration  response  completely,  it  is 
necessary  that  the  interior  support  consists  of  two  sets  of  impedances, 

each  optimized  for  its  own  set  of  standing  waves:  ilr(!  j)  for  the  waves 

traveling  to  the  right  and  originating  in  the  [0,  f  J  span,  and  z^  ^(f  j)  for 

the  waves  traveling  to  the  left  and  originating  in  the  [  JP  j sPan-  As  has 

been  suggested  in  the  preceding  chapter,  page  58  ,  optimized  interior 
supports  act  as  one-way  filters  in  that  they  let  incoming  incident  waves 
through,  but  block  the  reflected  waves  approaching  from  the  other  side, 
except  in  the  case  when  the  support  exhibits  impedance  characteristics 
which  satisfy  conditions  for  incoming  waves  from  both  the  left  and  right 
sides.  In  general,  it  is  not  easy  to  build  supports  meeting  such  conditions. 
Therefore,  it  is  suggested  that  the  interior  support  be  placed  as  close  to  an 
end  support  as  possible  so  that  it,  in  effect,  replaces  that  support.  In  the 
example  problems  of  this  chapter  it  would  mean  placing  the  interior  support 
in  close  proximity  of  the  right  end  support  since  Eq.  (62)  is  given  in  terms 

of  the  components  of 


It  follows  from  these  arguments  that  Eq.  (62)  can  be  applied  also  to 
problems  in  which  the  interior  support  is  not  placed  at  midspan,  but  placed 
in  close  proximity  of  the  right-hand  end  support.  In  this  case,  JP  ^  , 


and  inspection  of  Eq.  (62)  shows  that  the  effect 
essentially  eliminated  since  the  term  r2(ll)(fz) 

The  terms  in  Eq.  (62)  simplify  to: 


of  the  right  end  support  is 
must  vanish  for  l  ^  . 


(l+ej) 

Zlr(ll)  =  Zl(ll)  ~ 

-ele3+ie2 


e3(e2-iei) 

Zlr(l2)  =  Zlr(21)  =  ~ 

-ele3+ie2 


ie  1  e2^  ^  +e3) 

Zlr(22)  =  Z1  (22)  +  “ 

"ele3+ie2 
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which  together  with  E  (62)  are  the  equations  which  Nelson  derived  for 
single  span  shafts  (reference  14).  Thus,  the  results  for  the  symmetry 
cases  of  the  shaft  can  be  used  for  or  extended  to  the  problems  in  which 
the  interior  support  e  -'sentially  replaces  one  of  the  end  supports.  It  should 
be  noted  here  that  Eq.  (62)  and  associated  expressions  for  zir(i2)’ 

Zlr(21)'  an<^  zlr(22)  could  be  worked  into  relatively  simple  design  formulas 

in  terms  of  specific  support  conditions,  such  as  translational  supports  or 
rotational  supports,  especially  in  cases  for  which  the  shear  deformation, 
rotational  inertia,  and  gyroscopic  terms  are  neglected;  i.  e.  ,  when 

el  =  e2  =  e3  =  1  • 


EXAMPLES 


To  illustrate  the  existence  of  many  quasi-matched  interior  support 
impedances  yielding  a  "best"  vibration  response  behavior,  and  to  show  the 
manner  in  which  the  impedance  at  each  support  is  calculated,  some  specific 
examples  will  be  studied. 


Three-Support  Shaft  System  with  Specified  Interior  Support  Configuration 


As  a  first  example,  consider  a  shaft  with  both  ends  having  the  same 
configuration,  and  with  its  interior  support  consisting  of  both  a  translational 
and  a  rotational  spring  and  damper  unit  of  negligible  mass  as  indicated  by 

N  .  ni  ,  K  .  n  I  C.rn,  C 


This  shaft  system  is  shown  in  Figure  17.  Consider 
all  forces  acting  on  this  element 

a 

are  also  shown  in  Figure  17. 


IT  ’  “1R  ’  '"IT'  ~1R 
an  infinitesimal  shaft  element  at  x  =  St 


Since,  from  Figure  17,  all  forces  are  assumed  to  be  in  the  positive  sense, 
the  following  force  equilibrium  condition  exists  at  x  =  f  j  as 


F(f  r0)  -  F(f,+0)  =  -R(f  ,) 


and 


ClTYlt*tl)+KlTYl*f  1* 
ClRY2t<*l>+KlRY2»l> 


If  the  above  relationships  are  restated  in  Laplace  transform  form, 


F(fj-O)  -  F(f i+0)  =  -R(£j) 


and 
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R  *  R taction  of  the 
interior  support 


C»RY2t 

C|TY|t  KIT*, 


+  KIRY2 


Figure  17.  Three- support  Shaft  System  With  Specified 
Interior  Support  Configuration. 
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-R(f  !>  =  s 


If  this  expression  is  substituted  back  into  the  first  relationship  for  the 
force  equilibrium  conditions  at  x  =  (  .  and  if  this  expression  is  compared 
with  Eq.  (13),  which  may  be  written  as 

F(«j-0)  -  Fffj+O)  =  sZjlqjYdj)  , 
it  is  concluded  that 


C1  !,+  7  K1T 


C1R+  s  K1R 


Y(Ij) 


# 

Zl(^l)  = 


C1T+  s  K1T 


0 


0 


C1R+  a 


(63) 


which  is  the  interior  support  impedance. 


is  satisfied,  the  values  for  K 


IT 


K 


1R 


If  Eq.  (61)  is  used  and  if  Eq.  (62) 
j  T  ,  and  can  be  obtained; 


the  "best"  vibration  response  possible  for  the  shaft  supported  by  the  type 
of  supports  shown  in  Figure  17  will  result,  provided  that  the  impedance 
values  of  the  end  supports  are  known. 


Three-Support  Shaft  System  with  Floating  Ring  Damper  Assembly  as  the 
Interior  Support 

A  second  and  more  complicated  configuration  of  the  interior  support  is 
shown  in  Figure  18.  In  this  case,  the  impedance  of  the  interior  support  is 
as  follows  (see  the  mathematical  derivations  in  Apoendix  J): 

s4MaMhcts3MaMt,(Kb^Kc)ts[(MatMb)Ka(yKc)+MaKbKc] 

s3MbC+s2Mb(Kb+Kc)+s(Ka+Kc)C+(KaKb+KbKc+KcKa) 

0 


Zl<fl)  = 


+  s[  (Ma+Mb)Ka(Kb+Kc)+MaKbKc]  KaKbKc 

0 

0 
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3 


Ma,  Mb  ar«  mattts 

Ka,  Kb  ,  Kc  art  spring  constant* 


Figure  18.  Three  *  support  Shaft  System  With  Floating  Ring 
Damper  Assembly  at  the  Interior  Support. 
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This  rather  general  configuration  can  be  used  to  represent  a  wide  variety 
of  actual  support  configurations.  In  this  case,  it  was  selected  to  serve  as 
an  idealization  of  the  floating  ring  damper  bearing  assembly  (reference  8). 
Again,  Eos.  (61)  and  (62)  should  be  used  to  establish  the  "best"  values  for 
Ma,  ,  Ka,  Kb>  Kc  and  C. 


Three-Support  Shaft  System  with  Specified  Support  Configurations 


A  third  and  more  specific  shaft  system  is  the  one  shown  in  Figure  19.  For 
simplicity,  all  supports  are  assumed  to  have  the  same  configuration.  The 
support  impedances  which  can  be  derived  in  the  manner  shown  in  Appendix 
I  can  be  expressed  as  follows  (see  Eq.  (63)): 


- 

C0T+s  K0T 

# 

Z0(11)  0 

0 

Z  (o)  = 
o'  ' 

cor+7  kor_ 

_  0  Z0(22)_ 

0 

7 

0 

C,  _+  —  K,  _  0 

# 

1(11) 

IT  s  IT 

w 

0 

Zl(22) 

0  C1R+1K1R 

# 

Z2(ll) 

0 

C2T+7K2T  ° 

Z2<f2>  * 

0 

Z2(22) 

0  C2R+7  K2R 

Since 


>  •  (64) 


1  = 

Z11 

0 

. 

1  2, 

t—  11 

V  w 

0 

0 

Z22  _ 

0 

^  z22. 

1116,1  Zk(  1 1)  =  ^  Zk(ll)  ’  Zk(22)  =  ^  Zk(22) 


(65) 


where  k  =  0,  1,  2  and  s  =  ico  . 


By  substituting  Eqs.  (64)  and  (65)  into  the  quasi-matching  condition  as 
expressed  by  Eq.  (62),  the  "best"  values  for  the  support  parameters  for 
the  shaft  system  shown  in  Figure  19  can  be  obtained.  Thus, 
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L  ie1e2(1+e3)“(eie3+ie2)Q22](ClT"  uK1t)  + 

Vu) 

1+e3)-(e1e3+ie2)Ql  jJ  (Cir‘  uKir)  ■ 

(e!  e3+ie?)  l(C1T’  ”  Kj  T)(C1R-  ~K1r)]  + 

L  ( 1  ^ e3)( Q2 2  +  ie  1  e2^1 1 ) ~ 2 e3( eZ+ie  1  )^12+(el  e3+ie2^l2"^l  1^22^ 

(e  i  +ie2e^)  ]  =0  (66) 

where 


Q 


1 1 


Q-, 


Q 


12 


det 


(1^) 

det 

FT  2  2  T 

Cll+C.2r2<fl) 

*> 

ieje2(  i+ep 

det 

cntci2r2<ci> 

e3 

T  1  ¥  T 

det 

cn+ci2r2<V 

L 

#  A 

f  #  1 

r2(ll)^2^e 


-2i e^^f  2-f  j) 


-1 


-2ie  .) 

r2(ll)^2)e  +1 


-2ie  .*^5 1$  , ) 

■(e2+iel>r2(ll)<£2>e  -Ke2-iel) 


Cll+C12r2<£l> 


"ele3 


le 


- 2ie  i ) 

e 


r2(ll)^2)  +  1 


# 

*2(1 1)^2^ 


det 


■  |  # 

■C22+^2*£2*C12 


r2(ll)(£2)  -  1 


<1+C3>  ^C2R"K2R)+ 


lel'!2 

\fZ) 


(C2T~  oj  K2T^ 


(e1e3+ie2)(C2T“  u)K2t^C2R"  u  K2R^‘ 


(e1+ie2e3) 
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det 


~CZZ^Z^2)C\l 


Z 

1C3 


1+-tt(c2t“z:kzt) 


,+i^c2R-  w  K2r) 


+ 


le. 


1 "  —  (C2T"  co  K2t) 


n/co 


el^C2R-  w  K2r) 


As  can  be  seen  from  the  above  expressions,  the  only  unknowns  involved 
are  K._,  K1T5,  C._,  C.  _  of  the  interior  support,  provided  that  both  end 

support  configurations  are  given.  Hence,  Eq.  (66)  provides  a  means  of 
establishing  the  values  that  should  be  used  for  the  unknowns  ,  K^R, 

Cj,p  ,  and  CjR  in  order  to  satisfy  the  quasi-matching  condition  for  the 

interior  support  in  the  shaft  system  shown  in  Figure  19  The  interesting 
point  here  is  that  the  set  of  values  for  K^,  K^R  ,  C^,  C^R  required  to 

satisfy  Eq.  (66)  is  not  unique.  This  permits  a  wide  range  of  selection  of 
these  parameters  to  match  the  feasibility  and  availability  of  various  damper 
and  spring  materials  and  configurations. 

For  a  shaft  system  similar  to  the  one  shown  in  Figure  19,  but  which  has  no 
rotational  spring  and  damper  units  attached  to  its  supports,  the  quasi¬ 
matching  condition  of  the  interior  support  can  be  written  simply  by  setting 
=  Crr  =  0  and  k  =  0,  1,  2,  in  Eq.  (66): 

cit-^kit  +  q  =  0  <67> 


where 


Q  = 


iele2(1+e3)“(ele3+ie2>Q22 


( 1  +  e3)(Q22+iei  e2^i  i)-2e-j(e2+ie1)Q12  + 


(eie3+ie2)(Q12-Q^1Q22)-(e1+ie2e3) 


^1 1  ’  Q2 2 ’  Q12 ’  det 

Eq.  (66), 


have  the  same  expressions  as  in 


rz(  =  fa  :  1  =1 

+lz«z>ci; 


det 


22 


ie  ie2(  l+e3) 

4Z 


(C2T~  w  K2T^"^el+ie2e3) 


det 


C22+^2^2)C12 


e2^  ^ 


(C2T~  w  K2T^+^el"ie2e3^  ' 
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The  quantity  Q  is  a  given  constant  which  relates  to  end  configurations, 
shaft  characteristics,  and  rotating  speed.  Eq.  (67)  demonstrates  the 
usefulness  of  the  transmi ssion  line  analogy  by  quasi-matching;  C^p  and 

Kjp  are  the  only  unknowns  appearing  in  a  linear  equation  and  can  be 

determined  directly. 

In  the  cases  for  which  the  interior  support  is  placed  in  close  proximity  of 
the  right  end  support,  Eqs.  (66)  and  (67)  can  be  used  also  for  establishing 
the  optimum  impedance  terms  of  the  interior  support  by  letting  j  ^(f  2) 

Simplification  of  Quasi-Matching  Conditions  by  Neglecting  Rotary  Inertia, 
Shear,  and  Gyroscopic  Effects 

It  can  be  verified  that  the  differential  equation  accounting  for  the  bending 
effects  in  the  shaft  only  can  be  obtained  from  Eq.  (7)  simply  by  setting 
e^  =  e^  =  e^  =  1  .  In  practice,  if  the  rotating  speed  is  not  extremely  high, 

the  terms  for  rotational  inertia,  shear  deformation,  and  gyroscopic  effects 
can  be  neglected.  When  the  above  argument  is  used,  the  qua  si- matching 
equations  for  a  shaft  with  one  interior  support  may  be  simplified  accord¬ 
ingly.  With  reference  to  Eq.  (62)  and  Figure  1 6 , 


2^zlr(22)+izlr(ll)^+^1+i^zlr(12)zlr(21)‘zlr(ll)zlr(22)"zlr(l2)“zlr(21)“1^  0 

(68) 


where 


zlr(ll)  =  Zl(ll)+ 


lr(l2)  =  Z 1  r( 2 1 )  "  A 


r2(ll)^2^e 


-2i^f2-f1) 


l  L 


'(1+1)r2(l  1)^2^ 


-21^2-Ij) 


+(l-i) 


;lr(22)  =  z  1  ( 2  2 )  + 


2i 


r2(ll)^2^e 


-2i\ Qtz-flJ 


+  1 


A1  =- 


-2i j) 


r2(ll)^2)+1 


-  i 


■2i  N^(f2-£1) 


r2(ll)^2^  "  1 


r2(ll)^2^  /  2^Z2(22)+1Z2(11)^+^1+1^Z2(12)Z2(21)'Z2(11)Z2(22)' 

^2  L 


z2(  12)-z2(21)~  ^ 


83 


^2 


1  +  1Z2(11)+Z2(12) 


1  +  iz2(21)+Z2(22) 


1  ‘  iz2(21)+iz2(22) 


1  “  1Z2(11)  +  1Z2{1  2) 


For  the  specific  shaft  configuration  of  Figure  19,  the  corresponding  quasi¬ 
matching  condition  can  be  rewritten  in  the  following  simplified  form  (see 
Eq.  (66)): 


1 


2i-(l+i)Q,2 


<cit-STkit> 


2-(l+i)Q 


1 1 


^C1R"  co  K1R^ 


(1+1)(C1T“  w  KiT)(C1R-  co  K1R}  + 


2(Q22+iQll>  + 


(1+i)(Qj2-QIiQ22‘2Qi2-1)  =0 


(69) 


where 


'll 


1L 


-2WU  (K-J?  ) 

r2(n)W  2>e  -  1 


2i 


'22 


r2(ll)^2^e 


-2i4ZUz-6  {) 


+  1 


Q 


_  1 


12 


-<1+i)r2(ll)<*2>e 


- 2i\rco  (iz-(  j) 


*1 


■2i\ fZ  {iz-ix) 


r2(ll)(f2)  +  1 


-1 


+  (1-i) 


r2(ll)<C2) 


r2(ll)^£2^ 


A 


(1+i) 


^  ^C2R"  CO  K2R^  +  r-  ^C2T_coK2T^ 

V  co 


(C2T"  u  K2T^C2R‘  co  K2R^  +  1 


*2 


-l 


1  + 


-  IC2T‘  u  K2T* 


1 


For  the  simplified  shaft  configuration  in  which  no  rotational  springs  and 
dampers  are  being  used,  a  very  simple  qua  si- matching  equation  for 
determining  optimum  damping  and  spring  coefficients  in  the  interior 
support,  can  be  obtained.  With  reference  to  Eq.  (67), 

cit-^kit  +  q=0  <70> 

where 

2(Q22MQu)+(Hi)(Q^-QnQ,2-EQl2-l) 


Q 1 1  ,  Q->2  *  Q j ^  have  the  same  expressions  as  in  Eq.  (69) , 


(1+i) 


r2(ll)^2^ 


2i 
L  n/u 


r  (C2T‘  cj  K2T^ 


A,  =  -  -  LK?J)  +  (1-i) 


-IS 


2T  u  2T; 


Thus,  fairly  straightforward  solutions  have  been  obtained  which  permit  the 
designer  to  obtain  in  a  direct  manner  the  spring  and  damping  coefficients 
required  in  the  interior  support  for  the  minimum  vibration  response 
operation  of  a  hypercritical  shaft  on  three  flexible  supports.  To  suit  cases 
other  than  those  discussed  specifically  in  this  report,  Eq.  (68)  can  be 
manipulated.  Again,  it  should  be  noted  that  for  the  cases  in  which  the 
interior  support  is  placed  close  to  the  right  end  support,  Eqs.  (68),  (69), 

1  1)^2^  =  0  •  It  should  be 

pointed  out  also  that  the  effect  of  mass  in  the  interior  support  can  be 
included  directly  in  Eqs.  (66),  (67),  (69),  and  (70)  by  replacing  the  spring 
coefficients : 

K1T  =  (kit  "  “X)  ’ 

k2t  =  (k2t  “ w  ivy  ’ 

K1R  =  (kir  "  w  V  ’ 

K2R  =<K2R  -  ’ 

where  Mj  and  M2  are,  respectively,  the  mass  of  the  support  systems  at 
x  =  H  j  and  x  =  and  1^  and  are,  respectively,  the  mass  moments 
of  the  support  systems  about  their  center  of  gravity  at  x  =  f  ^  and  x  =  fL^ 
of  the  shaft. 


and  (70)  can  be  simplified  further  by  letting  r2( 
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Design  Formula  for  the  Shaft  with  Its  Interior  Support  Placed  in  Close 
Proximity  of  One  of  Its  End  Supports 


For  the  case  in  which  the  shaft  is  being  treated  essentially  as  a  two- support 
system  (i.  e.  ,  the  interior  support  is  placed  in  close  proximity  of  one  of  the 
end  supports),  it  is  possible  to  develop  a  very  convenient  formula  for  the 
determination  of  the  optimum  quantities  (in  terms  of  quasi-matching)  which 
the  interior  support  must  provide  for  a  minimum  vibration  response  of  the 
shaft.  In  this  case,  it  is  assumed  that  ej  =  =  1  .  If  the  impedance 

of  the  interior  support  is  considered  to  act  in  a  translational  fashion  only, 
the  impedance  matrix  of  the  interior  support  can  be  written  as  follows: 


# 

Zl(*l> 


IT 


"  co  (K1T_ 


co2Mj) 


0 


0 


Substitution  of  this  condition  in  Eq.  (70),  and  letting  r2(  1 1  )(f  2)  =  0  »  results, 
after  some  algebraic  work,  in  the  following  formula: 


c._  =  — <  U3  + 

1  1  2u) 


3/2  x  ?IXA  2 

CO  T  Z(MjC0 


K  ) 

IT* 


2  >  1/2 


(71) 


where 


IT 


.  ClTCs 
E  A 

y 


co 


CO  ' 


*b 


M, 


K 


M 


IT 


PARb 


K*ITRb 
E  A 

y 


c 

9  V  o 


*  *  * 

Also,  C  ^  ,  M  ^  ,  and  K  ^  are,  respectively,  the  dimensional  quantities 

for  the  damping  coefficient,  mass,  and  spring  coefficient  of  the  interior 
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support,  and  «*  is  the  di  mensional  quantity  for  the  angular  velocity  of  the 
shaft. 

Eq.  (71)  will  calculate  the  optimum  damping  coefficients  for  preselected 
sets  of  masses,  and  spring  coefficients  in  the  support  in  terms  of  the 
angular  velocity  of  the  shaft.  It  rhould  be  note'!  that  the  values  for 

are  independent  of  the  impedance  conditions  of  the  end  supports  and  that 
the  values  for  optimum  damping  for  two  interior  supports  placed  close  to 
each  end  support  are  the  same  as  that  for  one  interior  support  placed  close 
to  one  or  the  other  end  support. 
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support,  and  a>*  is  the  dimensional  quantity  for  the  angular  velocity  of  the 
shaft. 

Eq.  (71)  will  calculate  the  optimum  damping  coefficients  for  preselected 
sets  of  masses,  and  spring  coefficients  in  the  support  in  terms  of  the 
angular  velocity  of  the  shaft.  It  should  be  noted  that  the  values  for 

are  independent  of  the  impedance  conditions  of  the  end  supports  and  that 
the  values  for  optimum  damping  for  two  interior  supports  placed  close  to 
each  end  support  are  the  same  as  that  for  one  interior  support  placed  close 
to  one  or  the  other  end  support. 
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CHAPTER  5 


CONCLUSIONS 

EXTENSION  OF  TRANSMISSION  LINE  ANALOGY 

The  transmission  line  analogy  can  be  extended  to  shafts  having  any  number 
of  interior  supports.  However,  the  solution  is  considerably  more  compli¬ 
cated  for  shafts  having  more  than  one  interior  support. 

MATCHED  END  IMPEDANCES 


For  the  shaft  with  both  end  impedances  matched  with  the  characteristic 
impedance  of  the  shaft,  no  intermediate  support  is  needed  to  assist  in  the 
minimization  of  vibration  response. 

QUASI- MATCHING 


If  it  is  physically  impossible  or  impractical  to  terminate  a  shaft  in  its 
characteristic  impedance,  quasi-matched  end  impedances  or  quasi-matched 
interior  supports  (when  end  conditions  are  not  available  for  optimization) 
should  provide  good  performance.  Quasi-matching  involves  the  selection 
of  support  conditions  in  such  a  way  that  the  predominant  Tj  ^  term  of  the 

reflection  matrix  vanishes.  The  concept  of  qua  si -matching  is  based  on  the 
assumption  that  expf-e^v^T  «  1  ,  which  is  true  as  w  increases. 


ONE  OPr  I  MUM  INTERIOR  SUPPORT 

When  the  shaft  and  end  support  impedances  are  not  matched  and  when  only 
one  interior  support  is  used,  two  approaches  may  be  employed  to  assist  in 
the  minimization  of  vibration  response: 

1.  If  both  ends  have  the  same  configuration,  the  interior  support  may 
be  placed  at  the  midspan,  provided  the  conditions  of  this  interior 
support  obey 

rir(ll)^l)  5  ri£(ll)^l^  =  ° 

as  indicated  by  Eq.  (62). 

2.  If  one  end  support  is  different  from  the  other,  the  closer  the 
matched  intermediate  support  is  placed  to  one  of  the  ends,  the 
more  effectively  it  will  minimize  the  vibration  response  of  the 
shaft  with  provision  that: 


a. 


If  the  interior  support  is  placed  close  to  the  left  end,  it  should 
be  arranged  such  that  £(ll)^l)  =  ®  • 


b.  If  the  interior  support  is  placed  close 


li  IUC  UHCUUi  O  uppui  L  X  O  Uit  -C  UL  tU  UiC  XXgAH  Air 

should  be  arranged  such  that  r(  1 1 )(^  1 )  =  ^  ecluat^on 


89 


for  Fj 1 1 )( "^  1 )  =  ®  ^ave  the  same  form  as  Eq.  (62), 

except  that  f  ^  (see  Figure  4)  indicates  the  actual  location 
of  the  interior  support. 


TWO  OPTIMUM  INTERIOR  SUPPORTS 


When  the  shaft  and  end  support  impedances  are  not  matched,  the  U3e  of  two 
matched  intermediate  supports  placed  close  to  the  ends  of  the  shaft  is 
recommended.  This  is  equivalent  to  letting  n  =  3  for  the  shaft  system 
shown  by  Figure  4,  in  which  St .  ,  H  7  should  approach  0  ,  J P  -  ,  respectively, 

^  t  ^  |  *  w  J 

and  the  1  and  the  2  interior  supports  should  be  designed  in  terms  of 
configurations  which  satisfy  the  conditions 

rue(n)^i) =  0 


and 


r2r(  11)^2)  0  * 

In  other  words,  only  the  incident  wave  exists  on  [fj  ,  S  ,  which  is 
almost  equal  in  magnitude  to  [  0  ,  4^]  »  namely,  the  whole  portion  of  the 
shaft. 
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APPENDIX  A 


MATHEMATICAL  DERIVATIONS  OF  THE  ROTATORY 
AND  GYROSCOPIC  EFFECTS 


In  Figures  20,  21  and  22,  a  right-hand  screw  convention  is  used  to  indicate 
the  positive  directions  of  the  applying  torques  on  the  differential  mass 
element  of  a  rotating  shaft.  If  it  is  assumed  that  the  transverse  deflections 
(Y’l'j  and  Y*^)  of  the  shaft  during  vibration  are  very  small  and  that  the 

center  of  gravity,  o,  of  the  differential  mass  element  coincides  with  the 
axis  of  the  shaft,  the  position  of  this  element  will  be  completely  determined 
by  the  coordinates  Y|  ^  and  Y^  of  its  center  of  gravity,  o,  and  by  the 

angular  rotations,  Y|  ^  and 


The  conditions  assumed  here  correspond  to  the  case  of  a  vertical  shaft 
when  the  weight  of  the  differential  mass  element  does  not  affect  the 
deflections  of  the  shaft.  Under  these  conditions,  if  W  equals  the  weight 
of  the  differential  mass  element  and  if  only  the  elastic  reaction  of  the  shaft 
is  taken  into  consideration,  the  equations  of  motion  of  the  center  of  gravity 
of  the  differential  mass  element  are  as  follows  (see  Figures  21  and  22): 

In  -X^  direction, 

w 

Inertia  force  =  — 
g 

In  -X^  direction, 


Inertia  force  =  —  Y|^t 
g 


w 

where  —  =  pAdx  , 
g 


The  equations  of  relative  motion  of  the  differential  mass  element  with 
respect  to  its  center  of  gravity,  o,  will  now  be  obtained  by  using  the 
principle  of  angular  momentum,  which  states  that  the  rate  of  increase  of 
the  total  moment  of  momentum  of  any  moving  system  about  any  fixed  axis 
is  equal  to  the  total  moment  of  the  external  forces  about  this  axis.  In 
calculating  the  rate  of  change  of  the  angular  momentum  about  a  fixed  axis 
drawn  through  the  instantaneous  position  of  the  center  of  gravity,  o,  the 
relative  motion  alone  can  be  taken  into  consideration. 

In  calculating  the  components  of  the  angular  momentum,  the  principal  axes 
of  inertia  of  the  differential  mass  element  will  be  taken.  The  axis  of 
rotation,  oo,  is  one  of  these  axes.  The  two  other  axes,  oa  and  ob,  will 
be  any  two  perpendicular  diameters  of  the  element. 
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a  Differential  Mass  Element. 
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m 


m 


=  mass  moment  of  inertia  of  the  mass  element  about  the 
axis  oo, 

=  mass  moment  of  inertia  of  the  mass  element  about  a 
diameter. 


Observe  that 


m 


=  21 


m 


I'  u  =  component  of  angular  momentum  about  the  axis  oo, 

”^mYIlt  -  component  of  angular  momentum  about  the  axis  oa, 

I  =  component  of  angular  momentum  about  the  axis  ob. 

Since  Y^j  and  are  assumed  to  be  small,  then  Yjjj.^  and  Y?^  will  be 

approximate  values  of  the  angular  velocities  about  the  diameters  oa  and 
ob.  If  components  of  the  angular  momentum  on  the  fixed  axes  oX^  and 

oXj  are  projected  through  the  instantaneous  position  of  the  center  of 

gravity,  o,  the  resultants  of  the  angular  momentum  in  the  X^-  and  X-,- 

directions  can  be  obtained.  It  is  shown  that  cosYjj^  and  cos' 

that  sinY^  Yh  and  sinYii  *  Yh  • 


iY^i *  1  and 


In  the  X^- direction,  with  reference  to  X^X^-plane  projection  (see  Figure 

21), 


In  the  X-, -direction,  with  reference  to  X^X^ -plane  projection  (see  Figure 

22), 

(ImY?2,X‘:osY?l)-(I^o)<sinY!l)  "  1taY!2t-Im"oY!l  ■ 

Then,  from  the  principle  of  angular  momentum  and  since  I'  =11  the 

mm’ 

following  can  be  shown: 


In  X^-direction, 


F*.  =  I  Y*,  +21  u  Y* 

21  m  2 1  tt  m  o  22t 


In  X-,  -  direction, 

F?2=f 

at 

or 

F?2=ImYhtt-2Im“oY?lt  ' 

If  I  is  the  moment  of  inertia  of  the  cross  section  of  the  differential  mass 
element, 

I  =  pldx 
m  F 


or 


I'  =  21  =  2  pldx  . 

mm 


Then,  the  above  expressions  for  moment  of  momentum  may  be  rewritten 
as : 


In  X^-direction, 

Fh  *  piYSmdx+2I“o‘’YS2td* 


In  X^- direction, 


F$2  =  pIYl2ttd31'2I“opY!ltdx 
which  were  shown  in  Figures  2  and  3  in  Chapter  1. 
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F?1  =ImY!l..+2Im“oY?2t  ' 


In  X-,- direction, 

r?2  ^fiysu-wsi) 

at 

or 

F!2=ImY?2,t-2Im“oYil«  • 

If  I  is  the  moment  of  inertia  of  the  cross  section  of  the  differential  mass 
element, 

I  =  pldx 
m  r 


or 

I'  =  21  =  2  pldx  . 

m  m 

Then,  the  above  expressions  for  moment  of  momentum  may  be  rewritten 
as  : 


In  X^- direction, 


Fh  =  pIYhttdx+2L>oPYhtdx 

In  X^-direction, 


F?2  =  <>IYS2ttdx-2I“o'>YSl«dx 

which  were  shown  in  Figures  2  and  3  in  Chapter  1. 
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APPENDIX  B 


IMPORTANT  FUNCTIONS 


^e'-K 
wo  ^"T”) 


e2  = 


V  I+e'a)(2w  -u>) 


where  w  =  -is  and  u  =  u  for  steady- state  condition. 

o 
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-ie.e-  -ie_ 

13  2 
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0  -1 
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‘l  0 

f  3/2 

0  ' 

# 

# 

# 

CJ 

c_ 

-0  -1 

o 

Hi 

II 

.0 

U>  . 

where  p  is  an  arbitrary  real  number  which,  once  chosen,  is  the  same  for 
each  p. 


.  _i_ 

1 

0 

# 

z 

1 

0 

• "  vs 

_0 

sTu_ 

s 

_0 

n/gj  _ 
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Z  =  characteristic  impedance  of  the  rotating  shaft 

0 


#  *  #-i 

*.  =c2icn 


t  =(• 

a  ' 


l+e. 


2,24 

e2+elc3 


(*le3+ie2> 


2  2  2 
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/  ,.e2“cle3\ 

(ele2+1 - —  ) 
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APPENDIX  C 


MATHEMATICAL  DERIVATIONS  OF  APPLYING  BOUNDARY 
CONDITIONS  AT  SUPPORTS 

GENERAL  SOLUTION 


The  general  solution  of  the  equations  of  motion  is  shown  in  Eq.  (10)  as 


and  is  written  out  as 


#  r#  #  #  # 


Y(x)  =7  c  CwR;x)q+C12R(-x)r 


,  #  r#  #  #  # 

x)  =  C£  C21R(x)q+C22R(-x)7 


I1.  ;:ie  general  solution  is  applied  to  each  span,  the  following  set  of  _ 
equations  can  be  obtained  with  integration  constants  q^  ,  q2  ,  ....  q^_  ^ 

*  qjrJ  qj+l . Vl’  V  rl  ’  r2  *  *’•  ’  rj-l  ’  rjC  ’  rjr  ’  rj+l . 

n_  i  »  rn  (aee  Figure  4): 

For  [0,  ij]  . 

~  #r##  #  # 

Y(x)  =iCy|cnR(x)i1+Ci2R(-x)7iJ  , 

»  #[#»_## 

F(x)  =  Cf  C21R(x)q1+C22R(-x)r1 


Uj-2*  fj-l] 


#  r#  # 


#  # 


-<*>  =tCy  CnR(x)qj_i+cl2R(-x)rj.i 


*  #  r#  #  #  # 

F(x)  =  Cf  LC21R(x)q^_1+C22R(-x)r^_  1 


For  ^j-1’  ajk] 
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Y(xl  =  -  C 
'  s  y 


#  # 


#  # 


C1  lR(X)qjl  +C12R<"X)rjf 


F(x)  =  Cf 


'  #  #  #  # 
C2iR(x)cljf  +C22R^’x^rjj? 


For  [  ajk  ,  H.]  , 

-  !  # 
Y(x)  =  i  C 


#  # 


#  # 


C,  ,R(x)q^+C,7R(“x)r; 


s  y  11  '  /njr  12 


F(x)  =  Cf 


'#  #  #  # 

C21R(x>V+C22R<"x>7 


jr 


For  [*.  ,  £  j  +  1]  , 


«  ,  #  #  #  _  #  #  _  ’ 
Y(x)  Cy[cnR(x)q.^  ;i2R(*x)rj+1 

=*  #  r#  #  _  #  #  _ 

F(x)  ~  cf  C2lR^x)qj  +  l+C22R^“X^rj+l 


For  [ f  ,  ,  JP  ] 
L  n- 1  nJ 


-  1  * 
Y(x)  =  -  C 
'  '  s  y 


r#  # 


#  # 


CllR(x)qn+C12R("x)rn 


F(x)  =  Cf 


'»#_## 

C2iR(x)qn+ci2R(-x)rn 


APPLYING  THE  BOUNDARY  CONDITION  AT  x  =  0-0 

~  #  =* 

-F(0)  =  s  ZQ(0)  Y(0)  . 

If  x  =  0  is  substituted  in  the  general  solution  for  [  C  , 
then 


#  # 

and  if  R(0)  =  I  , 
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■*  1 

Y(°)=7Cy 


CllVC12rl 


F(0)  =  Cf 


C2lVC22ri 


and  then  substituting  back  into  the  boundary  condition, 


# 

•C, 


'#_#_■ 

C21ql+C22rl 


1 


#  r# 


=  sZft(0)  -  C 
O'  '  s  y 


Cllql+C12rl 


or 


'#  # 

C21ql+C22rl 


#  .  # 

=  VZ0(0)C 


#  r# 

y 


Cllql+Cl2rl 


By  direct  matrix  algebraic  operations, 


u  l  #  #  , ## 

4  =  — —  FZF  =  Cj  ZC 

^  1  ' 


F  = 


0 


Hence, 


‘#  _  # 

C2lql+C22rl 


■  4o<°> 


Cllql+C12rl 


or 


#  u  # 

C21+V°)C 


11 


V 


czz40(o)c 


12 


rl  =  0 


By  definition  (see  Appendix  B), 


-1 

0 


C22+i0<°> 

*n«>) 


"l 

o" 

*  1 

0 

-1 

C,2 

"l 

o" 

#  1 

0 

-1 

cn  ' 

/ 

By  premultiplying 


1  0 
0  -1 


C21  + 
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If  the  "reflection  matrix"  (looking  to  the  left)  at  x  =  0  is  defined  as 


then  the  preceding  equation  may  be  rewritten  as 

q,  =  r0(o)7[ 

which  forms  the  basis  for  Eqs.  (72).  If  this  expression  is  substituted  back 
into  the  original  general  solution  for  [0  ,  £ j]  , 

#r###  ## 

Y(x)  =  i  Cy  |c,  -x)7, 

or 

*  ,  #  r#  #  #  #  #i# 

Y(x)  =iCy  CnR(x)r0(0)R(x)+C12  R(-x)r1  . 


By  defining 

#  #  IV  # 

r0(x)  =  R(x)ro(0)R(x)  f o r  [  0  ,  (  j  ]  , 

then 

YW=iCy[cnr0M+cu _R{-*)7,  . 
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Similarly, 


F(x)  =  C. 


'U  U  U 

CZU^C12 


R(  ”x)  r  1 


APPLYING  BOUNDARY  CONDITIONS  AT  x=f.,  j  =  1,  2,  ....  n-1 

-  -  -  -  —  J-  -  -  -  - 


If  the  general  solution  is  substituted  directly  into 


Y(fj-O)  =  Y(« .+0)  , 

then 


1  # 
-  C 
s 


V 


‘  #  #  jj  # 

cnR':fj>Vci2R!'V?j 


i  » 

=  -  C 
s  y 


'  #  #  #  # 
CUR^j)qj+l+C12R(‘£  j)rj+l. 


or 


##  #  #  ##  ## 

C..R{f.)q.+C.7R(- It. )r.-C.  .  R{l  .)q. .  , -C.  ,R(-£.)7. .  . 
ll  "  j/nj  12  '  y  j  li  •  j/xj+i  12  s  j7  j+l 


=  0 


j  =  1,  2,  .  .  .  ,  n-1 

which  forms  the  basis  for  Eqs.  (73).  If  the  general  solution  is  substituted 
directly  into 


F(L-O)  -  F(C.+0)  =  sZ^Yd.) 


then 


■#  #  t  #  # 

# 

"#  #  #  *  # 

cfc2.-zj(h)CyCn. 

R(£.)q.+ 
J;  J 

CfC??-Z  .(£  .)C  c,7 

L  f  22  j'  j'  y  12  J 

#  #  #  #  #  # 

CfC21R({j)qj  +  rCfC22R(^j)rj+l  =  0 


#_1  #.j##  m 

If  C^.  is  premultiplied  and  if  Z  C^.  -  &  is  used  again,  then 

# 

R{-£  )r 

x  y  J 

a  h  #  # 

C21R,'tj)Vl'C22R(‘lj)7jU  =  ° 
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# 

1C21 


Rff  .)q.- 

'  F  J 


# 

-c 


22+l/*j)ClZ 


j  —  If  2,  •  .  •  i  1  • 


This  equation  forms  the  basis  for  Eqs.  (74). 

APPLYING  BOUNDARY  CONDITIONS  AT  x  =  f  +0 

n 


F(jf  )  =  sZ  (jf  )  Y(jf  )  . 

'  n7  n'  n7  '  n7 

If  x  =  Jf  is  substituted  into  the  general  solution  for  [  jf  .  ,  jf  1,  then 
n  6  1  n- 1  nJ 


1 


# 


Y(jf  )  =  -C 
'  n7  s  y 


'#  # 


C,,R(£  )q  +C.?R(-f  )r 
1 1  '  n7  m  12'  n7  n 


#  r#  # 


F(«„)  •  Cf 


#  # 


C21R<WC22RHn>rn 


If  these  expressions  are  substituted  back  into  the  boundary  condition 
equation,  then 


#  #  #  #  „ 

C21R'-‘n)VC22R(-{„)rn  =  VCn> 


'#  #  #  # 
CllR<fn)VC12R(-f„)?n 


which  may  be  rewritten  by  solving  for  r^  as: 


r  =  R{jf  ) 
n  '  n7 


-C 


,74  (i  )c, 

22  n'  n7  1< 


-1 


#  1  # 


'21 


R<*„>% 


If  the  reflection  matrix  (looking  to  the  right)  at  x  =  f  n  is  defined  as 


# 

r  in  #  1 

_  1 

r#  „  #  i 

r  (f  )  = 

n'  n7 

-C,,+t  (Jf  )C17 
22  n'  n7  12 

•  1 

C,,-l  (jf  )C., 
21  n'  n7  11 

then 


r  =  R(jf  )T  (jf  )R(jf  )  , 

n  v  n7  nv  n7  v  n7 


which  may  be  rewritten  as 


■  rJ°K 


by  defining 

#  #  #  # 

r  (x)  =  R(jf  -x)r  (Jf  )R(jf  -x) 
nv  7  '  n  7  n'  n7  v  n  7 
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# 

The  next  to  the  last  equation  forms  the  basis  for  Eqs.  (7  5).  If  =  T^Ojq^ 
is  substituted  back  into  the  original  general  solution  for  [^n_|  ,  ,  then 

Y«  =  TCy[Cll+C12rn(x>R;x)qn  • 

Similarly, 

-  #  r#  #  #  1  # 

F(*>  *  Cf  C21+C22rn(x)J  R(x)%  • 

Before  going  any  further,  all  the  equations  obtained  from  boundary 
conditions  of  all  supports  are  summarized  as  follows: 

# 

q,  =  r0(0)r,  (72-1) 

##  #  y  ##  ## 

CllR(Jl)5’l+C12R('l!l)7rCllR(!tl)VC12R<'fl)T2  =  0  (?3-l) 

^  ft  ^  ^  ^  ^  ft  ft 

CllR^2^q2+C12R'’^2^r2“Cl  lR'x^2^cl3~C12R^~^2^r3  =0  (73"2) 


##  ##  ##  ## 

CUR'tj-l^j-l+C12Rf‘tj-l)rj-rCllR^tj-l^jt'C12R^*fj-Prjl  *  0 

(73-(j-l)) 

CiiR(Jj)qj/CuR(-t.)7.r-CnR({j)qj+1-^2k(-f.)7jtl  =  0  (7J-J) 


V!,n.1IVltSR(-|,,)t„.,-CuR"».l)’n4^<n.l)'n  "  0 

(73-(n-l)) 
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#  u  # 

C214l(fl)C 


1 1 


‘  #  #  # 
-c224iUi)c 


12 


RH^r, 


#  #  #  # 

C21R^l^q2“C22R^“^l^r2  =  0 


(74-1) 


*  M  f 
C2l'  2^2^C1 1 


R(c2)q2 


■C22  +  Z2^2)C12 


RH2)r2‘ 


#  #  #  # 

C21R^2^q3"C22R^"^2^r3  =  0 


(74-2) 


#  1  # 


C2l4j-l(*j-l)CU  R(fj-l)qj-l*  'C22+^j-l^j-l)C12 


§  #  #  # 
C2lR^j-l)qj£"C22R("£j-Prj£  =  0 


R^Vi)rj-r 

(74-(j-l)) 


r#  u  #  • 

C21*VVC11 


R(f  .)q 

V  /V 


‘  #  „  #  • 
“C22+Vfj)C12 


R<“Vrjr' 


^  ^  fj 

C21R(f  j)qj+rC22R("f  j)rj+l  =  ° 


(74-j) 


r#  «  #  i 

C,.-f  ,(£  ,  )C .  . 

21  n-lv  n-1'  11 


*'Vi>vr 


#  #  #  # 

C21R(fn-l>%-C22R(-£n-1>7n  *  0  <74-(n-l)) 


(75-n) 


where 
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# 

r0<°>  = 


r  (f  )  = 

n'  n7 


-C22+V°)C12 


-1 


#  A  # 

C.,-^n(0)C 


'21 


11 


'  #  „  #  ‘ 

-C22+WC12 


»  #  ' 

C„-l  (f  )C„ 

21  n'  n7  11 


#  #  #  # 

r0(x)  =  R(x)r0(o)R(x) 

#  #  #  # 

r  (x)  =  r(jC  -x)r  (£  )r(£  -x) 
n'  7  '  n  7  n'  n7  '  n  7 


Now,  Eqs.  (72-1),  (73-1), and  (74- 1)  will  be  considered  first.  If  (72-1)  is 
substituted  into  (73-1), 


#  1 

r#  #  #  i 

1 

r#  #  #  #  1 

jo 

1 

►— 

►— * 

II 

Clir0^Cl^+C12 

L  j 

-  I 

C1  lR^  l)q2+C12R^"^  l^r2 

If  (72-1)  is  substituted  into  (74-1), 


#  #  # 

R{« ,)r0(o)R(f ,)- 


# 


#  #  #  # 

=  C2,R«i)q2+c22R(-£i)72  . 


Again,  if  the  underlined  terms  are  substituted. 


C22+C2ir0^  1^ 

-  - 


r#  #  #  i 

C12+CllV£l) 


X 


#  #  # 

C1 1R^  l^q2+C12R^"^  l^r2 


#  #  #  # 

=  C2iRUi)q2+c22R{^i)r2 


9 


or 


-1  0 

/ 

“#  #  # 

1  0 

1  0 

■#  #  #  1  ] 

_  o  1. 

< 

s 

C21+C22r0^1^ 

+ 

_0  -1 

0  -1 

— 1 

_C11+C12r0^l)_  j 

"#  #  # 

_  1 

1  0 

-1 

r#  #  #  # 

Cll+C12r0^P 

“  1 

_0  -1 

CUR(£l)q2+c12R(-f,)r2 

-1 

0~ 

#  # 

-1 

o" 

#  # 

0 

1 

C22R<b>V 

0 

1 

C21R^"^l^r2  ’ 
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so  that 


r#  #  #  i  r  #  #  # 


C21+C22r0<*l> 


cn+ci2ro<{i> 


'1+Vi> 


r#  *  »  * 

C12R«l>q2+CllR<-*l>r2 


#  #  '  #  # 

=  C22R(fl)q2  +  C21R(-£l)72 


It  will  be  shown  in  Appendix  D  that 


A  ,  v 

Vx> 


'#  #  # 

C21+C22roW 

.Cll+C12r0<*) 

or 


V*i>  * 


Hence,  by  substituting. 


Vi>+Vi> 


r#  # 


Cll+C12r0^1^ 


#  # 


-1 


#  #  #  # 


C12R<Cl>VCllR<-£l>r2 

By  observing  that  the  total  impedance  at  x  =  £^+0,  looking  to  the  left,  is 

*u<£i>  =  *i<£i>  + 
then 


*u«i> 


'#  #  #  # 

C12R(£1>VC11R(‘C1>72 


#  #  #  # 

=  C22R<f1^2+C21R('£l)T2 


Solving  for 


q2 


-C22+V£1>C12 


-1 


C2i"zlft£l)Cll 


R(-J,)r2 


By  defining  the  reflection  matrix  at  x  =  ,  looking  to  the  left,  as 


ri£^l^ 


‘C22+zl£^l^C12 


-1 


*  A  * 

c2r^£(*i>c 


11 
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and 


#  #  #  # 
rlf(x)  =R(x-tI)rlj((f1)R<x-i1)  , 

the  previous  equation  can  be  written  as 

# 

q2  =  rH<°)r2  •  <72‘2> 


If  Eq.  (72-2)  is  substituted  back  into  the  general  solution  for  [f  j  ,  , 


1  ” 

Y(>  )  =  -  C 
'  s  y 


'#  #  # 

Clirif(x)+C12 


# 

R(-x)r2 


Similarly, 


F(x) 


#  r#  #  # 

Cf  C2irif(x)+C22 


# 

R(-x)r2  . 


If  Eqs.  (72-2),  (73-2)  and  (74-2)  are  used  and  if  the  same  mathematical 
operations  as  before  are  followed,  the  following  set  of  similar  expres¬ 
sions  can  be  obtained: 


'#  #  # 

"#  #  # 

= 

_C21+C22rit<*>. 

Cll+C12ri£M 

2)  =  total  impedance  at  x  =  £2+0  looking  to  the  left 


■  fyf2>  +  <f2> 


# 

r2£^2) 


reflection  matrix  at  x  =  £2  looking  to  the  left 


# 


r2£(x) 


r  #  A  #  i 

"C22+^2£^2^C12 


C2r^2£^2^Cll 


#  #  # 

=  R(x-£2)r2£(£2)R(x-£2) 


# 

q3  =  r2f(0)r3 
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(72-3) 


1  n 

Y(x)  =  -  C 
v  '  s  y 


U«4' 


# 

R(-x)r3 


F(x)  =  Cf 


# 

r\ 


21  zn 


# 

(x)+C 


22 


# 

R(-x)r, 


Now,  if  the  same  mathematical  operations  are  followed  again,  Eqs.  (72-3), 
(73-3)  and  (74-3)  can  be  solved;  this  must  be  repeated  until  the  solutions 
of  Eqs.  (72-(j-l)),  (73-(j-l))  and  (74-(j-l))  are  obtained.  These  are 
egressions  associated  with  waves  traveling  to  the  left. 

Similar  mathematical  operations  can  be  used  in  the  derivations  of  expres¬ 
sions  associated  with  waves  traveling  to  the  right;  however,  it  is  suggested 
that  the  right  end  of  the  rotating  shaft  will  be  considered  first.  That  is,  if 
Eqs.  (75-n),  (73-(n- 1))  and  (74-(n-l))  are  used  and  if  Eq.  (75-n)  is  substi¬ 
tuted  into  Eq.  (73-(n-l)),  and  so  on.  The  important  relationships  derived 
in  this  appendix  are  summarized  in  Chapter  2,  Eqs.  (15)-(22). 
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APPENDIX  D 


MATHEMATICAL  DERIVATIONS  FOR  IMPEDANCES 
IN  TERMS  OF  REFLECTION  MATRICES 


The  general  solution  for  [0  ,  JP  ]  can  be  written  as 

-  #  r#  #  #1# 

Y(*)  =¥Cy[Cliro(X)+C12jR(‘X>7l  (76) 

> 

~  #  r#  #  #  1  # 

F(x)  =  Cf  C2ir0(x)+C22  R(-x)7,  .  (77) 

-  J 


x  -0  x  +0 


Figure  23.  Forces  at  *  on  [o, 


Since  F(x)  is  in  the  negative  direction  at  x  =  x-0,  the  impedance 
relationship  at  this  point  can  be  expressed  as  (see  Figure  23): 

—  #  — 

-F(x)  =  sZof(x)Y(x)  (78) 

where 

# 

ZQjp(x)  =  total  impedance  at  x  =  x-0  on  [  0  ,  f  j]  looking  to  the  left. 
From  Eq.  (76), 
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#  _  r#  #  #  i , #  .« 

&("*)*!  =  s  Cnr0(x)+Cl2  Cy  Y(x) 
and  substituting  back  into  Eq.  (77), 

«  #  r#  #  #  if#  # 

F(x)  =  sCf  |C2iro<X>+C22  CliroW+CuJ  V  Y(X)  ' 

If  this  expression  is  compared  with  Eq.  (7  8),  it  can  be  concluded  that 

#  #t#  #  #  ir#  #  #  i.,#., 

Z0fW=  -Cf[C2ir0<X)+C22j[Clir0(X>+C12j  Cy 

or 

n  f#  #  #  1  f#  #  # 

WX)  C2ir0^X^+C22  ~  “  C2ir0^x)+C22_  * 


If  the  matrix  relationships  in  Appendix  B  are  used,  then 


Then,  on  [  0  ,  l  , 

Vx>  =  [4l+<=22*o(X>] 


If  similar  mathematical  manipulations  are  followed,  the  remaining 
expressions  for  total  impedance  on  each  interval,  as  listed  in  Eqs.  (23) 
in  Chapter  2,  can  also  be  obtained. 

One  may  begin  with  the  right  end  span,  i.  e. ,  on  [  ft  .  .  ft  1  ,  first,  and 

1  n- l  nJ 
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recalling  that  the  general  solution  for  ^  »  ^n]  is 


Y(x)  =  ic 


8  VL  11 


#  § 


c.,+c 


12rn(X) 


# 


F(x)  =  Cf 


#  # 

C21+C22rn<X) 


# 

RMqn 


(79) 

(80) 


l<0 


x-0  x+0 


- h 


F  (x) 


1 -  _  -  _  - 

#  , 

*»<*„> 

Figure  24.  Force  at  x  on  [ln-|»  *^n] 


Since  F(x)  is  in  the  positive  direction  at  x  =  x+0  ,  the  impedance 
relationship  at  this  point  can  be  expressed  as  follows  (see  Figure  24): 


2*  # 

F(x)  =  .Znr( 

where 

# 

Znr(x)  =  total  impedance  at  x  =  x+0  on  [^n.  j 


(81) 

Jt  ]  looking  to  the  right. 


From  Eq.  (79), 

#  r#  #  # 


# 


and  substituting  back  into  Eq.  (80), 


«  # 

■#  #  »  ■ 

-  §  #  #  ' 

F(x)  =  sCf 

C21+C22r„W 

.  Cll+C12r„<1[)_ 

C  Y(x) 
V  ' 


When  this  expression  is  compared  with  Eq.  (81),  it  can  be  concluded  that 


#  # 

r#  #  «  i 

'  #  #  #  ‘ 

Z  (x)  =  c, 
nr'  '  f 

C21+C22r„lx) 

Cll+C12rn'x> 

-1  #-l 

C 

y 


or 


U  * 

I  (x)  =  C/ Z  (x)C  = 

nr'  '  f  nr'  '  y 


f#  #  #  1 

r#  #  #  i 

C^.+C,^  (x) 
21  22  n'  ' 

_cn4Cur„W_ 

f  nr' 

Now,  the  following  equation  can  be  written: 


-1 


u 

'###' 

'  #  »  * 

°"  l'n-1  '  fnl  ’ 

l  (x)  = 
nr'  ' 

C21+C22rJX^ 

_Cll+C12r„W 

-1 


If  similar  mathematical  steps  are  followed,  the  other  expressions  of  total 
impedance  on  each  interval,  as  listed  in  Eqs.  (24)  in  Chapter  2,  also  can 
be  obtained. 
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APPENDIX  E 


MATHEMATICAL  DERIVATIONS  FOR  THE  COMPLETE 
SOLUTION  IN  MATRIX  FORM 


The  general  solution  in  Appendix  C  may  be  rewritten  in  matrix  form  as 
follows : 


On  [0, l  J 


F 


1  # 
-  C 

s  y 

# 

0 


# 

0 

# 

c. 


# 

11 

C12 

# 

21 

C2  2 

J  L 


#  # 

r0(x)  *(-*)*! 
# 

R("x)r  x 


Onl^.2.  fj.,1  • 


F 


i  #  * 

-  c  0 

s  y 


# 

0 


# 

c. 


■"11 

cu 

# 

# 

■"21 

c22 

#  # 
r(  j-2)£(x)R(‘X)rj 
# 

R(-x)r 


L 


J 


'  #  #  ' 

-  C  0 

s  y 

■# 

# 

'#  # 

O 

3 

1 

P 

Ut. 

«• 

Y 

cn 

C12 

r(j-l)f(X)R('X)rkf 

F 

#  # 

0  c£ 

# 

C21 

# 

C22 

# 

R(’x)rkf_ 

°nlajk-  d1 


On[tj.  fj+1] 


Y 

[l  # 

-  c 
s  y 

# 

0 

'# 

cn 

#  '' 
ci2 

# 

R(x)qjr 

F 

.  J 

# 

0 

# 

cd 

# 

c2i 

« 

c22. 

#  # 
r.  (x)R(x)q. 

Y 

'  !  # 

-  C 
s  y 

t 

=*=  o 

’» 

cn 

# 

C12 

- 1 

X 

Jpl 

+ 

_ 1 

Oji 

F 

# 

0 

# 

cf 

* 

C21 

# 

C22 

m 

#  # 
r(j+l)r^x^R^X^j+l 

1 17 


°ntfn-l  ’  ’ 


Y 

,  #  # 

-  c  0 

s  y 

'# 

cn 

# 

C12 

» 

R(x)% 

1 

21  Uh 

-  -  « 

#  # 

.  °  Cfj 

# 

_C21 

# 

C22. 

- 1 

lof 

=Jfc  cxj 

X 

If  the  terms  in  the  column  matrices  in  Eqs.  (82)  are  examined,  the 
following  relationships  are  obtained: 

##_#### 

=  R(x)r0(°)R{*)R(-x)r1 

#  # 

=  R(x)r0(o)7j 


#  #  #  #  # 
ru(x)R(-x)72  =  Rfx-jtjJr^U^Rtx-fjJl^-x^ 

#  #  #  # 

=  R(x)R(-f1)ru(f1)R(-i1)72 

#  # 

=  R(x)ru(0)7z 


#  #  #  # 

r(j-l)£(x)R(‘x)rj£  =  R{x)r(j-l)f(0)rjf 

§  §  #  #  #  # 
rjr(x)R(x)^jr  =  R(ij-x)r.r(i  ,)R(f  ,-x)R(x)qjr 

■  R(-x)R(«  p^WjWpqj, 

#  # 

=  R(-x)r.r(0)qjr 

#  #  #  # 
r(j+i)r(x)R(x)qj+i  =  R(-x)r(j+1)r(o)qj+1 


#  #  #  # 

rn(x)R(x)\  =  R(-x)rn(o)^  . 


If  the  above-derived  relationships  are  used,  Eqs.  (82)  may  be  rewritten 
in  the  matrix  forms  given  by  Eqs.  (26)  in  the  text. 


The  integration  constants  of  the  matrix  relationships  given  by  Eqs.  (26)  are 
determined  as  follows:  thus  far,  the  only  boundary  conditions  that  have  not 
been  used  are  those  at  x  =  a^i  t*le  locations  of  the  driving  forces.  The 

boundary  conditions  at  these  points  are  as  follows  (see  Figure  8  in  Chapter 

2): 


-0)  ■  Y(ajk+0) 


F(a.k-°)  ♦  P(a.k) 


> 


F(a.k+0)  =  0 

* 


(83) 

(84) 
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If  x  =  a.,  and  is  substituted  into  the  corresponding  general  solution,  as 

JK 

given  by  Eqs.  (82),  then 


1 

at  x  =  a.,  -0  ,  Y(a..  )  =  —  C 

jk  v  jk'  s  y 


#  # 


Clir(j-l)f(ajk)+C12j  R^"ajk)rjjf 


«  ,  #  r#  #  # 

and  at  x  =  a..  +0  ,  Y(a..  )  =  —  C 

jk  '  jk'  s  y 


1# 


C. .+C.  ,r.  (a.,  ) 

11  12  jr'  jk' 


R(aJk>qjr 


F<ajk)  =  Ci 


*#  #  # 

|  ■fC t t r.  (a.,  ) 
21  22  jr'  jk' 


R(ajk)qjr 


^  #  r#  #  #  "  ^ 

F(a.k)  -  Cf  [c2ir(j.1)j((a.k)+C12J  R{-ajk)7jf  j 


(85) 


(86) 


If  Eqs.  (85)  and  (86)  are  substituted  into  Eq.  (83),  then 


"#  #  # 

# 

#  # 

.Clir(j-l)j?(V+C12. 

R('ajk)rjr 

C.,+C.,r.  (a.j 
L  11  12  jrv  jk'J 

R(Vq 


Jr 


(87) 


Premultiplying 


#  #  ,  #  # 
i+c7|  c1,r.  (a.j 

11  12  jr'  jk' 


u» 


lJR<ajk)qjr-Lr(j-l)£<ajk)+CllCl2JR(-Vr)f  =°  •  (88) 


If  Eqs.  (85)  and  (86)  are  substituted  into  Eq.  (84),  then 

#  r#  #  #i#  * 

Cf  .C2ir(j-l)^ajk^C22jR^"ajk)rjjf  +P^ajk^  " 


Cf  |_CZl+C22rjr^ajk^J  R(ajk)qjk  "  0 


Premultiplying  , 


■#  #  # 

# 

'#  #  #  “ 

C21+C22rjr(ajk) 

*<ajk>V 

C2ir(j-l)f(ajk)+C22_ 

#  #.,* 

R('aik)7ii  =  cf  p(aik>  • 


(89) 


-1 


Premultiplying 
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“#  #  ,#  #  1  # 
_I+C21CZ2rjr^ajk^  R^ajk^qjr" 


From  Eq.  (88)  and  Eq.  (90), 


#,  #  1# 
A  — .  I  r\ 


.  *  r>/ 


r(j'  1  )^ajk)+C2 1 C2 2 J  R{"ajk)rj*  ~C2  1  Cf  P(ajk^ 

(90) 


It  may  be  verified  by  direct  matrix  algebraic  operations  that: 


#  r#  #  #  j  i  , 
C+=  C2rC22Cl2CH 


#  r#  #  ,#  #  • 

C-  =  C21C1  1 C 1 2  "  C22 


Hence, 


jf  ft  ft  ft  ft  ~ 

R(-ajk)7jj =  vVR<ajk>Vc-cf’lp(V  • 


From  Cj*  (Eq.  (87)), 


#  ,#  # 


#_,#  # 


C;2Cn+rjr(VjR(ajk)qjr-[CUClir(j-l)«(ajk)+IJrjf  *  0 


From  (Ecl-  (89)), 


"#,#  #  #  _  #,#  #  ##  #,#,“ 

_C22C21+rjr^ajk^J  R^ajk^qjr"  [,C22C2  ir(j-  l).f  ^ajk^IJR^"ajk^rj£=C22Cf  P^ajk^* 


.  y"»  *  /n  *  t-i 


From  Eq.  (93)  and  Eq.  (94), 


#  _  #  #  r#  #  #  #  ]  # * 

R^ajk^qjr"r(j-l)£^ajk^R^"ajk^rji  =  C2rC22C12Cll  Cf  P^ajk^  * 

By  using  relationship  (91)  again,  then 

R<Vv  =ru-i)j(ajk)R('ajk)7jc+*+*f  lp(V  • 


u  ~ 
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Now,  Eqs.  (92)  and  (95)  can  be  solved  simultaneously  to  yield: 


rjt =  V 


#  # 


#  #  |#_L* 


LWc++c-Jcf  P(V 


(96) 


qjr  =  R("ajk) 


‘#  #  # 


-1 


#  # 


Lr(j-1)^(ajk)C-+C+J 


cf  P(ajk)  •  (97) 


If  Eqs.  (72),  (75)  are  again  considered, 


ql.  *  ro(°)rl 


q2  =  ri((°)r2 


Vl  =  r(j-l)i(°)rj-l 


qjf  rjf(°)rjf 


r.  =  r.  (0)q. 
jr  jr'  '  jr 


■>1  =  r(j+l)r(0)qj+l 


rn- 1  -  ^(n-ljr^^n-l 

■  h°K 


(72-1) 

(72-2) 


(72-(j-l)) 

(72-j£) 

(75-jr) 

(75-0+1)) 


(75-(n-l» 

(75-n) 
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and  subsequently  combined  with  Eqs.  (73),  then 

7i  =RU1)[ciir0Ui)+c12J-1[ciiR(41)q2+Ci2R(-t1)72j  (98-i) 


(98-(j-l)) 

Vi^-v[cn^»*U*>)r«v]‘1pii^j)V=«“(-*j)7Jr. 

(98-(j+l» 

•  •  • 

%  =  R(-fn-l)[C„+*,2*„(«»-l)]‘1[cuR(«n.1)Vi+C12R(-fn-1)7n-i  . 

(98-n) 

If  Eq.  (72-2)  is  substituted  into  Eq.  (98-1),  then 

#  r#  #  #  i.i r#  #  #  #  _  i 

rl  =  R^fl)  C1  ir0^1^+C12  C1  irU^  l)+C12R^"£  l^r2 

V-  J  L 

Similarly,  if  Eq.  (72-3)  is  substituted  into  Eq.  (98-2),  then 

#  r##  #  i  r#  #  #1# 

r2  =  R^2^  _Clirif^2^+Cl2.  ,Clir2£^2^+C12.  R^"£2^r3  * 

If  these  calculations  are  continued  until  Eqs.  (72  j£)  and  (98-(j-l))  have 
been  reached,  the  complete  set  of  relationships  between  integration 
constants,  r,  is  obtained. 

If  Eq  (75-(n-l))  is  substituted  into  Eq.  (98-n),  then 

#  r#  #  #  i  , r#  #  #  i# 

%  '  R^"^n- 1  ^  Cn+Ci2rn^n-1^  _C1 1+C12r(n-  l)r^n-l^.  ^n-l^-l  * 
Similarly,  by  substituting  Eq.  (75-(n-2))  into  Eq.  (98-(n-l))  and  continuing 
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until  Eqs.  (75-jr)  and  (98-{j+l))  have  been  reached,  the  complete  set  of 
relationships  between  integration  constants,  q,  can  be  obtained. 

By  combining  the  previously  derived  relationships  with  Eqs.  (96)  and  (97), 
the  complete  set  of  relationships  is  given  as  Eqs.  (26a)  in  the  text. 
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APPENDIX  F 


MATHEMATICAL  DERIVATIONS  FOR  THE  COMPLETE 
SOLUTION  IN  TRAVELING  WAVE  FORMS 

If  the  complete  solution  in  matrix  form  is  referenced,  i.  e.  ,  Eqs.  (26)  and 
(26a)  of  Chapter  2, and  if  the  following  matrix  relationships  are  used,  the 
solution  in  traveling  wave  forms,  can  he  rewritten. 


The  relationships  used  below  are: 


(99a) 

(99b) 


#  ft 

R"1^)  =  R(-a)  .  (99c) 


The  expressions  for 
i.  e.  , 


and  q. 

V 


will  be  considered  first  (see  Eqs. 


(26a)); 


Substitution  of  relations  of  Eqs.  (22)  into  this  equation  yields 

#  r#  #  #  #  # 

7jf  =  E(a.k)[l-R(-a.k)r.r(0)r(..1)J(0)R(a.k)J-‘ 
ft  ft  ft  ft  ft  1  ft  ~ 

x[R(-ajk)R(fj)rjr(fj)R(t.-ajk)C++C.jcf-lp(ajk)  . 

#r# # #i. ,  r#  ##i 

With  relationship  (99b),  A  I-BA  =  I-AB  A  ,  and  taking 

#  #  ft  ft  ft  ft 

A  =  R(a.k;  .  B  =  R(-a.k)rjr(0)r(..1)t(0)  . 

_  #  #  Ilf##  *  *  *  Hi#,* 

VrV0)r(j-l)«<0)J  '*  [R(fj)r|r(fj)R(Jj-ajk)C++R(a.k)C.j  C£-'p(a.k)  . 

(100) 

Also  from  Eqs.  (26a), 

^  ~ft  ft  ft  ~  ~ft  ft  ft  '  ft  ^ 

vR< v^-^-w'vvvrhj- Df(ajk)c-+c+J  ■ 
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Substitution  of  Eqs.  (21)  yields 


q.  =  R( -a  ) 
IT  '  jk; 


'#  # 


Li-R(a.k)r(j.i)t(0)rjr(°)R(-a.k) 


-1 


X 


LR(a.k)R(-tj.1)^j.i)f(tj.i)R(a.k-i..1)c.+c+jc-‘p(a.k) 


#  #  I# 


With  relationship  (99b),  A 
#  # 


# 

-1 

"#  ##' 

A 

I-BA 

I-AB 

-1 


A  ,  and  taking 


A=  R(-ajk)  , 


#  #  #  # 

B  =  R(Vr(j-l)«(°)rjr(0)  • 


qjr  = 


'#  # 


-1 


I-Ib-i)«(0)rjr(0)  R<-Vi)r(j-i)fuj-i)RVfj-i)c-' 


R(-ajk)C+J 


c;'?(ajk)  . 


(101) 


If  Eqs.  (100)  and  (101)  are  substituted  into  the  last  terms  of  the  matrix 

expressions  of  Eqs.  (26)  on  [£.  .  ,  a.,  ]  and  [a.,  ,  H  .],  the  following 
n  J-1  JK  J 

manipulations  can  be  shown: 


#  # 


#  # 


R(x)r(j_i)f(°)rjf  =R(x)r(j.1)£(°)  I-rjr(0)r(j.1)|(0) 


’#  # 


-l 


#  # 


LR(fj)rjr(£j)R(f.-ajk)C++R(ajk)C.J 


# 


cf'lp(ajk)  . 


■###i,  r#  ## 

I-BA  =  I-AB 

#  #  #  # 


With  relationship  (99b),  A 


-1 


A  ,  and  taking 


A  =  W<°>- 


b  =  rjr(o) . 


#  # 


# 


R<x)r(j-l)f(°)rjf  =  R(x)LI-R<-£j-l)r(j-l)£(£j-l)R(-£i-l)rir(°)J'lR(-fj-l> 


With  relationship  (99b), 


#  #  # 

I-AB 


-1 


A  =  A 


#  #  # 

I-BA 


,  and  taking 
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#  #  #  #  #  # 

A  =  R<-lj.1).  B  =  r,^  1)£(ij.1)R!-/!j.1)r.r(0)  . 


•  /  r  ^ 

U-1M'  J-A' 


(102) 


Also, 


#  r#  #  #  #  # 


R(-x)rj£=  R(-x)|I-R(ij)rjr(fj)R(lj)r(j.lK(0)J  ‘R^r.^.jR^.-a.^ 


#  # 


#  #  # 


#  l  #-l“ 


R(-x)^I-rjr(0)r(j.1)f(0)j  ‘R(ajk)C.jCf  P(a.k)  . 

###  #  #r###i 

With  relationship  (99b),  I-AB  *A  =  A  I-BA  1  ,  and  taking 


#  # 

A  =  R(f  j)  , 


#  #  #  #  #  # 


R(.x)rj£  =  |R(-x)R(£j)LI-rjr(fj)R(£j)r(j.1);(0)R(fj)J  rjr<*.) 

#  ##("##  #  l.J  #  ]  #_,* 

XR(£.-ajk)C++R(-x)[l-r.r(0)r(j_1)f(0)J  1R(ajk)C_  jCf  AP(ajk)  . 


##,###•#, 

With  relationship  (99a),  I-A  =  I  +  A  I  -  A  ,  and  taking 

^  =  V°ViK(0)' 

I  r'  #  r#  #  ft  1##  ##  # 

R(-x)7jt=  mi,-*)  J-Vtpry.^dp  •1rjr(£j)R(«ra.k)C++R(ajk-x)C.+ 


#  # 


#  # 


R(-x)R(fj)rjr(fj)R(£j)r(j.1)£(0)Li-R(9rjr(£.)R(fj)r(j.1)£(0) 


#  #  #  #  # 


#  #  If.!* 

XR(a  )C_  c/^a  ) 


#  #  #  # 


With  relationship  (99b),  A  I-BA  =  I-AB  A  ,  and  taking 


###.,# 
T  A  r>  *  A 
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A=r..(fj)i(ij)r(j-1)((°).  b=r((j), 


# 


f#  # 
,-x)  I  i-r. 


# 


"i «  LI-V  Wj.  i)4  ( VJ ' “WiW 


-1# 
*  T** 


r# 


# 


# 


# 


R(fj-^j-l)r(j-l)^(fj-l)R(ajk-pj-l>c-J  +R(ajk— )c-  Jcf  p(-jk)  •  (103) 


*  ^-1“ 


Also, 

#  # 
R(x)qjr=R(x) 


#  # 


# 


rr(j-l)«(0)rjr(°)|  r(j-l)f(0)R(ajJC-  + 

# 


.  #  # 
R(-VS 

With  relationship  (99a), 

(  n 

,  >  =<fe . 

-1 

#  # 
=  I  +  A 

"#  #“ 
I- A 

.-1 


-1 


,  and  taking 


#  # 
R(x)qjr=R(x) 


#  # 


- 1 


#  # 


LI-I(j-l)f(0)rjr<0>J  r(j-l)^0Wajk)C-+R<“ajk)C++ 


#  # 


r(M)l«0)Ijr<0)LI"r(j-l)l<0)V°)J  R^ajk)C+  Cf  P(ajk) 


#  '#-1* 


With  relationship  (99b),  A 

#  #  # 
A’r(j-l)l(0)Ijr(0) 


# 

"#  #  #1 

- 1 

#  #' 

A 

i-baJ 

I-AB 

-1 


A  ,  and  taking 


R(*)qjr=<i  R(x) 


'#  #  # 
I-r(j-l)l(0,rjr<0> 


-1 


a  it 

B  =  I  , 

*u-o«<0) 


R(t,)r 


# 


R(a.k)C-  + 


-1. 


+R(x-ajk)c+j.Cf‘p(a.k)  . 


jR' 


By  following  the  same  matrix  algebraic  operations  as  those  used  in  the 
evaluations  of  R(-x)r^  ,  the  preceding  equation  can  be  transformed  to: 
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R(*)V  R(x-fj.1)[”(j.1)^.-l)yf..1)j-‘r0.1)f(Cj.1) 


n-l» 


x[R(a.k-fj.1)C.+R({.-i..1)nr(fj)R(tj-a.k)C+ 

#  #  ]  #  ~ 

R(x-a  >c+  c£-‘^a  >  . 


H- 


(104) 


Also, 


# 


R(-x)rjr(o)qjr=R(.x)R(fj)rjr(^j)R(^)^r(j_1)f(o)R(iprjr(Cj)R(fpj 


'#  # 


#  # 


-l 


#  #  #  # 


*  1  !-i: 


With 


x|r(j.i)((o)R(ajk)c.+R{-aik)c+|cf  P(ajk)  . 

#  r#  #  #*]  ,  rii  #n  .# 

relationship  (99h),  A  Ll“BAj  -  lI-ABJ  A  ,  and  taking 


# 


# 


A  =  rjr(f j)R(Cj)  ■ 


R(-x)rjr(0)qjr=R(t.-x) 


#  # 


# 


i-r  iJfjr, 


OU 


jr'  y  (j~  i)*  r J 


-l 


)R(a  -I  C-‘P(a.k)  .  (105) 


*  I  #.-l“ 


Eqs.  (102),  ( 1 G 3 ) ,  (104),  and  (105)  are  the  complex  entries  of  the  right-hand 
column  matrix  of  Eqs.  (26), 


°n  i { j- 1  •  ajJ  • 


■  ■ 

1  # 

# 

# 

# 

’#  # 

Y 

-  c 

s  y 

0 

C11 

C12 

R(X)r(j-l)f(°)rjf 

“ 

# 

# 

# 

# 

F 

0 

cf 

C2  1 

C22 

R(-x)T.t  ^ 

On[a.k.  1.]  . 
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r 

1  # 

-  c 
s  y 

#  ' 

# 

# 

Y 

0 

cu 

C12 

R(x)qjr 

rsj 

# 

# 

# 

# 

#  # 

F 

0 

cf 

CZ1 

C22 

R(-x)rjr(0)qjr 

Also,  if  the  identity  of  the  various  parts  of  Eqs.  (102),  (103),  (104),  and 
(105)  is  observed,  the  following  relationships  can  be  obtained: 


#  #  #  # 

R(-x)7jjf  =  R(-x)r.r(0)q.r+R(a.k-x)C_Cf_1P(a.k) 


»  *.-i- 


#  #  ,«* 
>•  1  < 


#  #  #  # 

R(*)qjr  =  R(x)r(j-i)£<°)7ji  +R<x-ajk)c+cf‘p(ajk) 


(106) 

(107) 


From  the  above-derived  equations,  some  useful  functions  related  to  the 

concept  of  traveling  waves  can  be  defined.  On  [j£  ,  J?.]  ,  the  incident 

J  J 

£3t 

waves,  caused  by  the  single  driving  force  P(a.,),  traveling  to  the  right 

JK 

and  to  the  left,  respectively,  are  as  follows: 

>  # 


Ujjk(x)  -  R(x"ajk)C+Cf  P(ajk)[  H(x'aj)()'H(x'{j)) 

exists  only  on  [  a^k  ,  i 

<  # 


(108) 


*  *.-1“ 


Ujjk(x> =  R<ajk‘x>c-c£  p(ajk)tH(ajk'x)'H(fj-rx>] 

exists  only  on  [^j_j  »  ajjJ 


(109) 


If  the  expressions  for  incident  waves  are  used,  Eqs.  (102),  (103),  (104), 
and  (105)  can  be  rewritten  as  follows: 


°nlfj-i-  ajki 


'#  # 


-i, 


R(x>r( j- l)£(°)rjf  "R(x"f  j- 1>  LI_r(j- 1  )f(® j-l)rjr(fij- l^J  r(j-l)£({j-l) 

xVniVWVi-.» 


(110) 
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# 

R(’x)rjje=R^j"x) 


‘-Vv'u-dW. 


-1 


r.  (£.) 

jr'  y 


.VWhWh'V'j-i1 


+v(x) 


(111) 


on  f  ajk’  fjl  • 
#  * 


'#  * 


-1 


R(*)qjr!SR(*-«j.1)[i-r(j.1)£(ij.1)rjr(i..i)J  V^.,) 

x  i>  J + v*> 

#  #  #  r*  #  #  ■■ 


(112) 


R(-x)rjr(0)q.r=R(fj-x)  [I‘rjr(tj)r(j.1)f(fj)J  VV 

x[vw-|j-i4.i)i«j-i)VV,!I  • 


(113) 


From  an  examination  of  the  above  expressions,  the  reflected  waves, 
originated  from  the  single  driving  force  P(a.,  )  ,  traveling  to  the  right  and 

JK 

to  the  left,  respectively,  can  be  defined  as  follows: 

LI-rO-i)f(ej-i)V{i-i)J'1*(J-i)«(fri) 


[R(«j-fj.1)rjrUj)0jjk(«j)+4k(«j.i)JlH(x-{j.i)-H(x.£.)J 


(114) 


« 


#  # 


°JjW)  ^jr^j^-l^V 


-1 


r.  (£.) 

jrN  J7 


r>  #  #  _< 

u..,(£.)+R(£  •-£•  ,)r/.  lX,(£.  ,)u..,(£.  .) 

L  JJk  J  J  J’1  (j-lK  J_1/  JJk  3-1  J 


H(x-l._1)-H(x-4j) 


115) 


both  exist  only  on  [£^  ^  ,  ft 


When  the  expressions  for  reflected  waves  are  used,  Eqs.  (110),  (111), 
(112),  and  (1 13)  can  be  rewritten  in  much  more  compact  forms  as  follows : 
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°n[£j-l-  ajkl  • 


I  »  » 


R(*)r(j.i)J(0)r.f  =  ujjk(x) 

R(-x)7jt  =  Ujjk(x)  +  Ujjk(x)  . 

On  [a.k,  (.]  . 

»  »  > 


Rtx)q.r  =  Ujjk(x)  +  U.jk(x) 

#  #  « 


(116) 

(117) 


(118) 


R(-x)rjr(0)qjr  =  ujjk(x)  . 


(119) 


Now,  if  the  above  expressions  are  substituted  back  into  the  matrix  form, 

the  dynamic  response  on  [  ft .  ^  ,  f  .]  can  be  expressed  as  follows: 

J  J 


°n  [  *  j.  j  .  *jk]  • 


,  #  # 

— C  0 

s  y 

# 

# 

» 

Y 

1 

C11 

C12 

V(x) 

cy 

#  # 

» 

# 

«  < 

F 

0  Cf 

C2 1 

C22 

UijkW+v'*’ 

On[a.k,  i.)  , 


>v 

Y 

fi# 

~C 
s  y 

O 

'# 

cn 

1 

C12 

'  »  > 

TF.  (x)+TT. (x) 
jjk'  '  jjk'  ' 

F 

3 

# 

0 

# 

cf 

# 

C21 

# 

C22 

i 

"ST 

Vp^ 

i 

From  the  following  observations, 

> 

Ujjk(x)  exists  only  on  [a  k>  f .  ] 
< 


and 


Tfjjk(x) 


exists  only  on 


t(i-> 


(120) 


(121) 


132 


additional  terms  may  be  included  in  the  preceding  matrix  expressions 
without  changing  their  original  nature,  as  follows: 


°"  Uj-,.  ajk]  • 


r,# 

#  “ 

■# 

» 

>>  > 

Y 

i-C 
s  y 

0 

C11 

C12 

V(x)+V(x) 

a* 

= 

# 

# 

# 

# 

«  < 

F 

0 

cf 

i 

C21 

C22 

Vx)+u»*(It) 

°n  tajk’  V  • 


# 

#  ‘ 

"# 

If 

»  > 

Y 

i-c 

s  y 

0 

C11 

C12 

Ujjk(x)4U.jk(x) 

22 

# 

# 

f 

* 

«  < 

F 

0 

cf_ 

C21 

C22 

U.jk(x)lTJ .  jk(x) 

It  is  obvious  that  the  above  two  matrix  equations  are  completely 
identical;  hence,  they  may  be  combined  as  follows: 


(122) 


(123) 


On[tj.,.fj) 


Cs£ 

\  #  # ' 

•— C  o 

s  y 

’  # 

# 

»  > 

Y 

C11 

C12 

V(x)+Ujjk(x> 

#  # 

# 

# 

«  < 

F 

0  Cf 

1 

o 

C22 

Ujjk(x)+Ujjk<x) 

(124) 


This  matrix  equation  suggests  that  the  total  traveling  waves,  originated 

from  the  single  driving  force  P(ajjc)  »  traveling  to  the  right  and  to  the  left, 
respectively,  are  as  follows: 


»>  »  > 


'°jjk(x)  =  Ujjk<x>  +  Ujjk<x> 

(125) 

and 

«<  «  < 

vm{x)  °  ujjk(x)  +  ujjk<x)  • 

(126) 

Hence,  the  final  expression  of  the  dynamic  response  caused  by  the  single 
driving  force  P(a.^)  on  [£  ._^  ,  £  .]  is  as  follows: 
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(12?) 


# 

# 

# 

Y 

ic 

s  y 

0 

C11 

C12 

Jjk'  ' 

& 

# 

# 

# 

# 

<« 

F 

0 

cf 

(M 

u 

_ 1 

C22 

i _ 

If  the  superposition  technique  is  used,  the  dynamic  response  on[f . 


j-1 


caused  by  all  the  driving  forces, 
span  is  as  follows: 


P(ajk)  ,  k  =  1,  2,  . 


k^j)  ,  in  this 


Incident  waves  on  [  6.  ^  ,  £ .  ]  : 

J  .  J 

>  k(j) 

R(x-ajk)C+^- 1  S(ajk)  [H(x-ajk)-H(x-f.)| 

vx)  ■  E 

k=  1 

(128) 

and 

<  «i> 

R(ajk-x)C.Cf1P(ajk)[H(ajk-x)-H(lj.1-x)]  . 

Vx>  =  E 

k=  1 

(129) 

Reflected  waves  on 


» 


» 


KJ) 

vx)=E  vx) 


k=  1 


and 


« 

*0) 

« 

U..(x) 
Jj'  ' 

II 

M 

Vx) 

• 

k=  1 

Total 

wave  8 

on 

fjl : 

»> 

KJ) 

>i» 

Mi)  »  > 

TT..(x) 
JJV  ' 

=  J] 

J* 

II 

E  Ujjk(x)+TTjjk<x> 

k=l 

k=l 

and 


(130) 


(131) 


(132) 


134 


(133) 


«<  *(1)  <«  «  < 

uu(x)  *  E  v(x)=I  v<x)+v(x)  • 

k=l  k-1 

The  solution  of  the  dynamic  response  on  [jL_^  »  £j]  is  as  follows: 


Y 

=* 

F 


1# 

•ic 

s  y 

# 

0 


# 

0 

# 

c. 


# 

c 


11 


# 

c 


12 


# 

c 


21 


# 

C 


22 


JL 


l^> 

U..(x) 
JJ  7 

«< 

U..(x) 

jy  ’ 


(134) 


Now,  the  traveling  wave  forms  in  other  spans  of  the  shaft  will  be 
investigated. 


On  [*j_2»  *  j-il  (8ee  E<18-  (26)) 


2* 

Y 

Oi 

F 


1# 
i-C 
a  y 

# 

0 


# 

0 

# 

c. 


# 

c 

# 

c 


11 


21 


22 


C12  R(x)r(j-2)f(°)rj-l 

# 

C 


R(’x)rj.l 


From  Eq.  (26a),  the  following  equation  can  be  obtained  by  substitution: 

##  ##  r  #  r##  #1 

R(x)r(j-2)i(0)'j-i=R(lt)ro-2)£(0)|R(tj-i)[ciir(j-2)i(1j-i)+c1aJ"1 

“#  #  #  1  #  _  I 

X  [Clir{j-l)f^j-P+C12J  R<'Vl)rjf  j  * 

From  Eq.  (82), 


R(x)r(j-2)£(0)rj-l=R<x‘£j-2)r(j-2)£^j-2^fj-rfj-2) 


■#  #  # 

.  1 

'#  #  #  ■ 

X 

Clir(j-2)f^j-l^+C12 

•  1 

_Clir(j-l)£^j-l^+C12_ 

X  s 


#  #  # 


LCHr(j-l)£(fj-l)+Cl2J  Cy  Y(fM+0)  • 


-ll-l 
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#_1#  #.# 

From  Eq.  (134),  and  since  » 

#R<X)*(j-2)l(°)7j-l 

=R(-‘j.2)r(j.,){(<j.2)R(«j.1-<j.2)[Hc;;c12r(j.2)<(«j.1) 

#  ,  #  »>  «< 

C,  c17TT..(f.  ,)+Tr..(£.  ,) 

11  12  jj'  j-1'  jj'  j-1' 


-1 


X 


The  total  wave  traveling  to  the  right  on  [i.  ,  ,  f  .  ,]  ,  i.  e.  ,  the  (j-l)*^ 

j“£  j-i 

span,  can  be  defined  as  follows  : 


»> 


VOjW  =  R(X)*(j-2)«<°>rj-l 


or 


-R(x.|j.2)r(w)|(Ij.2)R(lj.1-Ij.2) 


'##  ,#  # 

I+Cll  C12r(j-2)f(fj-l) 


«K 


#  ,  #  >;>e>  ^ 
cnci2trjj(‘j.-i)+Vtj-i) 


-1 


(135) 


Similarly,  the  total  wave  traveling  to  the  left  on  (  i  .  ,  ,  i ,  .  ]  ,  i.  e.  ,  the 
(j-1)  span,  can  be  obtained  as  follows: 


«< 


TTu-i)j(x)  =  R(‘x)rj-i 

#  f#  #  ,#  # 


2^0-2, |(^.1)J"1[cnC,2V*j.1)«jj(|j.i) 


• (136) 


Hence,  the  dynamic  response  on  [jE.  ,,  £  }  ,  i.  e.,  the  (j-1)  span, 

J  J  * 

•2* 

caused  by  driving  forces  P(a.J  ,  k  =  1,  2,  .  .  .  ,  k(j)  ,  on  [  f .  ,  f .] 

th  J  J"1  J 


th 

i.  e.  ,  the  j  span,  can  be  expressed  by  the  following  matrix  equation: 

» 

11  ''1211  J(j-l)j( 


■U: 

#  ‘ 

Y 

0 

_ 

b  y 

# 

# 

F 

0 

Cfj 

. 

u 

C 


# 

c 


21 


# 

C 


22 


«< 

U 


(j-l)j 


.(*> 


(137) 


136 


With  the  same  mathematical  manipulations,  the  following  set  of  matrix 
equations  can  be  obtained: 


th 


0n  U  ;_3  •  &  :_2]  »  i-  e-  *  the  (j-2)  span, 

J  J 


Y 

F 


1# 
-iC 
s  y 

# 

0 


where 

«< 


# 


# 

0 

# 

c. 


11 


21 


# 

C 

# 

c 


12 


22 


»> 

«< 

Uu-i)/x) 


u(j.2)j<x)=R(« j.2-x)t  i+c,  ;cl2r(j.3)f(« j.2> 


-1 


’#,#>»>  <<< 


On 


[  ^ j-4  *  ^  j«3  ]  »  i-e.,  the  (j-3)th  span, 


U^j(x) «  R(« j-x)  I+Cj [c,2r2f (f  3) 


r#  #  »>  «< 

C  1 1 C1 2  3^+U4/f  3^_ 


On  [ij  ,  i2]  »  i-e.  ,  the  2nd  span, 


Y 

fl# 

ic 
s  y 

# 

0 

'l 

C11 

# 

C12 

»> 

U2j<X> 

F 

# 

0 

# 

cf 

» 

C21 

# 

C22 

«< 

U2j(x)  . 

where 


«<  #  r#  #  , 

U2j(x)=R(*2-x)  I+c;} 


#  # 

C12rif<f2) 


-1 


#  ,  #  >*•  «< 
CnC12U3/f2>+U3/t2). 


On  [  0  ,  J?  j  ]  ,  i.  e. 


the  l8t 


span, 
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where 

«<  #  r#  # ,  #  #  i  ,  r#  ,  #  «< 

Similar  mathematical  manipulations  can  be  applied  to  other  portions  of  the 
shaft. 

On  [  JP j  »  £j_  J  (see  Eq.  (26)), 


From  Eq.  (26a),  the  following  equation  can  be  obtained  by  substitution: 

#  #  r#  r#  #  #  i  . 

R(x)q.+  1=R(x)|R(-t.)^Cu+C12r(j+1)r(£j)J-1 

~U  #  U  1  # 

x  c  +c  2r.  (It.)  R(£.)q-  f  . 

L  11  12  jr'  j'J  v  j/njr J 

From  Eqs.  (82), 

#  #  r###  ir### 

r#  #  # 

x  C,,+Ci:>r.  (£.)  AC  1  Y(l  .-0)  . 

11  12  jr'  j'J  y  '  j  ' 

#_1#  #.j# 

From  Eq.  (134),  and  since  =  Cj  jC^  » 

R(x)q.+  1=R(x-lj)[cu+C12^+^dp]-1[^.(£.)4;^12^j(f^  . 


138 


The  total  wave  traveling  to  the  right  on  [j^  *  *•  e*  »  the  (j+l)**1  span, 

can  be  defined  as  follows: 


»>  # 

■°(j+i)j(x)  =  R(x>qj+i 


f 

=R(x-{j) 


'»  »  » 

Cll+C12r(j+l)r<V 


’>»  #  *  «< 
,V'j)+ciici2ujj({j) 


(138) 


Similarly,  the  total  wave  traveling  to  the  left  on  [jf .  ,  H  .]  ,  i.  e. ,  the 

fV.  J  J1"1 


(j+l)th  span,  is 

«< 


Vi)/x) =  R<-  4d-(#iVi 


r#  #  ,  #  # 


=R(«j+i-)r,j+i)r(|j+1,R((j+i-fJ)[i+c,;c,2r(j+1)r(r) 

'»>  #  ,#  «< 

xLTrjj(lj)+0iici2TfJj(|j)J  • 


-1 


(139) 


th 


Hence,  the  dynamic  response  on  [  JL  ,  ^ j+j]  »  i.  e.  ,  the  (j+1)  span, 

r>j 

caused  by  driving  forces  P(a ..  )  ,  k  =  1 ,  2,  .  .  .  ,  k(j)  ,  on  [  f .  ,  £ .  ]  , 

th  J  J 

i.  e.  ,  the  j  span,  can  be  expressed  by  the  following  matrix  equation: 


Oi 

Y 

fl# 

s  y 

#  ' 
0 

'# 

C11 

# 

C12 

u. 

+ 

*—• 

w. 

_ i 

F 

# 

0 

# 

cf 

# 

C21 

# 

C22 

<« 

V>W<x) 

.  J 

(140) 


With  similar  mathematical  manipulations,  the  following  set  of  matrix 
equations  can  be  obtained: 

On  [  f  j+1  ,  l.+z  ]  ,  i.  e.  ,  the  (j+2)th  span, 


a  >< 

I _ 

1  # 

s  y 

#  ' 
0 

’# 

C11 

# 

C12 

»> 

V2)/50 

a  U4 

- 1 

11 

# 

0 

# 

cf 

# 

C21 

# 

C22 

<« 

UU42)/X> 
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where 


<<<  #  #  # 

U(j+2)/X)=R^j+2"x)r(j+2)r^j+2)R^j+2"fj+l^ 


X 


#  jf  ft  “|  P»>  ft  ft  «< 

I+CllC12r(j+2)r(^j+l)  1  U(j+l)j^j+l)+  CllC12tT(j+l)j^j+l) 


th 


On  [  ,  f  j+3  ]  ,  i.  e.  r  the  (j+3)  span. 


«< 


"#  #  ,# 


-iP# 


»  .  #  «< 

^  *  ▼ 


LI+CllC12r(n-2)r^n-3^J  [_  U(n-3)j^n- 3^+Cl  lCl2U(n-3)j^n-3^J 


On[j^n_2,  fn-1]  ,  i.e.,  the  (n-l)th  span, 


r-J 

,  #  # 

•i-C  0 

s  y 

# 

# 

»> 

Y 

cu 

C12 

D(n-l)iW 

#  # 

# 

# 

«< 

F 

1 

- — 1 

U 

o 

C21 

C22 

U,  .  .  .(x) 
(n-l)j'  ' 

where 


«<  #  #  # 

U,  ,  N  .(x)=R(f  ,-x)r,  |  \  (f  ,)R(f  ,) 

(n-l)j'  7  'n-1  '  (n-l)r'  n-1'  '  n-1  n-2' 


X 


#  #  ,  #  # 

i+c‘  c,,r,  n  (f  ,) 

11  12  (n-  l)r'  n-2' 


>>>  #  .#  «< 

u,  ?\-U  ,)+c7  c,  tj,  ,) 

(n-Z)j'  n-2'  11  12  (n-2)jl  n-2' 


On  [  i  ^  ,  £n  ]  ,  i.  e.  ,  the  n^1  span. 


>v 

Y 

M 

s  y 

# 

0 

- 

# 

cn 

# 

C12 

»> 

U  .(x) 
nj'  ' 

F 

O  =*- 

# 

cf_ 

» 

C2 1 

f 

C22. 

«< 

TJ  (x) 

L  J 
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where 


^  ft  ft 

VVnWn-l) 


#  #  ,  #  # 
i+c:;c.,r  (i  ) 

11  1 2  n'  n' 


-1 


X 


»> 


#  ,#  <<< 

•  *  r  r 


U,  ,v(*  ,)+c,  C^U,  lX.(f  .) 

(n-l)j'  n-1'  11  12  (n-l)j'  n-1' 


Again,  if  superposition  techniques  are  applied,  the  total  dynamic  response 


on  any  span  due  to  all  the  driving  forces  P(a.,  )  ,  j  =  1,  2,  .  . 
k  =  1,  2,  .  .  .  ,  k(j)  ,  in  traveling  wave  forms  is  as  follows: 


,  n  ,  and 


th 

On  [^_i  »  ft  ,  i.  e.  ,  the  i  span. 


Y 

1  * 
--C 
s  y 

# 

0 

’# 

C11 

# 

C12 

- 

»> 

F 

# 

0 

# 

cf 

# 

C21 

# 

C22j 

«< 

n;(x) 

where 


j=l 


(141) 


Eq.  (141)  is  written  out  in  detail  in  Chapter  2. 

>  <  »  «  »>  «< 

The  subscripts,  i  and  j  ,  for  U. . ,  U. . ,  U. .  ,  U. .  ,  U. .  ,  U. .  denote  that: 

^  J  ij  ij  ij  ij  ij  ij 


If  i  =  j  ,  the  wave  is  originated  from  the  driving  forces  in  the  same 
span  as  where  the  wave  is  located. 

If  i  ^  j  ,  the  wave  is  originated  from  the  driving  forces  in  other  spans 
and  is  modified  when  the  wave  passes  through  supports  into 
the  span  under  investigation. 
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APPENDIX  G 


EVALUATION  OF  THE  REFLECTION  MATRIX  OF  A  FIXED  END 
For  a  fixed  end,  at  x  =  0-0  ,  the  boundary  condition  is 

at 

Y  =  0  . 

Hence,  from  Eq.  (10), 


“  1  # 
Y(0)  =  i  Cy 


'#  #  #  # 
CjjRfOjqj+C^Rf-O^ 


=  0 


or  #  ,  # 

q,  ■  (-cTtc 


ll“12)rl 


(142) 


For  an  end  support  which  is  not  specialized,  it  was  shown  in  Appendix  C 
that 

# 

qi  =  r0(0)7l  • 

Comparing  this  expression  with  Eq.  (142),  it  can  be  concluded  that  for  a 
fixed  end  at  x  =  0,  the  reflection  matrix  is 

#  #  .# 

r0<°>  =  -C11C12  • 


Similarly,  it  can  b  '  proved  also  that 

#  #  ,# 

r  (f  )  =  -cTJc,, 

n'  n7  11  12 


for  a  fixed  end  at  x  =  f  . 
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APPENDIX  H 


NOTES  ON  SUPPORT  IMPEDANCES 

#  # 

To  simplify  the  calculation  of  end  impedance  Zg(0)  or  ZR(£n)  »  where 
n  =  2  for  the  shaft  system  shown  in  Figure  16,  it  can  be  assumed  that  Aq 
and  Bn  are  rigid  bodies  (see  Figure  1).  It  has  been  shown  (reference  14) 

#  ^ 

that  if  Z  (x.)  is  the  j  concentrated  impedance,  where  f  <  x.  <  b  ,  the 
c  J  ^  ft  J 

concentrated  right  end  impedance  Zn(jf n)  may  be  calculated  by  the  following 

relationship: 

ff  ^  ^ 

Wl!  S,<V’‘j>Zc<*j>S<fn-xj>  •  xj  IV  bl  <143> 

j 

# 

Similarly,  the  concentrated  left  end  impedance  Zq(0)  may  be  calculated 
from 

^  t  If  n  a 

Z0iO)*YiS'(-*JZei*k)S(-*k).  xk  on  [  a  ,  0]  ,  (144) 

k 


# 


th 


where  Zc(xk)  is  the  k  concentrated  impedance,  a<  x^  <  0  ,  and 

1 
0 


# 

s(\)  = 


Y  is  real  , 


#  # 

and  S'(y)  is  the  transpose  of  S(y)  . 

#  4H 

If  both  ends  have  the  same  configuration,  i.  e.  ,  Zc(x^)  =  Zc(£n+  |x^  |  )  , 
where  a  <  <  0  ,  Eqs.  (143)  and  (144)  can  be  used  to  calculate  the  con¬ 

centrated  end  impedances  as  follows: 


Let  Zc(xk)=Zc({n+|xk|)= 


‘cl  1 


c22 


a  <  \  <  0 


If  Eq.  (143)  is  used  and  if  j x,  I  =  y  ,  since  x,  is  a  negative  value,  then 
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U  M  it  # 

Zn<£n>  =  S'(-V)ZCZ(-V) 


From  Eq.  (144), 


#  #  #  # 
Z0(0)  =S'(v)ZcS(v) 


If  Eq.  (145)  is  compared  with  Eq.  (146),  it  can  be  concluded  that 

A  # 

Z0<°>  *  Zn«n> 


(145) 


(146) 


(147) 
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provided  both  ends  have  the  same  configuration. 
With  reference  to  Eqs.  (57), 


rlt<«l)  - 

ff  A  If 

_‘C22+Z  ljf^l^C12. 

-1 

C21-l^.)ClJ 

J1  =•* 

1— • 

II 

"  #  #  " 
_"C22+^lr^l'C12. 

-1 

p21-,l*<il)C|l] 

A 
=  z, 


1<*1 


)  +  * 


Of^i) 


V*l> 

y*i> 


=  fi(ti)  +  i2r(ii) 


"#  #  # 

■#  #  # 

9 

-1 

“ 

C2l+C22r0^1^ 

cn+ci2ro<{i>. 

"###' 

'#  #  # 

-1 

4 

— 

_C21+C22r2^1^. 

[cn+ci2r2(fi)J 

#  #  #  # 

r yij)  =  R(*1)r0(o)R(jl1) 

#  #  #  # 

r2U!)  =Rd2-f1)r2(fi2)Rd2-£1) 


and 


# 

r0<°>  = 


#  ^ 

"C22+Z0^C12 


r#  a  #  i 
LC2l“Z0<°>CllJ 


# 

r2^2^ 


‘  #  a  # 
■C22+Z2^2^C12 


-1 


r#  a  # ' 

C2rz2^2^Cll 


> 


Eq.  (147)  implies  that  £n(0)  =  t2(£2)  ’  llence>  Eq-  (152) 


gives 


#  # 

r0(o)  =  r2(£2)  . 

If  and  if  Eq.  (153)  is  used,  then,  from  Eq,  (151), 


(148) 


(149) 


(150) 


(151) 


(152) 


(153) 
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#  # 
r0(*i)  =  r2Ui)  , 

which  also  implies  (see  Eq.  (150))  that 

Z0jf^l)  =  ^2r^l)  * 


# 


From  examination  of  the  matrix  form  for 

i.  e. 

•  # 

t 

zl<cl>^ 

zl(ll) 

0 

» 

0 

Zl<22) 

it  can  be  shown  that 

ii 

rH 

<]  nT 

1  0 

zl(ll) 

0 

1  0 

Zl(  1 1)  ° 

0  -1 

0 

Zl{22) 

0  -1 

°  Zl(22) 

or 

A  # 

zi<*i>  ■ 


which  implies  that  ^(f  ^)  =  ij(f  ^) 


From  the  above  relationships,  it  can  be  concluded  from  Eqs.  (149)  and  (148) 
that 

#  # 
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APPENDIX  I 


MATHEMATICAL  DERIVATIONS  FOR  THE  QUASI- MATCHED 
IMPEDANCE  AT  THE  INTERIOR  SUPPORT 

First,  all  related  equations  pertaining  to  the  upcoming  derivations  can  be 


summarized  as  follows  (see  Eqs.  (57)): 

#  r#  #~ir#  #’ 

rlr(fi)“  'C22+^ir(q)C12  C2l‘^lr(f  l)Cl  1  (154) 

Jlr<q)^l(fl)  +  <155> 

~ q  #  §  ir#  ft  ft  1 

^2r^l)=  C21+C22I2^1^  Cll+C12r2^1^  (15^) 

ft  #  ft  # 

r2Ui)=R(*2-fi)r2(*2)R(Vpi>  O^7) 

ft  I"  #  „  #  1  .r#  u  # 

r2^2^=  _’C22+^2^2^C12  C2 1  ‘^2^2^C1 1_  (158) 


where 

#  fZ2( 11)  Z2( 1 2) 

z2(^2)  = 

Z2(2 1 )  Z2(2  2) 

ft  r2^2hl  r2^2^l2 

r2^2^  = 

_r2^2^21  r2^2^22 


If  direct  matrix  algebraic  operations  are  used  in  Eq.  (158), 

ft 

r2^2^ 

e2e3_Xe2Z2(21)+ie3Z2(22)  "  1+ie2Z2(  1  l)"le3  Z2(  1 2) 

X 

el+iele3Z2(21)  +  Z2(22)  *le3"iele3Z2(ll)‘Z2(12) 
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X 


-ie3+ie1e3z2(1  ^“^(U)  1  +  ie2Z2(l  l)+ie3Z2(12) 


el+iele3Z2(21)_Z2(22) 


"e2e3+ie2Z2(21)+ie3Z2(22) 


Examination  rf  Eqs.  (156)  and  (157)  shows  that 

#  #  #  # 
r2^1^  =  R^2-^  l^r2^  2^R^2"^  1^ 


-ie^(f2-j?1)  -ie1^(f2-f1) 

;  r2^2hle 


“  i e  j  ( ■£  2  ”  ^  1  ^ 


r2^2)2l' 


~eZ^  ^  2~^1^ 


-e?NC  (^2_^i)  (^?""^i) 

1  r2«2>.2* 

-e2*X7(£2-£1)  -e2^(£2-£1) 

:  r2U2J22e 


(159) 


The  only  element  not  being  modified  by  an  exponential  decaying  function  in 
the  above  matrix  is  the  one  in  the  first  row  and  the  first  column.  For 
qua  si- matching,  only  the  element  not  being  modified  by  an  exponential 
decaying  term  needs  to  be  considered;  i.  e.  , 


#  #  #  # 

r2(£1)  =  R({2-{1)r2(f2)R(£i-«1) 


W 


-2ie1^7  U2"£  j) 


2'x  2'  1 1 


0 

0 


From  Eq.  (159). 


r2^2^lf 


1 


det 


1  I" 

"C22+^2^2^C12 


r[(I+e3)( 


z2(22)+lele2z2(ll))  + 


2.. 


(ele3+ie2^z2(12)z2(2l)"z2(ll)z2(22)^  ” 
e3(e2+ie1)(Z2(lZ)  +  z2(21)>-<el+ieZe3>]  • 


(160) 
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where 


det 


-c„+L(f,)c, 


22  2V  2'  12J 

=  "(le3+ie1e3z2(ll)+z2(12)^e2e3"ie2Z2(21)+ie3Z2(22)^  + 

(e1+ie1e3z2(21)+z2(22)^(l“ie2z2(1 1)+le3z2(12))  * 

With  reference  to  Eq.  (156), 


'  #  #  # 
C21+C22r2^l) 


■ 

■ 

" 

‘ 

"  ie3 

i 

+ 

ie3 

-1 

.  el 

‘e2e3. 

.  el 

”e2e3 

-Zie^sTu)  (^2“^  1^ 


r2^  *2^ 1 1 


ie. 


-Zie^sfu  ( £2-£  j, ) 


r2^2hl  "  1 


■Zie^'lu  (f  2-f  j) 


r2^2^11  +  1 


e2e3 


Similarly, 


#  #  # 

LGll+C12r2<£l)J 


* 

’xele3 

-ie2 

+ 

"iele3 

”ie2 

1 

"ie3  j 

-1 

ie3 

-2ie1*/w  (f  2-f  j) 


r2<£2>ll 


~lele3 


-2ie.*/w  (i9-S. ,) 

(•  r2<*2>ll  +  1 

-ie2 

/  -2ie  (£2-£.) 

(e 

-u3 

or 
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det 


WuW 


ele3 


I  -tie^  »;>-*!> 


r2^2^  1 1  +  1 


ie. 


r2^2hl  "  1 


then, 


Wuiu.i 


det 


n — n — : 

Cll  +  C12r2^l) 


-ie. 


ie- 


\e 


-2ie j*/ui  (f 2”£ j) 


-ie  j  e 


-2ie^7cj  (J^2-f  ^) 


Wu+l 


Hence, 


*2^l)  = 


det 


ri — FI — : 

Cll+C12r2^l) 


(1+®?) 


-2ie1's/u)  (^2-fj) 


r2<*2>ll  “  1 


“2ie1N/T3(i2-f }) 

L-e3(e2+ie1)e  r2^2^11  +  e3^e2”iel^ 

-2ie1^{je2-£1) 

e3(e2+iej)e  r 2^2^  1 1  +  e3^e2~^el^ 


iej  e2(  l+e3) 


■Zie^  (j lz-l{) 


r2U2)ll+  1 


From  Eq.  (155), 


i,rUi)  =  S,(ii)  +  t2ru,)  . 


or 


Zlr(l  1)  zlr(  12) 


z  1  r(  2 1 )  zlr(22) 
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7  4-  --  ■ 

(i+4> 

1(H) 

FI  1  ?  7 

det 

Cll+C12r2^1> 

-2ie  flu 


r2^2^11  ‘  ll 


det 


71 — m — : 

cn+cnr2<V 


-2ie.\^  (jL-£.) 

"(e2+ie1)e  rZ^  2>  1 1 e2  “ie  1  > 


det 


fi — m — : 

Cll  +  Cl2r2^l) 


-2ie.N fc(L-L) 

-(e^+ie^e  r2^  2^  1  l+^e2~iel^ 


Zl(22)  + 


ieje^l  +  e^)  /  -2ie  j'/u  i.^) 


'#  #  # 

cn+cur2«i) 


det 


r2^2^11  + 


This  relationship  gives  the  expression  for  each  element  of  ^rUj)  >  since 

two  matrices  are  equal  to  each  other  if  and  only  if  each  corresponding 
element  of  these  two  matrices  is  identical.  Thus, 


(i+4) 


'lr(  1 1  )-zi(  11)  rl  1  1 


det 


Cll+C12r2(fl) 
2 


-2iejN/w  (jf  2"Jf  j) 


Wnl 


ie1e2(  1  +e3) 

Zlr(22)=zl(22)+  F|  |  | 


det 


cn+ci2r2«i> 


-Zie^  (£2-£  j) 


r2^2>ll+1 


'lr(12)=Zlr(2l) 


det 


rj - # — i - : 

cn+ci2r2<*i> 


-Zie.\H2  (£2-f ,) 

-(e^+ie  j  )e  2^11  +  ^e2"iel^ 


(161) 

) 


Similar  to  the  way  Eqs.  (158)  and  (159)  are  solved,  Eq.  (154)  can  be 
expanded  in  complete  matrix  form, as  follows: 
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# 

rir<*l>  = 

e2e3"ie2zlr(21)+ie3Zlr(22)  '  1+ie2zlr(  1  l)'le3Zlr(  12) 

X 

e1+ie1e3zlr(21)+zlr(22)  -le3"iele3Zlr(ll)"Zlr(12) 

"ie3+ie1e3zlr(11)"zlr(12)  1+ie2zlr(ll)+ie3Zlr(12) 

X 

e1+ie1e3zir(21)'zlr(22)  "e2e3+ie2Zlr(21)+ie3Zlr(22) 

Hence,  the  quasi-matching  condition  =  0  gives  Eq.  (62)  in  the 

text  where  zlr^11j.  z  i  r(  1 2)  *  Zlr(2l)’  Zlr(22)  can  be  calculated  from 
Eq.  (161). 
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APPENDIX  J 


MATHEMATICAL  DERIVATIONS  OF  THE  IMPEDANCE  OF  A  FLOATING 
RING  DAMPER  ASSEMBLY  AS  THE  INTERIOR  SUPPORT 

With  reference  to  Figure  18,  a  free-body  diagram  can  be  drawn  for  an 
infinitesimal  mass  element  at  x  =  ,  as  shown  in  Figure  25.  When  the 

equilibrium  condition  at  point  (I)  is  considered, 

K  Y  =  C(Y,  -  Y  J  +  K,  (Y.  -  Y  )  , 

c  c  '  bt  ct;  b'  b  d 

of  which  the  Laplace  transform  is 


KCYC  =  sC(Yk  -  Y,)  +  Kk(Yk  -  Y J 


b'  b 


or 


~  sC+K,  ~ 
Y  = - Y. 

C  _  /■«  I  t  t  r  D 


1 0+1^+!^ 


(162) 


When  the  equilibrium  condition  at  point  (z)  is  considered, 


Ka(Y1  -  YJ  =  MKYKft  +  C(Y 


b  btt 


bt  •  Yct)  +  VYb  -  Yc>  • 


of  which  the  Laplace  transform  is 


<v  fsj 


rsj  rsj 


Ka<Yl  •  Yb>  =  8  MbYb  +  sC<Yb  -  Yc>  +  VYb  *  Yc> 


If  Eq.  (162)  is  substituted  into  the  above  expression, 


Yb  = 


Ka[  sC+(Kb+Kc)] 


s3MbC+.2Mb(Kb+Kc)+sC(Ka+Kc)+(KaKb+KbK<.+KcKa) 


Y,  .  (163) 


When  the  equilibrium  condition  at  point  (J)  is  considered, 

M  Y_+K  (Y.-YJ 
a  ltt  a'  1  V 


of  which  the  Laplace  transform  is 
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T(b  0) 


Figure  25.  Fret  *  tody  Diagram  for  Forces  Existing 
at  x  of  a  Floating  Ring 

Damper  Assembly. 
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rss 


■Rtfj)  = 


(S  Ma+K,)Y1-K4Yb 


If  Eq.  (163)  is  substituted  into  the  above  expression, 


»5MaMbC+.4MaMb<KbtKc)t*!lMa(VWa]C 


■R(ij)  = 


s3MbC+s2Mb(Kb+Kc)+s(Ka+Kc)C+(KaKb+KbKc+K  K  ) 


[(MatMb)Ka(VKc>+MaKbKclH'sKaKcC+KaKbKc  ~ 


(164) 


When  the  equilibrium  condition  at  point  (4)  is  considered, 
F({,-0)  -  F^+O)  =  -Rff,)  . 

If  Eq.  (164)  is  substituted  into  the  above  equation, 

F(«  j-0)  -  FUj+0) 

B4MA&,iMaMli(KbtKc)i.Va(VKc)tMbKalC 

s3MbC+s2Mb(Kb+Kc)+s(Ka+Kc)C+(KaKb+KbKc+KcKa) 

=  s  ' 


+s[(Ma+Mb)Ka(Kb+Kc)+MaKbKc]  i  KaKfaKc 


Y(C  j) 
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If  the  above  expression  is  compared  with 


F(£,-0)  -  F({j+0)  =  .ZjUjJYU,)  , 

it  can  be  concluded  that  the  support  impedance  at  x  =  f  is 

# 

s4MaMbC+s3MaMb(Kb+Kc)+s2[  Ma(Ka+Kc)+MbKj  C 
s3MbC+82Mb(Kb+Kc)+8(Ka+Kc)C+(KaKb+KbKc+KcKa) 


0 

.[(Ma+Mb)Ka(Kb+Kc)+MaKbKc]+KaKcC+  i  K^K,. 

'  —  • . —  . . ~  o 


0  J 

(165) 
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